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ABSTRACT 

(Distribution  Limitation  Statement  A) 

Research  is  described  on  a  pulsed  chemical  laser  operating  on 
the  reaction  between  hydrogen  and  fluorine  at  pressures  above  the  second 
explosion  limit.  A  laser  energy  density  of  80  joules/liter-atm  was  demon¬ 
strated  for  a  mixture  of  10%  F2/l0%  H2/80%  He  at  1.  i  atm  total  pressure. 
This  corresponds  to  a  chemical  efficiency  of  8%  and  an  overall  electrical 
efficiency  of  1.  3%  for  our  particular  photolysis  initiation.  The  highest 
vibrational  level  observed  to  give  laser  action  is  v  =  6,  and  it  is  concluded 
that  rapid  deactivation  of  the  upper  vibrational  levels  is  occurring.  It  was 
also  experimentally  determined  that  a  high  energy  electron  beam  could 
dissociate  F2with  an  expenditure  of  about  12  eV  of  energy  per  F-atom 
produced.  When  this  number  is  coupled  with  the  initial  F-atom  concentration 
required  for  good  laser  performance,  overall  electrical  efficiencies  of 
greater  than  100%  are  predicted. 

Attempts  have  also  been  made  to  develop  (1)  fast  yet  accurate 
computer  programs  which  contain  the  essential  physics  for  the  calculation 
of  the  vibration-relaxation  cross  section  in  an  atom  molecule  collision,  and 
(2)  theoretical  model  of  the  propagation  of  a  collimated  pulsed  laser  beam  in 
an  absorbing  medium.  The  energy  absorbed  by  the  medium  causes  heating, 
leading  to  an  index-of-refraction  change  which  causes  self-defocusing. 
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SECTION  I 


PULSED  HF  CHEMICAL  LASER 


1.  I  INTRODUCTION 


,  „  the  f  7'St  demonstration  by  Pimentel,  (0  chemical  lasers 

have  offered  the  promise  of  being  the  ultimate  laser  device.  Because  they 
use  the  chemical  ener<,y  of  reaction  to  produce  the  population  inversion, 
high  specific  power  is  possible  with  chemical  laser  devices.  By  employing 
a  chain  reaction,  hi«h  efficiency  is  also  possible.  Finally,  bv  pro-  ** 
ducing  laser  action  in  species  such  as  OF  or  HC1,  transitions  within  the 
region  ol  the  4  /i  atmospheric  window  are  possible,  thus  assuring  good 
atmospheric  propagation.  All  of  the  above  advantages,  while  significant 
are  accompanied  by  various  disadvantages  compared  to  CO?  or  CO  lasers 
For  example,  the  most  promising  chemical  lasers  utilize  halogen  or  halogen 
containing  compounds  with  all  the  attendant  handling  problems,  corrosive¬ 
ness,  and  toxicity  of  such  species.  Chemical  lasers  in  general  are  high 
gain  systems  and  undoubtedly  this  fact  will  be  a  major  technological  pro¬ 
blem  in  the  development  of  large  size,  high  power  devices  for  military 
applications.  Finally,  and  probably  most  important,  chemical  lasers  due 
to  their  very  nature,  necessitate  more  kinetic  information  and  detail  than 
™0.' St°th°*  laSfr  syste™s.  Frequently,  it  is  the  lack  of  these  scientific 
data  that  has  slowed  the  development  of  chemical  lasers  and  forced  this 
field  of  research  to  be  fairly  empirical  in  practice. 

ol  a  .u  ■  T?eSe,  disadv-itages  of  chemical  laser  systems  have  undoubtedly 
slowed  their  development  and  enthusiastic  support,  but  recent  advances 
appear  to  have  overcome  some  of  these  factors.  It  is  the  purpose  of  this 
report  to  describe  the  results  of  this  research  which  has  advanced  the 
development  of  chemical  lasers  for  military  applications. 


The  work  described  herein  is  the  second  year  of  a  two-year 
program  to  study  the  pulsed  H£/F2  laser.  The  mo.ivition  for  this  work 
was  the  belief  that  an  efficient,  high  energy  laser  could  be  developed  from 
the  rapid  H2 /F^  reaction  (this  is  often  referred  to  as  an  explosion  laser). 
The  Soviets  were  the  first  to  demonstrate  that  mixtures  of  H2  and  F2 
could  be  initiated  rapidly  and  made  to  lase.  It  is  intriguing  to  realize  that 
after  this  early  publication,  no  further  reports  on  this  type  of  laser  have 
been  noted  from  the  Soviet  Union. 


The  rationale  for  this  type  of  laser  is  straightforward.  There 
is  considerable  exothe rmicity  in  the  reaction  of  H2  with  Fz  to  form  HF  - 
about  6  eV/mole  of  reactants.  It  is  known  from  the  work  of  Poly  a  ni  <3  ) 
that  more  than  50%  of  this  exothermicity  ends  up  in  vibrational/rotational 
excitation  of  the  product  HF ,  and,  therefore,  in  principle  is  available  for 


extraction  as  laser  energy.  Furthermore,  the  reaction  mechanism  involves 
a  propagating  chain;  the  production  of  a  few  F-or  H-aioms  can  cause  the 
entire  mixture  to  react  to  completion.  A  chain  reaction  is  important  if  one 
is  to  obtain  high  overall  electrical  efficiency,  e.  g.  ,  greater  than  100%  conver¬ 
sion  of  electrical  energy  into  laser  photons.  By  replacing  H2  with  D2,  lasing 
on  DF  can  be  obtained  in  the  spectral  region  of  4  ji,  a  good  atmospheric  win¬ 
dow.  It  is  known  that  H2  and  can  react  spontaneously  (i.  e.  ,  explosively) 
upon  mixing.  There  are  regions  of  pressure  and  temperature  where  "stable" 
gas  mixtures  can  be  prepared.  In  our  investigations  we  have  chosen  to 
work  at  high  pressures,  above  the  second  explosion  limit.  This  requires 
that  the  gases  be  premixed  before  initiation  of  reaction  in  order  to  avoid 
a  mixing  limited  situation.  However,  the  advantage  to  be  gained  by  high 
prc-ssur"  operation  is  a  great  savings  in  overall  system  weight  due  to  the 
minimal  pumping  requirements.  Hence,  the  final  objective  of  this  research 
would  be  a  chemical  laser  device  opeiating  on  D2/F2  gas  premixed  at  about 
atmospheric  pressure  and  volumetrically  initiated  for  scalability.  For 
high  power  applications,  such  a  device  would  have  to  be  repetitively  pulsed, 
and  hence  would  involve  flow  and  mixing  of  the  gases  upstream  of  the  laser 
cavity. 

During  the  first  year's  research  effort  the  major  accomplish¬ 
ments  were  (a)  demonstration  of  the  ability  to  produce  stable  mixture  of 
H^  +  F^  +  various  diluents  at  room  temperature  and  pressure  around  one 
atmosphere,  and  (b)  the  demonstration  of  laser  action  in  some  of  these 
mixtures  initiated  by  flash  photolysis  of  the  F^.  The  overall  laser  energy 
and  efficiency  obtained  were  small.  It  was  found  that  more  reproducible 
stable  mixtures  could  be  produced  in  a  slow  flow,  mixing  device.  The  data 
obtained  in  the  first  year's  study  indicated  that  increased  laser  efficiency 
and  energy  should  be  obtainable  by  increasing  the  photolysis  of  the  F2  to 
produce  more  F-atoms  in  order  to  speed  up  the  chain  reaction  and  by 
increasing  the  optical  cavity  length  to  increase  the  rate  of  stimulated 
emission. 

Thus,  the  second  year's  program  started  with  the  objective  of 
improving  the  overall  energy  density  and  efficiency  of  the  pulsed,  atmo¬ 
spheric  pressure  H2/F2  laser.  A  larger  laser  cavity  device  was  con¬ 
structed  using  coaxial  flash  lamps  for  photolysis  with  improved  optical 
coupling  to  the  cavity.  With  this  device  a  parametric  study  of  laser  per¬ 
formance  was  investigated.  To  understand  and  correlate  the  trend  of  the 
data  and  to  provide  overall  insight  into  the  laser  phenomenology  a  theoret¬ 
ical  modeling  effort  was  also  performed  and  compared  where  possible  to 
the  experimental  data.  In  support  of  the  modeling  effort,  a  theoretical 
calculation  of  important  kinetic  cross  sections  was  undertaken.  Besides 
the  study  of  the  laser  itself,  an  experimental  investigation  was  also 
carried  ort  of  the  efficiency  of  a  high  energy  electron  beam  as  a  possible 
volumetric  initiation  scheme. 

The  major  conclusion  of  this  second  year's  research  program 
are  the  following: 
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A  laser  energy  density  of  80  joules/liter-atm  ’/as  demonstrated 
for  a  mixture  of  10%  F2/10%  H2/80 %  He  at  1.  1  atm  total  pressure.  This 
corresponds  to  a  chemical  efficiency  of  8%  and  an  overall  electrical  effi¬ 
ciency  of  1.  3%  for  our  particular  photolysis  initiation. 

The  main  features  of  the  laser  performance  appear  to  be  under¬ 
standable  and  predictable  via  our  theoretical  model.  However,  the  agree¬ 
ment  with  data  is  only  qualitative  and  certain  important  effects,  such  as  the 
influence  of  rotational  non-equilibration  a. id  optical  output  coupling,  remain 
to  be  understood. 

It  was  experimentally  determined  that  a  high  energy  electron 
beam  could  dissociate  F2  with  an  expenditure  of  about  12  eV  of  energy  per 
F-atom  produced.  When  this  number  is  coupled  with  the  initial  F-atom 
concentration  required  for  good  laser  performance,  overall  electrical 
efficiencies  of  greater  than  100%  are  predicted.  These  resu  ts  are  very 
significant  because  they  indicate  that  a  high  efficiency,  scalable  chemical 
laser  operating  at  atmospheric  pressure  is  close  to  reality. 

The  following  sections  present  the  significant  results  of  this 
investigation  of  the  pulsed,  atmospheric  pressure,  HF  laser.  Section  1.  2 
describes  the  experimental  details  including  the  various  diagnostics  used 
to  study  laser  performance  and  summarizes  the  experimental  data  obtained 
on  the  laser  device.  Section  1.  3  describes  the  analytic  model  on  laser 
performance  and  comparisons  with  data.  Section  1.4  summarizes  the 
conclusions  on  laser  performance  and  understanding.  Section  1.  5  describes 
a  theoretical  approach  to  calculate  vibrational  energy  exchange  for  species 
such  as  HF.  Section  1.  6  presents  the  experimental  data  on  laser  initiation 
techniques  including  a  pin  discharge  and  a  high  energy  electron  beam. 
Section  1.  7  presents  some  overall  comments  on  laser  scalability  and  future 
device  development. 
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1.2  EXPERIMENTAL 


1 .  2 .  I  Desc rip ,~ion  of  Laser  Experiments 

1.2.1  I  Experimental  Setup  and  Laser  Operations 

The  Avco  Everett  Research  Laboratory,  Inc.  ,  experiments 
with  the  pulsed  H2/F2  chemical  laser  had  shown  (5)  that  the  require¬ 
ments  for  improved  chemical  efficiency  were,  (i)  a  longer  optical 
path  length  to  increase  stimulated  emission,  and  (L)  a  higher  concentration 
of  fluorine  atoms  produced  by  the  flash  lamps  to  trigger  the  reaction. 

Thus,  under  the  present  contract  the  apparatus  was  modified  to  increase 
the  optical  path  length.  In  an  attempt  to  increase  the  fluorine  atom  con¬ 
centration,  coaxial  flash  lamps  were  also  usjd.  Since  these  coaxial  lamps 
surround  the  gas  that  is  to  be  photolysed,  they  offer  improved  optical 
coupling  to  the  cavity,  whereas  the  geometry  of  the  linear  lamps  used  pre¬ 
viously  was  such  that  only  a  small  fraction  of  their  output  was  coupled  into 
th’  gas.  Due  to  their  lower  inductance,  the  coaxial  lamps  also  have  a  much 
shorter  pulse  length.  It  appears  from  the  experimental  results  that  a  higher 
fluorine  atom  concentration  is  indeed  obt  ained;  the  efficiency  achieved  is 
about  two  times  higher  than  before,  and  larger  output  energies  are  also  ob¬ 
served. 


A  diagram  of  the  laser  apparatus  is  given  in  Tig.  1.  Fluorine 
and  hydrogen  are  each  premixed  with  diluent  helium  in  separate  storage 
tanks  at  pressures  above  ten  and  five  atmospheres  respectively.  During  a 
run,  fast  acting,  electrically  operated  valves  open  to  admit  these  gases  to 
a  mixing  manifold  which  distributes  the  gases  in  the  mixing  region  in  alter¬ 
nate  sheets  2  mm  thick.  The  cross  section  of  this  mixing  region  is  12  x 
10  cm.  The  storage  tanks,  lines  and  valves  are  of  stainless  steel,  the 
mixing  region  is  made  of  aluminum  and  Teflon  sheets.  The  mixed  gases 
pass  into  two  coaxial  Xe  flash  lamps  and  then  exhaust  into  a  4-liter 
dump  tank.  A  Kistler  pressure  gauge  mounted  upstream  of  the  optical 
cavity  is  used  to  confirm  that  the  xplosion  originated  in  the  cavity 
(see  Fig.  1).  Figure  2  shows  a  typical  oscillogram  of  the  pressure 
gauge  signal.  Note  that  the  abrupt  change  in  pressure,  indicating 
reaction,  occurs  coincident  with  the  flash  lamp  firing.  On  completion 
of  a  run,  the  fast  acting  valves  close  and  other  valves  open  to  admit 
helium  to  the  manifolds  and  to  exhaust  the  dump  tank.  The  system 
is  thus  flushed  automatically  at  the  end  of  each  run. 

The  optical  cavity  is  composed  of  an  8-m  radius  confocal  copper 
mirror  and  various  other  elements  such  as  sapphire,  salt,  germanium,  or 
silicone  plate  reflectors.  Both  plate  and  mirror  are  installed  in  aluminum 
housings  at  the  ends  of  the  flash  lamps,  and  are  protected  from  the  corro¬ 
sive  gases  by  a  helium  bleed  flow.  Figure  3  shows  the  diagram  of  the  porous 
wall  bleed  system  used  in  the  experiment.  The  flash  lamps  have  an  internal 
diameter  of  17  mm,  and  an  active  length  of  25  cm;  the  active  length  of  the 
cavity  is  therefore  about  50  cm.  The  lamps  have  a  duraction  (full  width  at 
half  height)  of  about  l/l  sec,  and  are  each  energized  by  a  1.5  /if  capacitor 
at  20  kV,  giving  300  joules  per  lamp.  Figure  4  snows  a  typical  flash  signal 
monitored  by  an  IP28  photomultiplier. 
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Fig.  1  Schematic  of  Flash  Photolysis  Laser  Experiment 
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Diagram  of  the  Bleed  System  Used  to  Protect  the  Mirrors 
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Fig.  4  The  Pulse  Shape  of  the  UV  Signal  Emitted  by  the  Flash  Lamp 
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Measurements  of  the  laser  output  energy  were  made  with  a 
calorimeter.  The  details  of  construction  of  this  calorimeter  are  described 
later  The  temporal  variation  of  the  laser  duration  was  monitored  with  an 
InSb  PEM  detector  viewing  scattered  radiation  from  the  calorimeter  surface. 
For  spectral  analysis,  an  0.  5  meter  Jarrell-Ash  monochromator  fitted  with 
a  gold-cioped  germanium  detector  was  used.  The  signals  from  the  detectors 
are  amplified,  displayed  on  an  oscilloscope  and  photographed. 

Matheson  prepurified  grade  H2  (99.95%),  high  purity  grade  He 
(99.95%),  and  high  purity  grade  O2  (99.95%)  were  used  without  further 
purification.  The  reason  for  using  O2  in  the  mixtures  is  given  in  Section 

1 . 2.  2.  1  In  order  to  remove  any  TF  contamination  (  - 1%)  from  the  F 2  tank, 
F2  gas  (98%),  supplied  from  Matheson,  was  passed  through  an  NaF2 
adsorber  prior  to  the  preparation  of  F2/He  mixtures.  It  was  found  that 
more  than  90%  of  the  HF  in  F2  could  be  removed  by  this  technique.  The 
temperature  of  both  the  laser  cavity  and  mixing  bottles  was  room  tempera¬ 
ture  in  all  experiments. 

1.2.  1.2  Diagnostic  Experiments 

1 . 2 .  1 . 2 .  1  Measurement  of  the  Percent  of  F?  Dissociation 
Produced  by  the  Flash  Lamp 

Since  the  rate  of  the  H2/F2  reaction  depends  upon  the  value  of 
the  concentration  of  F  atoms  [  F]  initially  produced,  a  precise  measure¬ 
ment  of  [  F]  or  %  F2  dissociation  is  very  important.  Two  efforts  were 
made  to  measure  the  initial  fluorine  atom  conce  itration.  First,  an  ultra¬ 
violet  He-Cd  laser  (Coherence  Radiation  Corp.  )  operating  at  3250  A  was 
*  used  to  measure  the  fluorine  molecule  concentration  with  the  object  of 

determining  the  percentage  dissociation  of  F2  caused  by  the  flash  lamps. 

This  experiment  was  unsuccessful  due  to  (i)  the  small  initial  percentage 
dissociation  (<1%)  and  hence  the  small  change  in  3250  R  absorption  and 
(ii)  the  large  inherent  fluctuation  in  intensity  of  the  laser  signal  (>10%). 
Clearly  a  much  more  stable  He-Cd  laser  is  needed  for  this  measurement. 
Due  to  lack  of  resources,  no  effort  was  expended  in  attempting  to  stabilize 
the  He-Cd  laser. 

The  second  technique  was  to  react  the  F  atoms  initially  pro¬ 
duced  by  the  flash  with  HCf  to  form  HF 

F  +  HCf  -*  HF  +  Cf  (1) 

and  measure  the  concentration  of  HF  with  a  probe  HF  laser.  In  this  tech¬ 
nique  no  H2  is  used  in  the  gas  mixture,  so  that  in  principle  no  chain  reaction 
occurs.  The  quantity  of  HF  produced  is  then  exactly  the  number  of  Huorine 
atoms  present  initially.  The  (  HF]  can  be  measured  after  sufficient  delay 
for  vibrational  relaxation,  by  absorption  of  an  HF  laser  pulse. 

Figure  5  shows  the  experimental  arrangement  used.  The  HF 
IP4  transition  was  generated  by  an  AVCO  C950  N2  laser  operated  on  a 
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Fig.  5 


Diagram  of  Apparatus  Used  for  Determination  of  Concentration 
of  HF  Produced  by  Reaction  of  F  +  HCf  — HF  +  C£ 
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mixture  of  SFfc  and  H2.  The  BaF2  beam  splitter  reflects  about  10%  of  the 
IP4  line  into  Ge-Au  detector  #1  and  transmits  90%  through  the  test  cell  and 
thereafter  into  the  monochromator  and  Ge-Au  detector  if 2.  The  test  cell 
is  the  actual  laser  cavity  with  two  BaF2  windows.  The  experimental  se¬ 
quence  is  identical  to  that  of  the  laser  experiment  except  that  H2  is  replaced 
by  HC£  in  one  mixing  tank.  A  time  delay  generator  was  used  to  produce  a 
*  desired  delay  for  firing  the  probe  laser  after  the  flash  lamp  initiation.  Thus, 

detector  ft  l  monitors  the  value  of  IQ  (initial  signal)  and  detector  §2  monitors 
the  intensity  I  with  absorption  simultaneously.  For  mixtures  containing 
0.  1  atm  F2  an  optical  path  length  of  about  50  cm  is  used.  For  higher  F2 
pressures,  an  optical  path  of  5  cm  was  used.  An  averaged  absorption  co¬ 
efficient  for  the  IP4  HF  line  determined  previously(  -  )  was  used;  this  value 
is  3  cm-  c-rr"1.  The  results  of  this  experiment  are  discussed  in  Section 
1. 2.  2.  5. 

1 . 2 .  1 . 2 .  2  Measurement  of  the  Rale  of  the  H 2  +  F 2  Reaction 
by  Monitoring  the  Transient  F?  Concentration 

There  has  been  a  controversy  in  the  literature  (6  )  concerning 
the  possibility  of  chain  branching  reactions  in  the  H2/F2  explosion.  Un¬ 
fortunately,  no  detailed  information  on  the  rate  constants  of  the  branching 
reaction(s)  is  known.  The  fact  that  a  second  explosion  limit  exists  in  the 
H2/F2  system  indicates  that  chain  branching  does  participate  in  the  re¬ 
action.  Russian  scientists  O  )  have  claimed  that  this  chain  branching  oc¬ 
curs  through  the  reaction 

HF  (v>4)  +  F2  -*  HF  +  2F  (2) 

<  which  indicates  that  HF  vibrationally  excited  in  states  greater  than  v  =  3 

can  participate  in  dissociating  F£.  It  is  not  obvious  whether  such  a  re¬ 
action  should  increase  or  decrease  laser  efficiency:  by  speeding  up  the 
reaction,  it  should  increase  the  efficiency,  however,  the  increased  re¬ 
action  rate  is  accomplished  by  removing  vibrationally  excited  HF,  which 
should  decrease  efficiency.  However,  if  branching  did  occur  with  a  rea¬ 
sonably  fast  rate,  the  H2/F2  reaction  could  proceed  much  more  rapidly 
than  without  branching.  The  temperature  of  the  reacting  gas  would  also 
increase  more  rapidly,  and,  hence,  the  laser  duration  would  be  shorter. 

In  order  to  experimentally  estimate  the  overall  rate  of  reaction  in  our 
H2/F2  system,  the  transient  concentration  of  F2  was  measured. 

Figure  6  shows  the  experimental  set  up.  An  AVCO  C950  N2 
laser,  producing  2-kW  pulses  at  3371  R  was  used  as  a  stable  light  source. 
The  reflected  3371  A  beam  from  the  front  surface  of  the  back  sapphire 
window  was  monitored  by  the  1P28  photomultiplier  equipped  with  a  100  R 
bandpass  UV  transmission  filter  centered  at  3371  J?.  The  detector  is 
located  11  meters  away  from  the  flash  lamp.  To  avoid  any  flash  lamp 
light  scattering  into  the  detector,  a  10-meter  long  black  box  with  small 
(5  mm  diameter)  entrance  and  exit  holes  was  placed  between  the  photo¬ 
multiplier  and  the  flash  lamp.  The  overall  response  time  of  the  detector 
and  amplifier  are  less  than  0.2jisec.  A  time  delay  box  was  used  to  control 
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Diagram  of  Apparatus  Used  for  Determination  of  the  Transient 
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Fig.  6 


the  firing  tine  of  UV  laser  after  the  flash  lamp  initiation.  The  transient 
percentage  absorption  of  the  3371  A  line  by  F2  molecules  in  the  laser 
cavity  is  repeatedly  measured  with  varying  time  delays  after  flash  initia¬ 
tion.  Based  on  these  data,  the  fraction  of  F2  reacted  at  various  times 
and  hence  the  rate  of  the  H2/F2  reaction  ie  calculated.  The  results  are 
given  in  Section  1.2.  2.  7. 

1.2.  1.2.3  Measurement  of  Laser  Output  Spectra 

For  a  better  understanding  of  the  detailed  reaction  kinetics  in¬ 
volved  in  the  H2/F2  system,  we  have  investigated  the  time -dependent 
spectra  of  our  HF  laser.  Figure  7  shows  the  optical  set  up  for  measuring 
both  total  laser  output  energy  and  laser  transitions  simultaneously.  An 
InSb  PEM  detector  was  used  to  monitor  HF  IR  laser  emission.  The  op¬ 
erating  time  constant  was  less  than  0.  ljisec.  Transition  wavelengths  were 
identified  on  a  0.5-meter  Jarrell-Ash  mon,?:hromator  (Model  82-000) 
equipped  with  a  4|j.-blaze  grating  and  calibrated  with  a  visible  He-Ne  laser. 
Laser  pulse  energies  were  measured  with  a  calibrated  AVCO  made  energy 
meter  (see  next  section).  Limited  experiments  on  D2/F2  laser  action 
were  also  performed.  The  results  are  given  in  Section  1. 2.  2.4  and  1.  2.  2.  5. 

1.2.  1.2.4  Construction  and  Calibration  of  an  Energy  Meter 

The  energy  meter  used  previously  was  a  commercial  product 
(TRG  model  100)  which  had  an  input  aperture  of  only  1-cm  diameter. 

If  the  output  beam  were  considerably  larger  than  this,  clearly  erroneous 
readings  would  be  obtained.  It  was  decided  that  a  large  size  energy  meter 
was  required  for  this  experiment,  placed  as  close  to  the  output  mirror  as 
possible.  Such  energy  meters  are  available  commercially,  but  only  with 
a  long  delivery  time,  and  at  quite  high  cost.  These  commercial  meters 
employ  a  polished  cone  to  collect  the  radiation,  and  although  the  manu¬ 
facturers  claim  that  they  are  equally  responsive  to  all  wavelengths,  it  was 
not  clear  that  this  claim  had  been  verified.  Moreover,  at  high  input 
energies,  a  cone  will  focus  the  radiation,  creating  an  air  breakdown  which 
prevents  the  meter  from  receiving  all  the  energy.  It  was  decided,  there¬ 
fore,  that  it  would  be  better  to  build  an  energy  meter. 

The  general  arrangement  of  the  resulting  meter  is  shown  in  the 
top  half  of  Fig.  8.  Radiation  is  absorbed  by  the  grooved  anodized  face  of 
a  thin  piece  of  aluminum  5 -cm  diameter.  On  the  back  of  the  aluminum  are 
attached  19  thermocouple  junctions,  and  19  junctions  are  also  attached  to 
a  similar  reference  slug  of  aluminum,  wired  in  series  as  indicated  in 
Fig.  8  to  increase  the  output  signal.  The  meter  is  supported  by  lucite 
stand-offs  inside  a  copper  shield  to  prevent  air  currents  from  affecting 
the  signal.  The  output  signal  is  fed  to  a  micro  voltmeter.  This  calori¬ 
meter  is  capable  of  detecting  a  millijoule  of  energy. 


It  was  calibrated  by  both  calculations  of  the  heat  rise  experi¬ 
enced  by  the  aluminum  on  absorbing  a  given  quantity  of  energy,  and,  in 
addition,  a  set  of  HF  laser  pulses  from  an  AVCO  C950  laser  operating  on 
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Diagram  of  Energy  Meter  and  the  Calibration  Tests  of  the 
Energy  Meter 


SF6  and  H2  was  directed  onto  the  meter  after  being  focused  down  and 
passed  through  an  aper'.ure  (see  lower  half  of  Fig.  8).  This  aperture  was 
smaller  than  the  >  ntrance  aperture  of  two  available  TRG  meters,  and  these 
were  also  used  tj  measure  the  energy  of  the  HF  pulses.  The  energy  of  a 
single  HF  laser  pulse  as  measured  by  each  instrument  is  given  in  Fig.  8, 
and  shows  good  agreement  for  all  3  meters. 

The  19  thermocouples  on  the  back  of  the  aluminum  slug  were 
spaced  so  they  each  served  an  approximately  equal  area.  In  this  case  it 
can  be  shown  that  the  signal  is  independent  of  whether  the  energy  is  ap¬ 
plied  at  a  single  point,  or  over  the  whole  surface  of  the  detector.  This 
was  verified  experimentally  by  moving  the  energy  meter  away  from  the 
aperture  so  that  the  pulse  spread  to  fill  the  entire  aluminum  slug.  Pro¬ 
vided  the  size  of  the  laser  spot  was  not  greater  than  the  meter  face,  the 
signal  was  identical  whether  the  meter  was  close  to  the  aperture,  with  the 
pulse  only  falling  on  a  small  region  of  the  meter,  or  far  away. 

Before  actually  building  the  energy  meter,  tests  were  made  of 
the  reflectivity  of  various  surfaces  at  various  wavelengths,  with  the  ob¬ 
jective  of  finding  a  surface  that  was  fully  absorbing  at  all  wavelengths.  It 
was  known  that  black  anodized  aluminum  was  an  excellent  absorber  at 
10.6/i,  so  a  flat  sample  of  this  was  chosen  as  the  reference  point.  Two 
other  surfaces  with  circular  grooves  were  made,  one  with  a  30  half-angle 
groove,  the  other  with  a  i5°  half-angle  groove,  and  both  were  black 
anodized.  These  two  surfaces,  plus  the  flat  surface,  were  measured  with 
a  reflecting  attachment  on  a  Perkin  Elmer  spectrophotometer  between  2 
and  11)1.  Unfortunately,  we  oossess  no  instrument  for  making  total 
reflectivity  measurements.  Nor,  despite  the  abundance  of  manufacturers 
of  optical  accessories  in  the  Boston  area,  was  it  easy  to  find  anyone  else 
who  had  one.  Finally,  one  was  located  at  EG&G,  and  this  firm  kindly 
agreed  to  run  the  samples  for  us.  Their  instrument,  however,  could  only 
measure  between  1  and  2.7  fi.  However,  the  traces  taken  showed  a  re¬ 
markable  independence  of  wavelength  within  their  measuring  range,  and  it 
was  felt  safe  to  extrapolate  this  result  to  longer  wavelengths.  The  results 
are  presented  in  Fig.  9.  Where  N.M.  appears,  it  means  no  measurement 
was  made.  The  measurements  on  the  flat  sample  showed  that  there  is  a 
considerable  diffuse  reflectivity  which  is  reduced  somewhat  by  the  30  half¬ 
angle  grooves,  and  considerably  reduced  by  the  15  half-angle^rooves . 

On  the  basis  of  these  tests,  the  black  anodized  surface  with  15  half-angle 
grooves  was  considered  to  be  equivalent  to  a  fully  absorbing  surface,  and 
was  used  as  the  collecting  surface  of  the  energy  meter  that  we  constructed. 

1.2.2  Results 


1 .  .  2 .  1  Gas  Mixture  Stabilization 

During  the  course  of  our  experimentation,  it  had  been  noticed 
that  different  results  were  obtained  with  different  bottles  of  fluorine  -- 
some  bottles  showing  spontaneous  reaction  on  mixing,  others  being  more 
stable.  Since  the  work  of  Levy  and  Copeland(®  )had  shown  that  oxygen 
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inhibits  the  spontaneous  reaction  ol  H2/F 2,  it  seemed  natural  to  suspect 
that  varying  amounts  of  oxygen  in  the  fluorine  tanks  were  causing  this  ef¬ 
fect.  The  manufacturer  claimed  a  maximum  oxygen  level  of  2%  in  the 
fluorine.  Consequently,  amounts  of  oxygen  up  to  2%  of  the  fluorine  were 
added  to  a  mixture  that  exhibited  spontaneous  reaction.  The  reaction  was 
monitored  by  the  Kistler  pressure  gauge  as  the  gases  flowed  through  the 
cavity.  The  results  were: 


i)  Additions  of  oxygen  to  the  hydrogen  mixing  tank  did  not 
affect  the  mixture  stability. 

ii)  Addition  of  oxygen  to  the  fluorine  mixing  tank  enhanced 
the  mixture  stability,  increasing  the  time  to  spontaneous 
explosion  on  mixing,  until  at  0.4%  of  the  total  mixture, 
the  gas  was  completely  stable,  i.  e.  ,  it  did  not  react  at 
all  unless  triggered. 

Next,  the  effect  of  oxygen  addition  on  laser  output  was  studied. 
Using  a  gas  mixture  composition  of  H2/F2/He  =  0.  08/0.  04 /0.  88,  various 
amounts  of  oxygen  were  added.  The  actual  technique  was  to  use  helium  to 
which  oxygen  had  been  admixed.  The  laser  energy  was  measured  using  a 
stable  semi-confocal  cavity  with  a  barium  fluoride  output  coupling  flat. 

The  results  are  shown  in  Fig.  10,  indicating  a  steep  rise  in  energy  as  the 
mixture  becomes  less  stable,  i.e.,  less  added  O2.  With  a  fairly  unstable 
mixture  an  energy  of  1. 8  joules  was  obtained.  However  above  0.  2%  of 
added  O2  the  mixture  is  very  stable,  and  the  output  energy  is  independent 
of  oxygen  concentration.  Since  the  reproducibility  is  the  most  important 
thing  in  any  systematic  study,  we  have  used  gas  mixtures  with  0.2  or  0.4% 
of  added  O2  for  all  our  subsequent  experiments. 


Possible  reasons  for  the  stabilization  of  H2/F2  mixture  by 
molecular  oxygen  may  be  postulated  as  follows: 


F2  +M 


"►  2F  +  M  (where  M  is  surface) 


(3) 


F  +  Oz  +  M' 


FOz  +  M'  (where  M1  could  also  be  a  surface)  (4) 


followed  by 


fo2+F 


-Fj  +°2 


For  steady-state  conditions,  the  following  relations  hold: 

HP  =  \  kl  J  FZ^  [  MJ  -  k2  [  F]  [  °21  [  M'  J  -  k3  (  F02]  [  F] 


=  0 


(5) 


(6) 


f 
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ENERGY  (JOULE) 


s .  s . 


(7) 


m 


i  kj  [F2]  (M) 
k2  [02]  [  M1  ]  +  K3  [F02] 


For  the  case  of  k2  [  O2]  [  M'  ]  «  k3  [  FO2) ,  [  F]  s .  s  _  does  not  depend  on  the 
amount  of  02  present.  But  for  the  case  of  k2  [  O]  [M'J  »k3[F02],  [  F]  s 
is  inversely  proportional  to  the  amount  of  O2  added.  As  02  is  added  to  the 
gas,  the  F-atom  concentration  decreases  to  a  value  where  it  can  no  longer 
initiate  the  chain.  Further  addition  of  02  does  not  affect  the  mixture  stabil¬ 
ity.  The  fact  that  laser  output  energy  decreases  as  02  is  increased  (  and 
the  mixture  becomes  noticeably  more  stable)  could  be  due  to  the  presence 
of  the  small  steady-state  concentration  of  F02  produced  by  Reaction  (4) 
which  could  act  as  a  chain  terminator  via  Reaction  (5).  02  could  also  affect 

laser  action  by  the  formation  of  H20  which  could  deactivate  vibrationally 
excited  HF.  However,  calculations  indicate  that  the  formation  of  H2O  is 
too  slow  under  our  conditions  to  have  an  effect  on  laser  output.  This  is 
consistent  with  the  observation  that  further  addition  of  O2,  beyond  that 
necessary  for  stabilization,  does  not  appear  to  reduce  laser  output  energy 
(see  Fig.  10). 

Unfortunately,  none  of  the  above-mentioned  rate  constants  by 
which  a  more  quantitative  discussion  could  be  made  are  known.  At  the 
present  time,  we  can  only  speculate  that  this  may  be  the  reaction  scheme 
involved  in  the  stability  of  H2/F2  mixtures.  However,  it  is  of  considerable 
interest  to  note  that  additional  laser  energy  might  be  available  if  stability 
could  be  achieved  by  other  than  addition  of  O2,  e.  g.  ,  by  cooling  the  gas. 

1 . 2 .  2 .  2  Output  Energy  Measurements 

Experience  with  the  earlier  version  of  the  pulsed,  flowing  laser 
had  shown  that  nitrogen  is  inferior  to  helium  as  a  diluent.  This  is  probably 
due  to  the  fact  that  helium  is  less  effective  than  nitrogen  in  deactivating 
vibrationally  excited  hydrogen  fluoride/9)  We  have  therefore  used  only 
helium  as  the  diluent  in  our  experiments.  The  optimum  ratio  found  pre¬ 
viously,  namely  H2/F2/He  =  0.04/0.04/0.92  was  normally  used  as  a 
standard  mixture.  Runs  were  also  made  with  various  other  mixture  ratios, 
with  deuterium  replacing  hydrogen,  and  with  other  mixed  halogen  fluoride 
compounds  in  place  of  fluorine.  HF  and  DF  laser  radiation  was  systemat¬ 
ically  studied  in  several  different  ways: 

1.  Varying  the  flash  lamp  energy, 

2.  Varying  the  F2  concentration  with  the  H2  concentration 
fixed, 

3.  Changing  the  initial  F-atom  concentration  by  introducing 
additive  s, 

4.  Changing  both  the  H2  and  F2  concentrations  simultaneously, 


1 


* 
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5.  Varying  the  output  coupling,  and 

6.  D2/F2/He  mixture  results. 

The  results  can  be  summarized  as  follows: 

1.2.  2.  2.1  The  Effect  of  Varying  Flash  Lamp  Energy 

The  laser  pulse  energy  was  measured  at  various  flash  lamp 
energies  (and  hence  at  various  initial  F-atom  concentrations)  for  a  gas 
mixture  of  H2/F2/He  =  0.  12/0.  04/0.  84  at  a  total  pressure  of  1 .  1  atm  and 
90%  output  coupling  obtained  with  a  total  reflecting  copper  mirror  at  one 
end  and  a  BaF2  plate  at  the  other.  The  H2  enriched  mixture  was  used  to 
assure  that  all  the  F-atoms  produced  by  the  flash  lamp  would  react  with 
the  excess  H2  molecules.  In  this  way  a  direct  correlation  between  the 
lamp  energy  and  laser  output  energy  can  be  obtained.  In  Fig.  11  the 
measured  laser  pulse  energy  is  plotted  against  the  energy  of  the  flash  lamp 
capacitor  bank.  As  one  can  see,  up  to  the  maximum  flash  energies  avail¬ 
able  with  this  apparatus  (700  J)  the  laser  energy  is  still  increasing  with 
flash  energy  or  equivalently  initial  F-atom  concentration. 

j.2.2.2.2  The  Effect  of  Varying  F?  Concentration 

Another  means  of  varying  the  F-atom  concentration  is  to 
change  the  F2  concentration  at  fixed  flash  energy.  These  measurements 
were  done  at  fixed  H2  concentration  and  variable  F2  ratios.  The  experi¬ 
mental  conditions  are  exactly  the  same  as  described  above,  namely  a  90% 
output  coupling  and  a  total  pressure  of  1 .  1  atm.  Figure  12  shows  the 
dependence  of  output  energy  as  a  function  of  F2  concentrations.  An  increase 
t  jn  Iaser  energy  as  F2  concentration  increases  is  observed.  The  slope 

gives  a  value  of  40  joules  of  laser  photons/atm  of  F2  used  in  an  active  vol¬ 
ume  of  0.  1  liter  at  the  described  experimental  conditions. 

1 . 2 .  2 .  2 .  3  The  Effect  of  Increasing  F-atom  Concentration  by 
Introducing  Additives  such  as  MoF^ 

Molecular  MoFfc  has  been  noted  to  have  a  large  absorption  co¬ 
efficient  in  the  UV  spectral  region  where  the  xenon  flash  lamps  emit. 

Hess  has  recently  reported^0)  that  an  enhancement  in  HF  laser  output 
energy  was  observed  as  MoF^  was  added  to  his  H2/F2  system.  It  is 
believed  that  this  enhancement  was  due  to  the  additional  initial  F-atoms 
generated  by  the  M0F6  flash  photolysis.  Since  both  of  our  experiments 
above  indicate  that  the  output  energy  will  increase  with  [  F]j,  a  trace  of 
MoF6  (0.  1%)  was  added  to  the  0.  12/0.  04 /0.  84  mixture.  Runs  were  made 
at  various  values  of  capacitor  flash  lamps  energy.  The  result  is  shown  in 
Fig.  13.  Again,  it  is  found  that  the  laser  output  increases  with  the  flash 
lamp  energy,  but  less  energy  is  observed  with  MoF(,  than  without  it.  At 
low  flash  lamp  energy,  the  opposite  result  is  obtained.  The  addition  of 
MoFfc  should  increase  the  [F]i(  as  should  an  increase  in  the  flash  lamp 
energy.  If  the  extra  fluorine  atoms  produced  by  MoFfc  were  actually 
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Fig.  12  The  Effect  of  Changing  F2  Concentration  to  the  Laser  Output 
Energy  Observed  in  the  H-^/F-j/He  Mixture 
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The  Effect  of  Changing  MoF^  Concentration  to  the  Laser  Output 
Energy  Observed  in  the  H?/F-,/He  Mixtures 


deactivating  the  HF,  one  would  expect  the  laser  energy  to  fall  at  high  flash 
lamp  energy.  It  is  not  likely  that  MoF^  has  extremely  high  quenching  ef¬ 
ficiency  for  vibrationally  excited  HF  molecules. 


The  interpretation  of  this  result  is  not  very  clear.  One  possi¬ 
ble  explanation  is  that  some  kind  of  chain  terminating  reaction  is  introduced 
as  MoF^  is  added  to  the  system.  For  example,  radicals  like  M0F5  gen¬ 
erated  from  flash  photolysis  of  MoF^ 

fla  sh 

MoF^  »  MoF^  +  F  (8) 

could  serve  as  a  chain  terminating  reactant 

MoF5  +  F  -  MoF4  +  F2  (9) 


This  could  result  in  a  lower  laser  energy  although  [  F)  j  might  be  higher. 
Another  possibility  is  that  M0F5  could  act  as  an  effective  quencher  for 
rotational  deactivation  of  HF.  Process  like 


HF  (v=n,  J=jp )  +  MoF,  (v.=o) 

o  1 

-  HF  (van,  J=f  -1  or  (- 2 )  +  MoF,  (v.=l) 

0  1 


(10) 


could  be  quite  fast  due  to  the  close  resonant  rotation-to-vibrational  energy 
transfer  between  M0F5  and  HF.  As  will  be  apparent  later,  rotational 
equilibration  effects  can  be  important  in  the  H2/F2  laser. 

If  rotational  relaxation  of  HF  by  M0F5  is  responsible  for  the 
observed  decrement  in  laser  output  energy,  a  similar  effect  should  also  be 
expected  for  H2/F2  mixtures  with  SF^  additive.  Since  both  SF^  and  MoF^ 
have  low  lying  vibrational  energy  levels,  both  should  be  capable  of  remov¬ 
ing  rotational  energy  of  HF  by  the  process  of  near  resonant  rotation-to- 
vibration  energy  transfer.  However,  SF^  is  not  dissociated  via  the  flash 
and  should  not  act  as  a  chain  terminator.  Figure  14  shows  the  effect  of 
SF^  additive  to  the  H2/F2  laser  output  energy.  A  drop  in  output  energy 
was  observed  as  moie  and  more  SF^  was  added.  Molecular  SF^  has  long 
been  noted  to  be  rather  inefficient  in  the  vibrational  deactivation  of  HF.(H^) 
This  drop  in  output  energy  can  only  be  explained  by  a  possible  R -*  V 
process  which  changes  the  HF  rotational  distribution  and  hence  the  laser 
kinetics . 


1.2.  2.  2. 4  The  Effect  of  Changing  Both  H?  and  F? 
Concentrations  Simultaneously 

The  laser  pulse  energy  for  gas  mixtures  with  both  H2  and  F2 
varied  at  a  fixod  flash  lamp  energy  was  measured.  The  laser  energy  is 
coupled  out  with  a  sapphire  plate  with  a  copper  mirror  at  the  other  end 
The  results  are  shown  in  Fig.  15.  It  is  found  that,  in  general,  the  output 
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Fig.  15  The  Effect  of  Changing  Both  H2  and  F?  Concentration  to  the 
Laser  Output  Energy  Observed  in  the  H^/F^/lie  Mixture 
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energy  is  limited  by  the  smaller  constituent  of  the  mixture.  The  ratio  of 
[  laser  energy]  /[  F2  pressure]  in  an  active  volume  of  0.  1  liter  is  about 
80  joules/atm  F2,  which  is  about  twice  of  what  we  had  obtained  in  the  BaF2* 
copper  optical  coupling  experiment.  Since  both  BaF2  and  sapphire  windows 
ha  ve  about  the  same  transmission  and  reflectivity  at  2.7ji  region,  it  is 
certainly  surprising  to  note  this  difference  in  measured  laser  energy. 

More  discussion  of  the  effect  of  output  coupling  on  laser  energy  will  follow 
in  the  next  section. 


In  order  to  determine  if  the  difference  noted  above  also  exists 
for  other  window  materials  which  have  about  the  same  transmission  and 
reflectivity  at  2.7p  laser  energy  measurements  were  made  for  a  gas  mix¬ 
ture  of  H2/F2  He  =  0.  04/0.  04/0.92  with  a  copper  mirror  at  one  end  and  a 
variety  of  partial  transmitting  windows  at  the  other.  Table  I  lists  the 
averaged  laser  output  energy  observed  for  these  optical  cavities  using 
either  KBr,  NaCf,  LiF,  CaF2,  BaF2,  quartz  or  sapphire  as  the  partial 
transmitting  window.  All  of  these  cavities  have  about  90%  output  coupling. 
A  difference  in  output  energy  between  Ge,  IRTRAN  IV,  and  Si  flats  as 
partial  transmitting  windows  (which  give  about  50%  output  coupling)  was 
also  observed.  It  was  found  that  Si  is  the  best  of  the  latter  materials  and 
among  the  90%  output  coupling  windows,  sapphire  is  the  best.  However, 
in  comparing  sapphire  and  silicon  windows,  for  the  same  gas  mixture  and 
every  other  condition  the  same,  sapphire  couples  out  twice  as  much 
energy  as  that  of  silicon. 

The  performance  of  an  unstable  resonator  was  also  tested.  An 
unstable  cavity  should  have  the  advantage  that  the  mode  can  completely  fill 
the  active  medium,  thereby  minimizing  the  tendency  of  parasitic  oscilla¬ 
tion  and  also  allowing  efficient  energy  extraction  from  the  medium.  The 
unstable  cavity  consisted  of  a  2.  5  cm  diameter  concave  mirror,  and  a  1.  3 
cm  diameter  convex  mirror  mounted  on  a  barium  fluoride  output  window. 
This  cavity  had  14%  coupling  and  gave  a  pulse  energy  of  0.9  joules. 

Figure  16  summarizes  these  data  on  the  relationship  between 
laser  energy  and  output  coupling  observed  in  the  0.  04/0.  04,  0.  92  mixture. 
A  full  discussion  of  this  phenomena  will  be  given  in  later  sections. 


In  replacing  H2  by  D2,  the  laser  output  energy  was  slightly  de¬ 
creased  but  the  pulse  duration  was  found  to  be  either  unchanged  or  longer. 
Figure  17  shows  a  typical  DF  laser  pulse  observed  in  a  mixture  of  D7/F2/ 
He  =  0.  04/0.  04/0.  92  and  using  a  sapphire  output  coupling  window.  The  laser 
energy  monitored  by  the  AVCO  energy  meter  was  found  to  be  1.6  joule. 

Since  the  rates  of  vibrational  relaxation  of  DF  are  considerably  slower 
than  that  of  HF,  the  pulse  duration  is  expected  to  be  longer.  However,  the 
energy  of  1 . 6  joule  appears  to  be  too  low  in  comparison  with  the  value  of 
3.4  joule  obtained  in  the  H2/F2/He  mixture  using  an  identical  experimental 
condition.  If  equal  number  of  photons  were  emitted  from  both  gases,  a  DF 
laser  energy  of  3. 4/1.4  =  2.4  joules  should  be  obtained. 
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TABLE  I 


LASER  OUTPUT  ENERGY  FOR  VARIOUS  COUPLING  SCHEMES 


Material 

%  Transmission 
at  2.  7/1 

Measured  Laser  Output  Energy 
For  Cavity  Length,  L,  cm 

L  =  50 

L  =  25 

Sapphire 

87 

3.  6_+.  4 

- 

KBr 

92 

1 .  8  +.  3 

- 

NaCf 

92 

2.  4+.  3 

- 

BaF2 

93 
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Fig.  16  Laser  Output  Energy  as  a  Function  of  Ouiout  Coupling  11-14 
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Fig.  17  Time  Dependence  of  DF  Laser  Output  Pulse  Observed  in 
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The  transient  behavior  of  laser  photons  has  been  systematically 
studied  in  two  ways; 

i)  Comparing  the  total  pulse  shape  for  various  gas 
mixtures, 

ii)  Comparing  the  total  pulse  shape  for  various  out¬ 
put  couplings. 

The  results  are  summarized  as  follows: 

1.2.  2.  3.1  The  Effect  of  Gas  Composition  on  Pulse  Duration 

A  typical  oscillogram  of  the  HF  laser  pulse  (overall  transitions) 
observed  in  mixture  of  H2/F 2/He  =  0.04/0. 04 /0.  92  and  sapphire  output 
coupling  is  shown  in  Fig.  18.  The  upper  trace  is  the  flash  lamp  output, 
showing  a  pulse  length  less  than  2jj,secs.  The  laser  output,  shown  in  the 
lower  trace,  lasts  for  4gsecs.  The  fact  that  the  laser  pulse  is  longer  than 
the  flash  lamp  indicates  that  the  chain  reactions  are  operative  in  this  laser. 
From  the  overall  time  of  the  individual  chain  reactions  it  can  be  estimated 
that  approximately  10  chain  steps  are  possible  during  laser  action.  Figure 
17  shows  the  shape  of  a  typical  DF  laser  pulse  observed  from  a  D2/F2, 

He  =  0.  04/0.  04/0. 92  mixture  using  a  sapphire  output  coupling  window.  The 
observed  laser  duration  is  somewhat  longer  than  for  HF  -  about  6p,sec. 

Figure  19  shows  representative  HF  laser  pulses  observed  in 
various  gas  compositions.  It  is  clear  that  the  duration  of  laser  pulse  de¬ 
creases  with  increasing  H2/F2/He  =  0.  12/0.  08/0.  80  the  laser  duration  is 
only  2/j.sec  long,  but  for  mixtures  with  0.02/0.02/0.96  composition,  the 
duration  is  longer  than  5gsecs. 

1.2.  2.  3.  2  The  Effect  of  Output  Coupling  on  Pulse  Duration 

Figure  20  shows  the  shapes  of  HF  laser  pulses  obtained  in  mix¬ 
ture  of  H2/F2/He  =  0.  04/0.04/0.92  using  either  sapphire  or  silicon  output 
coupling  windows.  It  is  found  that  both  pulses  h:  ve  about  the  same  shape 
and  duration  regardless  of  their  differences  in  optical  transmission,  reflec¬ 
tion,  end  therefore  output  coupling. 

1.2.  2. 4  Spectra  Analysis  of  Pulsed  HF  Laser 

The  spectral  transitions  of  the  HF  laser  pulses  were  investigated 
in  two  different  optical  regions: 

i)  The  short  wavelength  infrared  region  between 
2.  5-3.  5  jj.,  and 


ii)  The  far  infrared  region  between  15-30  p.. 
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Fig.  20  The  Shapes  of  HF  Laser  Pulses  Observed  in  H?/F?/He  Using 
Either  Sapphire  or  Silicon  Output  Coupling  Windows 


The  results  obtained  by  the  experimental  technique  discussed  in  Section 

1 . 2,  1 . 2.  3  are  the  following: 

1.2.  2.4.1  Short  Wavelength  Infrared  Laser  Transitions 

Table  II  lists  the  various  P  and  R  branch  transitions  observed 
in  the  standard  H2/F2/He  mixture  using  a  sapphire  output  coupling  window. 
It  is  seen  that  vibrational  transitions  are  observed  for  V  =  6— 5,  5—4,  4  —  3, 
3—2,  2—1  and  1—0.  The  fact  that  6—5,  5—4  and  4  —  3  are  present  indicates 
that  the  hot  reaction  is  operative.  Furthermore,  the  observation  of  R 
branch  transitions  indicates  an  absolute  population  inversion  at  least  for 
certain  vibrational/ rotational  states. 

Figures  21,  22,  23,  24,  25,  26  and  27  show  the  temporal  be¬ 
havior  of  the  IP,  2P,  3P,  4P,  5P,  6P  and  R  branch  transitions  observed  in 
the  standard  mixture  of  H2/F2/He  =  0.04/0  04/0.92  using  a  sapphire  out¬ 
put  coupling  window.  There  are  two  features  that  are  noteworthy; 

i)  Among  all  the  P  branch  transitions,  there  are 

definitely  two  noticeable  trends  that  reach  lasing 
threshold.  For  transitions  involving  low  J  states 
(below  J=8),  the  lowest  J  transition  alw'ays 
reaches  threshold  first  and  higher  J  transitions 
occur  later.  For  transitions  involving  high  J 
states  (above  J=8),  there  are  definitely  two  pulses 
in  each  transition.  A  pulse  occurs  at  early  time, 
and  the  second  pulse  arrives  later  (e.  g.  ,  see 
Fig.  21,  lPg  transition).  The  first  pulses  reach 
threshold  almost  as  early  as  the  lowest  transi¬ 
tion  for  the  low  J  transition  (below  J=8).  The 
second  pulse  appears  to  follow  the  temporal  trend 
of  the  low  J  transitions. 


ii)  The  R  branch  transitions  reach  threshold  quite 
early  and  only  lase  on  the  low  J  states  where  the 
P  branch  lasing  action  does  not  occur.  This  is 
similar  to  what  Prof.  John  Polanyi  predicted  six 
years  ago.  (Hk) 

1.2.  2.4. 2  Far  Infrared  Laser  Transitions 

HF  laser  radiation  in  the  far-IR.  region  has  been  reported  from 
both  the  SFfj  +  H2  and  CF4  +  H2  reactions/^)  Since  (a)  the  H2/F2  reaction 
favors  the  population  of  high  rotational  states  for  lower  vibrational  levels 
and  (b)  the  rate  rotational  relaxation  of  the  high  J  states  of  HF  could  be 
very  long,  i.e.,  comparable  to  the  vibrational  relaxation  time,  it  is  quite 
possible  that  population  iiversion  may  also  exist  between  rotational  energy 
levels.  In  order  to  determine  if  there  was  any  far-IR  radiation  in  our 
H2/F2  laser  system,  a  black  polyethylene  cloth  was  used  to  cover  the 
energy  meter,  and  the  transmitted  energy  was  measured.  It  was  found 
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TABLE  II 


OBSERVED  LASER  TRANSITIONS 


BAND 

P  BRANCH 

R  BRANCH 

l  -  0 

P  (3)  (4)  (5)  (|)  (7)  (8)  (9) 

NONE 

2  -  1 

P  (3)  (4)  (5)  (6)  (7)  (8)  (±)  (10)  (li) 
(12)  Tl3)  (14)  (15)  (16) 

R(0) 

3-2 

P  (4)  (5)  (6)  (7)  (8) 

NONE 

4-3 

P  (4)  (5)  (6)  (7) 

R  (1),  R  (4) 

5-4 

P  (3)  (4)  (5)  (6) 

R  ( 0 ) 

6-5 

P  (1)  (2)  (3)  (4)  (5) 

NONE 

7-6 

NONE 

NONE 

8-7 

NONE 

NONE 

*  Underlines  indicate  the  strongest  transition  observed  in  the  band 
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that  for  a  total  1.  5  joule  HF  laser  pulse,  about  0.  3  joules  of  energy  were 
transmitted  through  the  black  polyethylene.  Since  black  polyethylene  only 
starts  to  transmit  above  15(i,  these  0.  3  joult  s  of  energy  are  the  contribu¬ 
tion  of  laser  photons  with  a  wavelength  greater  than  15|iand  mos*  likely 
above  the  20  /j.  region.  Unfortunately,  without  a  spectroscopic  study  in  the 
far-IR  region,  it  is  not  possible  to  conclude  which  particular  rotational 
levels  are  involved,  and  to  determine  if  the  rotational  laser  action  has 
proceeded  before,  during  or  after  the  vibration- rotational  laser  transitions 
at  shorter  wavelength.  This  spectroscopic  study  was  not  performed. 

1.2.  2.  5  F?  Dissociation  Produced  by  Flash  Photolysis 

The  result  obtained  by  the  titration  technique  discussed  in 
Section  1 . 2.  1 . 2.  1  is  shown  in  Fig.  28.  This  plot  of  F/F2  ratio  vs  F2 
pressure  in  the  HCf/F2/He  mixture  shows  a  constant  value  over  the  F2 
pressure  jtudied.  HCf  enriched  mixtures  were  used  to  assure  that  all  the 
F  atoms  photolized  by  the  flash  reacted  with  HCf  to  form  molecular  HF. 

The  data  show  a  dissociation  of  F2  of  about  1%  independent  of  F2  partial 
pressure.  This  lack  of  F2  concentration  dependence  is  interpreted  as  a 
result  of  the  gas  being  optically  thin  under  these  experimental  conditions. 

In  order  to  determine  if  side  reactions  like 

Cf  +  F2  -  CfF  +  F  (11) 

also  occur  and  produce  additional  HF  molecules,  measurements  with  both 
HCf  and  F2  enriched  mixtures  were  tested.  It  is  found  that  the  amount  of 
HF  formed  does  not  increase  linearly  or  quadratically  with  the  concentra¬ 
tion  of  HCf  present.  Figure  29  shows  the  fraction  of  F-atoms  produced 
as  a  function  of  HCf  pressure.  Above  2  torr  of  HCf,  the  data  flatten  out 
and  yield  a  constant  result,  namely,  1%  F2  dissociation  or  F/F2  =  (8  +  2)  x 
10‘3.  Based  on  our  previous  photylysis  results  using  the  linear  flash 
lamps  and  less  efficient  optical  coupling,  the  coaxial  lamps  appear  to  yield 
about  a  factor  of  two  more  F-atoms. 

1  .  2 .  2 .  6  The  Rate  of  H 2  +  F 2  —  2HF  Overall  Reaction 

The  results  obtained  from  the  experiment  described  in  Section 
1 . 2.  1 . 2.  2  are  shown  in  Fig.  30  where  the  ratio  of 

(fF2l  -|F])/|F] 

\  t=o  t/  t=0 

is  plotted  against  the  time  t  after  flash  lamp  initiation.  By  adopting  the 
F /F 2  ratio  previously  determined  and  the  rate  constants  of  both  "hot"  and 
"cola"  reactions  reported  in  the  lite  ratur  e ,  ( *  3 )  the  rate  of  the  overall 
H2  +  ^2  -»  2HF  reaction  is  computed  and  plotted  along  with  the  experimental 
values.  It  is  clearly  evident  that  the  consumption  of  F2  molecules  in  the 
cavity  proceeds  with  a  speed  comparable  to  that  predicted  from  the 
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theoretical  computations  not  using  a  chain  branching  mechanism.  The 
calculated  value  using  the  chain  branching  Reaction  (2)  with  k2  =  10"  12 
cm3/sec  is  also  plotted  in  Fig.  30.  It  is  seen  that  the  overall  H2  +  F2 
reaction  would  proceed  with  a  much  faster  speed  than  that  observed  ex¬ 
perimentally  if  a  fast  chain  branching  reaction  is  included  in  the  mechanism. 


1.2.  2. 7  Determination  of  Beam  Distortion  Inside  the 
Laser  Cavity 

Since  (a)  the  H2/F2  reaction  liberates  significant  amounts  of 
energy  (-35%)  into  the  translational  degree  of  freedom,  and,  (b)  the  rate 
of  HF  vibrational  relaxation  is  extremely  fast,  the  temperature  of  the 
cavity  gas  must  increase  rapidly  with  time  after  flash  lamp  initiation.  It 
is  possible  that  the  photons  traveling  within  the  cavity  between  the  window 
and  mirror  could  be  seriously  scattered  because  of  possible  density 
gradients  in  the  gas,  and  thus  the  beam  may  not  be  regenerated  inside  the 
optical  cavity.  To  test  this  hypothesis  a  simple  light  scattering  experiment 
was  performed.  Figure  31  shows  the  experimental  arrangement.  The 
beam  from  an  He-Ne  laser  operating  at  6328  R  was  aligned  within  the 
cavity  and  reflected  between  the  two  windows.  After  two  passes  within  the 
cavity,  the  beam  was  reflected  onto  a  narrow  slit,  and  the  intensity  was 
monitored  by  a  photomultiplier  behind  the  slit.  The  purpose  of  the  exper¬ 
iment  was  to  determine  if  serious  density  gradients  could  exist  within  the 
cavity  gas  during  the  time  when  laser  action  had  been  seen.  The  mono¬ 
chrometer  was  used  before  the  PM  to  reduce  flash  lamp  scattered  light. 

Figure  32  shows  the  temporal  intensity  variation  of  the  re¬ 
flected  beam  after  the  flash  initiation  of  a  standard  H2/F2  mixture.  A 
periodic  oscillation  of  the  beam  was  observed.  By  changing  the  gas  com¬ 
position  from  H2/F2/Hc  to  H2/F2/Ar,  the  oscillating  period  was  found  to 
be  decreased  roughly  from  1 9  /J.fcc /cycle  to  36  gtsec /cycle  .  Since  the 
sound  speed  of  He  at  1  atm  pressjre  is  about  1  mm/jnsec  and  the  flash  lamp 
diameter  is  about  17  mm,  the  value  of  19psec/cycle  seems  comparable 
with  the  period  of  sound  wave  oscillation  between  walls  of  flash  lamp.  The 
observation  that  the  oscillations  are  slower  for  the  Ar  mixture  is  also 
consistent  with  the  reduced  sound  speed  in  that  gas.  Further  discussion  of 
the  possible  significance  of  this  effect  on  laser  performance  will  be  pre¬ 
sented  later. 


1 . 2.  2.  8  Examination  of  the  Laser  Burn  Pattern 

In  order  to  determine  (i)  some  information  concerning  the  mode 
structure  of  this  H2/F2  laser,  and  (ii)  to  determine  if  parasitic  oscillations 
were  involved  in  our  laser  system,  the  burn  pattern  of  the  laser  output  was 
taken  and  examined.  Figure  33  shows  the  laser  radiation  intensity  distribu¬ 
tion,  obtained  at  a  distance  of  40  cm  from  the  output  sapphire  window.  As 
one  can  clearly  see,  the  burn  mark  is  extremely  bright  at  the  center  of  the 
pattern.  The  diameter  of  this  circular  pattern  is  about  3  mm  smaller  than 
the  inside  diameter  of  our  laser  tube.  Since  our  optical  cavity  is  composed 
of  an  8  m  radius  concave  mirror  and  flat  reflector,  the  Fresnel  number,  N 
for  this  semi-confocal  resonator  can  be  calculated  as 
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He-Ne  LASER 


Fig.  31  The  Optical  Set  Up  Used  for  Determination  of  Beam  Distortion 
Int  de  the  Laser  Cavity 


/-SCATTERED  LIGHT  FROM  FlASh 


t 


GASES  FLOWING  BUT  NO  FLASH  FLASH  L A M P  FIRED  INTO  I  ATf/ He 

5^.SF.C/CM 


DISTORTED  SIGNAL  FOR  1/1/23  MIXTURE  INITIATED  BY 


0136’ 


20/xSEiVCV 

Dl  STORTED  SIGNAL  DlS  ’LACED  WITH  A 
LONGER  TIME  SCALE 


Fig.  32  The  Temporal  Intensity  Variation  of  the  Reflected  Beam  after 
the  Flash  Initiation  of  H^/F ^  Mixture 
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where  r  is  the  radius  of  the  mirror  (2.  5  cm)  and  L  is  the  distance  between 
reflectors  (75  cm).  The  radius  of  the  spot  size  can  be  estimated  as 
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where  b  is  the  radius  of  curvature  of  the  concave  mirror.  Sj  and  S2  are 
the  spot  sizes  on  the  two  mirror  surfaces.  The  diameter  of  a  multimode 
laser  spot  should  be  about  17-1.5  =  15.5  mm  which  is  precisely  what 
we  have  found  (see  Fig.  33).  If  parasitic  modes  were  dominating  our  laser 
system,  the  size  and  the  diameter  of  the  burn  mark  should  be  considerably 
larger  due  to  beam  divergence,  which  is  clearly  not  what  we  have  seen. 


1.2.  2.  9  Estimate  of  the  Rate  of  Rotational  Relaxation  of  HF 

In  our  last  contract  report(  5  )  we  had  clearly  pointed  out  that, 
during  the  duration  of  the  HF  vibration- rotational  laser  action,  the  rota¬ 
tional  degree  of  freedom  of  the  HF  molecules  may  not  have  reached  equi¬ 
librium  with  the  translational  mode.  The  fact  that  far-IR  laser  radiation 
was  observed  in  our  H2/F2  system  as  described  in  Section  1. 2.  2.4  strong¬ 
ly  suggests  that  a  measurement  of  the  rate  of  rotational  relaxation  of  HF 
is  important  to  understand  and  model  this  laser  system. 

Figure  34  shows  the  experimental  arrangement  of  a  double - 
resonance  absorption  experiment  that  we  assembled  to  attempt  to  measure 
rotational  relaxation  of  HF.  Two  electric  HF  pulsed  lasers  utilizing  the 
SF(j/H2  reaction  were  used.  Laser  #  1  puts  out  5  mj/cm^  with  30%  of  this 
energy  in  the  v  =  1  -*  0  transitions  (10%  IP4,  10%  IP5,  6%  1P(,  and  4%  IP7). 
The  IP7  photons  were  chosen  to  excite  the  HF  molecules  in  the  cell.  Laser 
#2,  fired  with  a  variable  delayed  time  after  Laser  #1,  was  operated  on  the 
2P&  transition,  and  was  also  positioned  to  pass  through  the  same  volume  of 
gas  that  was  excited  by  Laser  #1.  The  amount  of  2Pfc  absorption  was 
monitored. 


Unfortunately,  due  to  the  limitations  of  time  and  money,  this 
measurement  was  not  completed.  Some  delay  between  the  IP7  pump  and 
2P(,  probe  pulses  were  observed,  but  with  the  very  low  pressures  of  HF 
used  (<0.  1  torr)  it  was  difficult  to  determine  the  absolute  concentration  of 
HF  and  hence  deduce  a  rate.  However,  with  our  crude  data,  we  can  con¬ 
clude  that  the  relaxation  of  the  HF  (v  =  1 )  state  by  processes  like 


X 
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HF  LASER  #l 


Fig.  34  The  Experimental  Arrangement  of  a  Double  Resonance  Absorption 
Experiment  Used  for  Measuring  the  Rotational  Relaxation  of  HF 
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HF  (v  =  l,  J  =7)  +HF  (v  =0,  J  =  2) 

-  HF  (v  =  0,  J  =  7)  +  HF  (v  =  1,  J  =  2)  (15) 

+  A  E  a  0 


and 


HF  (v  =  1,  J  =  1)  +  HF  (v  =  0,  J  =  2) 

—  HF  (v  =  1,  J  =  6  or  5)  +  HF  (v  =  0,  J  =  2,  3,  or  4)  (16) 

+  AE  =  small  but  larger  than  zero 

can  be  extremely  fast  and  takes  only  a  few  collisions  or  could  even  be  as 
fast  as  gas  kinetic.  It  should  be  informative  to  measure  the  relaxation 
time  of  the  HF  (v  =  0,  J  >  10)  states  where  the  exact  v  —  v  energy  transfer 
will  not  be  present  and  only  R  -*  R  or  R  -►  T  could  partic  ipate  in  the  relaxa¬ 
tion. 
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1.  3  MODEL  CALCULATION  FOR  THE  PULSED  HF  CHEMICAL  LASER 


1.3.1.  Basic  Philosophy 

Theoretical  modeling  of  the  pulsed  H2/F2  laser  was  performed 
for  the  following  reasons: 

(i)  To  see  if  a  direct  comparison  between  the  experimental 
results  (such  as  pulse  shape,  output  spectra  and  output 
energy),  and  predications  from  a  simplified  theoretical 
model  could  lead  to  sufficient  understanding  of  laser 
performance  so  that  the  potentiality  of  a  scaled-up  laser 
system  can  be  assessed  without  involving  too  much 
additional  expe  riT.iental  effort. 

(ii)  Since  this  is  one  of  the  few  laser  systems  in  which  the 
rates  of  a)  stimulated  emission,  b)  chemical  reactions, 
c)  rotational  relaxation  and  d)  vibrational  relaxation  can 
be  comparable,  a  parametric  or  systematic  study  on 
their  effects  relative  to  the  theoretical  laser  performance 
should  be  of  crucial  importance  and  value  to  the  under¬ 
standing  of  fundamental  laser  physics. 

1.3.2.  Reaction  Kinetics 


The  major  processes  involved  in  H2/F2  laser  reaction  can  be 
summarized  as;  (a)  the  photodissociation  reaction  used  for  initiation,  (b) 
ihe  chain  reactions  producing  vibrationally  excited  HF  molecules,  (c)  the 
various  vibrational  relaxation,  rotational  relaxation  and  energy  transfer 
processes  of  the  excited  HF  molecules  and  (d)  the  stimulated  emission 
producing  the  laser  actions.  Table  III  lists  these  major  reactions.  Ihe 
rate  constants  o'  most  of  the  important  reactions  are  known,  but  some 
critical  rates  are  not  known.  5ince  a  modeling  calculation  excluding  these 
unknown  but  important  reactions  is  already  quite  complicated,  the  compu¬ 
tation  was  first  performed  without  consideration  of  reactions  (VIII)  (IX)  (X). 
Later  on,  these  effects  are  s'parately  tested. 

1.3,3  Basic  Assumptions 

In  order  to  simplify  the  calculations  a  few  assumptions  were 

made: 


(i)  For  a  given  vibrational  energy  level,  two  extreme  rota¬ 
tional  distributions  are  corr pared; 

a)  A  completely  equilibrated  Boltzmann  distribution  at 
the  translational  temperature  which  assumes  that 
rotational  relaxation  is  extremely  fast. 


b)  A  completely  un^elaxed  rotational  distribution,  assum¬ 
ing  that  rotations  1  relaxation  is  very  slow  compa  red 


TABLE  III 

REACTIONS  INCLUDED  IN  KINETICS  MODEL 


(I) 

Photodis  sociation 

F2  +  hw  -  2F 

(II) 

Chemical  Production 

F  +  H2  -  HF*  (v)  +  H 

H  +  F2  —  HF*  ✓)  +  F 

(III) 

Vibrational  Relaxation 

HF  (v)  +  M  —  HF  (v  -  1)  +  M 

M  =  HF,  H2,  F2,  N2,  He,  F,  H. 

H2(v)  +  M  -*  H2  (v  -  1)  4-  M 

(IV) 

Vibrational  Energy 

HF  (v)  +  HF  (v')  — *  HF  (v  +  1)  +  HF  (v’ 

-  1) 

T  ransfer 

HF  (v)  +  H2  (v1)  — ♦  HF  (v  +  1)  +  H2  (v* 
H2  (1)  +  H2  (1)  *-*  H2  (2)  +  H2  (0) 

-  1) 

(V) 

Recombination 

H+F+M-HF+M 

f+f+m-f2  +  m 

H  +  H  +  M  -*  H2  +  M 

(VI) 

Branching 

HF  (v)  +  F2  -  2F  +  HF  (0) 

(VII) 

Stimulated  Emission 

HF  (v)  +  hv  -*  HF  (v  -  1)  +  2hv 

RATES  UNKNOWN  BUT  IMPORTANT 

(VIII)  Effect  of  Vibrational  HF  (v)  +  HF  (v')  —  HF  (v  -  1)  +  HF  (v') 
Excitation  upon  the  Rates 
of  Vibrational 
Relaxation 

(IX)  Rotational  Relaxation  HF  (v,  J)  +  M  —  HF  (v  J')  +  M 

(X)  Vibration-to-Rotational  HF  (v,  J)  +  M  —  HF  ( v -  1 ,  J ' )  +M 
Energy  T  ransfer 
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with  the  rate  of  stimulated  emission  and  during  the 
lasing  period  the  molecules  still  maintain  their  original 
rotational  distribution. 

(ii)  For  each  vibrational  level  only  one  rotational  line  (that  with 
a  maximum  gain  i.e.  a  P-branch  transition)  is  capable  of 
laser  action  at  any  one  time.  However,  /arious  vibrational 
levels  may  lase  simultaneously. 

(iii)  A  homogeneous  photodissociation  of  F£  molecules. 

(iv)  Both  stimulated  emission  and  vibrational  relaxation  are 
limited  to  single  quantum  transitions. 

(v)  It  is  assumed  that  the  gain-equal- loss  condition  holds  during 
laser  oscillation. 

(vi)  It  is  further  assumed  that  the  line  width  of  the  laser  transi¬ 
tions  is  mainly  due  to  collisions  with  He,  HF,  H2  molecules 
and  hole  burning  within  the  Doppler  line  width  is  not  con¬ 
sidered. 

Among  all  these  assumptions,  Nos.  (iii)  and  (iv)  are  reasonably 
close  to  reality.  Number  (ii)  is  in  obvious  contradiction  with  the  experi¬ 
mental  results  which  show  that  often  more  than  two  vibration- rotation  lines 
oscillate  simultaneously.  However,  this  approximation  will  not  appreciably 
alter  the  total  laser  energy  or  laser  duration.  Assumption  No.  (i)  is  un¬ 
avoidable  because  of  the  lack  of  knowledge  of  rotational  relaxation  rates  in 
HF,  but  should  provide  limits  by  which  to  compare  the  experimental  data. 
Number  (v)  is  not  a  bad  assumption  for  small  gain  lasers.  However,  for 
laser  transitions  with  extremely  high  gain  or  supe rradiation ,  the  results 
estimated  under  this  condition  can  only  be  considered  as  a  lower  limit. 
Number  (vi)  may  be  correct,  unless  atoms  like  H  or  F  have  extremely  large 
collisional  broadening  cross  sections.  For  high  pressure,  viscous  gas  flow 
(which  is  the  condition  of  our  HF  laser),  hole  burning  should  not  be  important 
at  all. 

1.3.4  Description  of  the  Model 

The  theoretical  model  is  basically  that  used  by  Airey^^for 
the  pulsed  HCf  laser  but  extended  to  many  vibrational  energy  levels  for  the 
HF  system.  The  model  assumes  that  for  a  particular  vibrational  level  only 
that  P-branch  transitions  with  the  greatest  gain  can  lase  and  that  gain  can 
not  become  greater  than  loss. 

As  in  the  experiments,  the  reaction  is  initiated  by  ultraviolet 
light  which  photodissociate s  some  of  the  fluorine  molecules  present,  i.e. 
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+  hv  —  2F  (I)* 

Following  which  th',  reaction  proceeds  via  the  chain 

F  +  H?  —  HF  (v)  +  H  -  AH  =  32  kcal/mole 

(II) 

H  +  F2-HF  (v)+  F  -AH  =  98  kcal/mole 

The  first  reaction,  called  the  "cold"  reaction,  has  sufficient  energy  to  pop¬ 
ulate  vibrational  levels  of  HF  up  to  v  =  3,  and  the  second,  or  "hot"  reaction 
may  populate  up  to  v  =  8.  Roughly  half  of  the  reaction  energy  appears  as 
vibrational  energy.  This  energy  may  be  redistributed  through  (v-v)  colli¬ 
sions,  e.  g.  , 


HF  (v)  +  HF  (v')-HF  (v+1)  +  HF  (v'  -  1)  (IV) 

or  may  be  converted  to  translational  energy  (or  more  likely  rotational 
energy)  through  collisions,  e.  g.  , 

HF  (v)  +  M  —  HF  (v  -  1)  +  M  (HI) 

Obviously  the  vibrational  energy  may  also  be  removed  as  laser  energy. 


For  each  vibrational  level,  the  population  is  determined  by 
integrating  with  respect  to  time  an  expression  of  the  form 


dN 

v 

dt 


where 


V  ,  -  A 

2J  -  1) 
V _ 


N 

v 

2J  + 


*  1. 

iy 


■) 


jp 


v-l 


B 


v-l 

V 


production,  vibra¬ 
tional  relaxation 


stimulated  emission  from 
v  level 


v+1. J1-  1 
2 J'  -1 


N 


v,J' 


2J'  +  1 


J'p 


v-l 


,  V-  1 


v' - 

stimulated  emission  from 
level  v  +  1 


J 


(17) 


Nv  is  the  population  of  the  v^  vibrational  level 
k^  represents  rate  constants  of  various  reactions 

♦  These  Roman  numerals  relate  to  the  reactions  as  listed  in  Table  HI 
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X  and  Y.  are  molecules  involved  in  reaction  i 
1  i 

N  is  the  population  of  the  rotational  level  of  the  v1^  vibrational 

vJ  level 


Pv 


v”  ^  is  the  radiation  density  from  level  v  to  v-1  (ergs- sec/cm  ) 

B* is  that  part  of  the  Einstein  coefficient  which  is  J  independent 

V 

The  gain  of  a  particular  P-branch  transition  can  be  written  as: 


,  v- 1 


N 


G(v,J)  = 


'  v-1.  2 

Av  t(2ttu  ) 
v,J  v 


v,  J-l 
2  J-  1 


(18) 


where 

Av  ,  is  the  full  line  width  at  half  maximum  (sec 
v,  J 

wV”  1  is  the  wave  number  of  the  line  center  (cm  ) 
v 

AV"  1  is  the  Einstein  probability  of  spontaneous  emission  corres- 
v  ponding  to  the  integrated  band  intensity 

The  Einstein  constants  B  and  A  are  correlated  as 

Cl  ■  3)  (.9, 

During  laser  action,  the  gain  is  set  equal  to  the  cavity  loss,  i.e. 

R  (1  -  Lc)  exp  (2  G  C)  =  1  (20) 

where 


Lc  is  the  output  coupling  fraction 

R  is  the  "effective"  reflectivity  (v'hich  includes  all  useless  losses) 
t  is  the  active  length  of  the  cavity. 

The  calculational  procedure  is  simplified  as  follows: 

(a)  The  quantities  1^  j  and  defined  as  follows: 


N 


v,  J-  1 


2J  -  1 


N 


v-1,  J 


2  J  +  1 


(21) 
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and 


a  /->  v- 1 .  2 

Av  (2  7Tu>  ) 

'  V  ' _ 

21  Av" 1 

V 

are  calculated.  The  line  width,  A  v ,  is  taken  to  be  independent 
of  quantum  number. 

(b)  At  any  step  in  the  calculation,  the  maximum  value  of  the  quantity 
I  T  is  found  for  each  vibrational  level,  i.e.  , 

C  =  {k.j}  <23» 

J=  1 

It  is  assumed  that  lasing  action  from  a  particular  vibrational 
state,  v,  is  only  possible  if 


M 


=  Cn  £- 


R  (1  - 


T7 


Im>M  (24) 

V  V 

Since  the  quantity  loss  per  path  is  a  constant  value,  the 
condition  that  gain  may  not  become  greater  than  loss  leads 
to  the  derivation  that  wherever  the  above  condition  (Eq.  (24)) 
is  met; 


dG(v,  J) 


d  I 


m 


dt 


0,  or 


dt 


=  0 


dN 


The  expression  of  — ^ -  can  be  rewritten  as 


dN 


dN 


dt 

where 


dt 


N  =  L 


J=o 


,  jm  v  _v 

v+1  Pv+1  v+1 


.  I 


m  v- 1  „v-l 


v  ^  v 


B 


N 


v,  J 


n„  E  (J); 


V  —  v 
J=o 


dN 

_ v 

dt 


(25) 


(26) 


(27) 


is  the  contribution  to  the  rate  of  change  of  vibrational 
level  v  due  to  collisional  processes  and  spontaneous 
emission. 
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(c)  Combining  Eqs.  (25)  and  (26)  with  the  assumption  that  p^'  *  =  o 

gives  a  set  of  algebraic  equations  for  p^"  *  and  dNy/dt  if  condition 
(Eq.  (24))  is  not  met.  Solving  these  simultaneous  equations  give 
us  pU-1. 


(i)  For  the  case  of  the  equilibrium  distribution  of  rotational 

n  la),  fv(J)  of 

exp  (  -  J  (J+ 1] 


energy  levels  (assumption  la),  f  (J)  of  Eq.  (27)  can  be 
written  as: 


fv(J)  =  (2J 


...  I 


(-J(j+h/Qr  v)/ 


'Rv 


! 


where  Q^v  i-8  the  rotational  partition  function: 


(28) 


FT 


Rv 


hci3 


(29) 


with  Bv  the  rotational  constant  of  vibrational  state  v. 
This  gives  the  following  expression  for  Eq.  (25): 


9  I 


m 
v 
a  N 


dN 


m 


dt 


9  I 

v 

In 


rm 


v-1 


dN  .  9  1 

v- 1  ,  v  dT  «  /ic  v 

— dt  +~irr  sr  =  0  (25a) 


(ii)  For  ilie  case  of  frozen  rotation  (assumption  (ib)),  the 

rotational  distribution  function  can  be  determined  by  the 
detailed  rate  constants  of  the  chemical  production  reactions; 


F4HZ  lv(J)>  HF  (v,  J)  +  H 
H+f  — 2W1L>  HF  (v,  J)  +  F 


(II) 


The  production  of  a  particular  vibration- rotation  level  is 


d  [HF  (v.  ,  [F]e,  j  (J|  t[H][F2]  k.  (J) 


(30) 


and  the  distribution  function  f  (J)  of  Eq.  (27)  is  then 
assumed  to  be 


f  a) ,  jChf  <y.j?3 


v'J> - ar 


L 

J 


d  [hF(v,  J)] 


dt 


[F][H2]klv  (J)4[h][F2]  k2v(J) 

£  {[f][h2]  klv  (J)  +  [h][f2]  k2v(j)| 


(31) 
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Since  the  characteristic  time  for  removal  of  H  atoms  by 
the  hot  reaction  (the  slower  of  the  two  reactions)  is  l/4jj. 
sec  for  the  conditions  of  interest,  and  this  time  is  less 
than  the  time  for  production  of  HF.  (The  production  of  HF 
is  limited  by  the  production  of  F-atoms  by  the  flash  lamp.  ) 
A  steady  state  is  set  up  between  the  two  reactions  such  that 

^1=  £  [f][h2]  klv(J)  -  £  [h][f2]  k2v(J)*0(32) 

V,  J  v,J 


which  leads  to: 

k 


f  (J)  = 

v'  ' 


w«j»a1t  +k2v<J|/t 


2T 


‘1v/l 


+k 


"1 T 


(33) 


2T 


where 


lv 


=  £  klv(J), 


L  kiv 


(34) 


TT  Lj  *Tv 

V 


This  distribution  functi  in  is  thus  independent  of  time  and 
temperature  and  the  last  term  in  Eq.  (25a)  drops  out 
because 


3  iJP 
St 


=  0 


(35) 


(d)  Throughout  the  integration  the  gas  temperature  T  is  calculated 
by  integrating 


In  which  Cv  is  the  specific  heat  at  constant  volume;  E^  and  E2 
are  those  portions  of  the  reaction  energies  for  the  cold  and  hot 
reactions,  respectively,  which  go  directly  into  translation;  kj'p 
and  ^2T  are  ^e  rate  constants  for  the  cold  and  hot  reactions, 
respectively,  E^  is  the  portion  of  the  absorbed  photon  energy 
which  appears  as  translational  energy  (photon  energy  minus 
dissociation  energy);  and  the  final  term  is  the  contribution  due 
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to  deactivation  of  vibrational  energy.  This  translational  temper¬ 
ature  is  used  in  calculating  those  rate  constants  which  depend  on 
temperature,  and  in  calculating  the  gain. 

(e)  Finally,  for  converting  the  resulting  values  of  into  the  power 

output  in  watts/cm^,  a  multipliative  factor  Cc  is  used,  where 


C 

c 


10"7L  Av  3  x  108 
c 

U 


The  total  energy  output  in  joules/cm  is  then 


(37) 


(38) 


1.3.5  Major  Rate  Constants 

Clearly  the  validity  and  usefulness  of  the  computation  depends 
heavily  on  the  accuracy  of  the  rates  used.  The  values  of  the  major  rate 
constants  listed  in  Table  III  are  given  below  in  units  of  cm^/particle- sec. 

(a)  Production  (II) 


F+H 

k 


2 

10 


k 


:k 


— >  HF(v)  +  H 

1 1  ;ki  2:ki3;ki4  =  0:0.29:1.0:0.47:0:0: 
4 

L  k  =  2.  6  x  10'10  exp  (-800/T) 


*\ 

H+F2  — >  HF  (v)  +  F 

k20:k21:k22: . k29  =  0:0.  13:0.  16:0.  29:0.  35:0.76:1.0:0.  38:0.  25:0 

10 

£  k2v  =  2.  2  x  1  O'  1  °exp  (-  1  200/T) 
k3v 

H+HF(v)— >  F  +  H2 

k3v  =  °’  forvs3 

k_  =  5  x  10  for  v>4 
3v 


-66- 


(b)  Non-reactive  processes 


(i)  Vibration  to  translational  energy  transfer  rates  (III) 

KF(v)  +  F  -  HF(v-l)  1-  F  ,  kv  =  2.  5x10"  12  v11 

“12  n 

HF(v)  +  H  —  HF(v-l)  +  H  ,  ky  =  2.  5x10  v 

HF(v)  +  F2  —  HF(v-l)  +  F2  ,  ky  =  2.  5xlO"15  vn 

HF(v)  +  H2  —  HF(v-l)  +  H2  ,  kv  =  2.  5xl0‘15  vn 

HF(v)  +  HF  —  HF(v-l)  +  HF  ,  ky  =  8x10"  13  exr>(785(T“  -  .  105)2)  vn 

where  n  =  1  or  2 

(ii)  Vibration-vibrational  energy  transfer  rates  (IV) 

k+ 

HF(v)  +  HF(v'  )  =  HF(v- 1)  +  HF(v'  +  1)  for  v  >  v'  +  l 

k" 

Two  sets  of  rate  constants  have  been  tried:  1)  quantum 
number  dependent  rates  and  2)  energy  discrepancy 
dependent  rate  constants: 

1)  k'(v,v')  =  2.5  x  10  12v(v'  +  l) 

or 

2.5xl0"12,  AE<  300  cm_1 

2)  k“(v,  v')=  2.5xl0-13,  300<  AE  <600  cm-1 

2.  5xl0'14,  600  <  AE 

and 

k+(v,  v' )  =  k"exp(-AE/kT) 
k+ 

HF(v)+H,(v')-  HF(v+l)  +  H?(v»  - 1) 

2  k" 
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where  k“(v,  v'  )  =  k+(v,  v' )  exp  (-  AE/kT) 

(c)  Other  Rate  Constants 

The  reactions  described  above  are  those  considered  to  be  the 
most  important  for  the  H2/F2  chain  reaction.  However,  the 
other  reactions  photodissociation  of  F2(I),  recombination 
reactions  of  atoms  (V),  and  chain  branching  reactions  of  HF 
by  F?(VI1  are  definitely  involved  and  also  important.  The 
following  rate  constants  for  these  other  reactions  have  been 
adopted  in  the  calculation: 

F2  +  hv  -  2F 

e  =  absorption  coefficient  =  6  litre/mole-cm 
H  +  F  +  M  -  HF  (0)  +  M  ,  k=  7.2x10"  /T 
F  +  F+  M-F2  +  M,  k=3.  0x10"32/T 

H+H+M-H2+M,  k =  6.6x10'30/T 

HF(v)  +  F2  -  2F  +  HF  (0)  . 

k(v  <  4)  =  10"1  'T  1//2  exp(Ev-E)/RT 
k(v  >  4)  =  10"1?T  1//2 

where  E  =  37  kcal/mole  and  E  =  vibrational  energy 
level  of  v  state 

1,3.6  Systematic  Studies  on  Various  Parameters 
1.  3.  6.  1  Line  Width 

As  one  can  clearly  see  from  both  the  gain  Eq.  (18)  and 
output  power  expression  (38),  the  line  width  Avvj  has  direct  influence  on 
the  calculated  gain  and  laser  output  energy.  Consequently,  in  the 
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computation,  the  value  of  Avyj  has  to  be  carefully  chosen.  Meredith 
et  al^ 5)h? «  '  recently  published  their  calculated  line  widths  for  HF-HF 
collisions.  At  400°K,  a  peak  line  width  of  0.42  cm"  /atm  for  J=2  was 
found.  The  value  drops  off  sharply  as  J  increases.  For  example,  a 
value  of  0.08  was  found  for  J  =  8  state.  For  our  standard  gas  mixture, 
with  a  maximum  final  HF  pressure  of  61  torr,  the  maximum  broadening 
should  be  around  0.04  cm"*,  For  reactions  with  33%  completion  (the 
approximate  end  of  laser  action),  the  maximum  HF  pressure  will  be  less 
than  20  torr,  and  the  maximum  line  width  broadened  by  HF=HF  collisions 
will  be  less  than  0.  013  cm”1. 


Meredith  et  al  have  also  calculated  the  line  width  values  for 
collisions  of  HF-H2  ,  HF-N2  and  HF-He.  It  was  found  that  y  (line  width) 
values  decreases  only  slowly  with  the  rotational  quantum  number  J.  For 
HF-H2  collisions  aline  broadening  value  of  y  =  0.092  cm'Vatmwas  found 
for  J=  1,  y=  0.  062  for  J=4  and  y=  0.  05  for  J=8.  The  J  dependence  for  HF -N~ 
is  about  the  same  but  with  y  10%  larger.  The  values  for  HF-He  also  have*" 
about  the  same  J  dependence,  but  with  a  maximum  y  of  about  0.  02  to 
0.  03  cm"  Vatm  for  J=  1.  For  higher  J  state  transitions,  the  y  for  HF-He 
value  is  probably  about  0.  01  to  0.  02  cm"  ^/a±m. 


Therefore,  the  total  line  width  broadened  by  combined 
collisional  effects  of  HF-HF,  HF-H^,  HF-He,  and  HF-F2  will  probably 
be  around  0.  03  .  01  cm”  for  low  J  state  transitions  and  0.  015  +_  .  0005  cm* 

for  high  J  transitions.  Most  of  the  broadening  is  due  to  HF-He  collisions. 
As  the  temperature  of  the  gases  increase,  the  difference  in  y  between  high 
and  low  J  state  transitions  become  smaller,  and  the  average  line  width 
gets  broader,  which  decreases  the  gain.  Figure  35  shows  the  calculated 
output  energy  as  a  functinn  of  line  width  Av  .  For  simplicity  a  line  width 
average  over  all  J  state  transitions  was  used.  It  was  found  that  for  an 
averaged  line  width  of  0.  02  cm"  an  output  energy  of  3.  5  joule/100  cc  was 
calculated  for  90%  output  coupling  for  a  0.04/0.  04/0.  92  mixture  which 
is  in  rough  agreement  to  the  experimental  observations  (see  Fig.  16). 


1.  3.  6.  2  The  Effect  of  Rotational  Relaxation 


As  we  have  discussed  earlier  the  assumption  of  a  Boltzmann 
distribution  among  rotational  energy  levels  is  quite  questionable.  But 
with  present  lack  of  knowledge  of  rotational  relaxation  rates  in  HF,  it  is 
not  possible  to  make  a  detailed  calculation  of  its  effects.  Since  the  effect 
of  rotational  distribution  on  the  laser  performance  is  not  obvious  from 
the  equations,  we  shall  compare  the  calculated  output  energy  and  laser 
duration  only.  Figure  36  shows  the  calculated  laser  output  power  of  both 
frozen  rotational  and  equilibrated  rotational  distributions  for  the  standard 
gas  mixture  of  0.  04/0.  04/0.  92  and  output  coupling  of  90%.  It  is  clear 
that  for  the  case  of  the  frozen  rotational  distribution,  the  program  predicts 
a  much  larger  output  energy  and  longer  pulse  duration  than  that  of 
equilibrated  rotational  distributions.  Since  the  program  allows 
only  one  vibration- rotational  lasing  transition  on  each  vibrational 
energy  level,  the  frozen  rotation  case  still  maintains  its 
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original  rotational  distribution,  but  the  equilibrated  rotation  case  gives  a 
peak  intensity  spectra  close  to  that  of  a  Boltzmann  distribution.  The 
experimental  observation  falls  right  in  between.  Further  discussions  of 
these  rotational  effects  is  given  in  Section  1.4.2. 

1 .  3 .  6 .  3  The  Effect  of  Vibrational  Self  Relaxation  of  HF 

As  we  have  previously  noted  one  of  the  major  deactivation 
processes  in  the  HF  laser  system  is  the  vibrational  self  relaxation  of 
HF.  Airey  and  Smith(*-°)  have  recently  reported  their  experimental 
results  on  deactivation  of  HF  (v)  by  HF.  Table  IV  shows  the  data.  As 
the  authors  have  pointed  out,  for  v=2  the  measured  rate  constant 
5.5  x  105  sec"l  torr"!  is  believed  to  be  composed  of  contributions  from 
both  v-*R.  T  and  v-*v  rea  Tons: 


HF  (v=2)  +  HF  (v— o )  -  ~*R'  T>  I1F  (v  =  l )  +  HF  (v=o) 

,  v— R  .  T  .  Q  5  -1  -1 

=18x10  sec  torr 


(39) 


and 


HF  (v=2 )  +  HF  (v=o) 


2HF  (v=l) 


k^j  V"*’  =  3.7  x  10^  sec  ^  torr’* 


(40) 


Thus,  the  vibration-to -vibrational  energy  exchange  rate  is  only  two  times 
faster  than  that  of  vibration-to-translational  or  rotational  energy  transfer 
rates.  The  rate  constant  for  the  process: 


HF  ( v  —  1 )  +  HF  (v=o)  2HF  (v=o) 
,  v-*T .  R 


TO 


=  1.0+0.3x10^  sec  *  torr  * 


(41) 


is  about  two  times  smaller  than  that  of  k 


v-T.  R 
21 


If  one  assumes  that 


v_*t  R 

then  the  contribution  from  k  '  will  become 


.  v-*T.  R  ^  ,v-»T .  R 

k  °c  v  k. 

n,  n-1  10  v^T 

the  dominate  process  for  the  v=4  or  5  states.  If  k  is  indeed 

r  n,  n  -  1 


v ^  T  *1  A  J 

proportional  to  v,  then  kY.  '  should  be  5  x  1  x  1 0  ^  =  5x10^  sec  torr'*  and 

v^T  j  ^  5  1  1 

k^j  ‘  should  be  4  x  1  x  10  =4x10  sec  torr  ,  both  of  which  are 

considerably  smaller  than  the  observed  combined  rate  constants  (see 

Table  IV).  However,  if  one  assumes  a  v  ^  dependence,  then  k^  R*  = 

25x  1  x  10  =  2.5  x  lO^and  k^^'  R  =  16  x  1  x  10^  =  1.6x10^,  both  value  s 

of  which  agree  well  with  the  data.  Table  V  shows  the  effects  of  both  v—v 
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TABLE  IV 

RATE  COEFFICIENTS  FOR  DEACTIVATION  OF  HF  (v)  BY  HF 


,  HF,  HF 

k  . 

V,  v- 1 

(uncorr) 

V 

(torr"1. 
sec  "  1) 

z  HF,  HF 

V,  v-  1 

1 

9.  4  X  104 

70 

2 

5.  5  x  105 

12 

3 

6  x  10'1 

11 

4 

>  1.  7  x  I0b 

<  4 

5 

>  2.  6  x  106 

<  3 

6 

? 

? 

7 

? 

? 

8 

? 

? 
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and  v—T.R  rates  on  the  calculated  output  energy  and  laser  duration.  It 
clearly  indicates  that  the  v—v  rate  does  not  have  too  much  influence  on  the 
laser  output  energy  but  the  v-*T,R  rate  does.  Computations  using 

k  *  .  oc  v  predicts  twice  as  much  output  energy  as  that  using  kV~*  \  oc 

n,  n- 1  r  r  6/  b  n  n_i 

The  former  rates  also  predict  a  larger  laser  duration  than  that  using 

an  n^  dependence. 


1.3.  6.4  The  Effect  of  Chain  Branching  Reactions 

Russian  authors  ^  have  claimed  that  chain  branching  occurs 
through  the  reaction 


HF  (v  >  4)  +  F2 - »  HF  +  2F 


(VI) 


If  this  indeed  occurs  in  the  H2/F2  laser  system,  the  rate  of  the  overall 
reaction  will  increase  sharply  (exponentially),  and  the  duration  of  the  laser 
pulse  should  also  decrease  sharply  due  to  the  rise  in  gas  temperature. 
Unfortunately,  the  rate  constants  for  the  above  reactions  are  not  known 
with  any  degree  of  certainty.  If  one  adopts  rate  constants  of 

k  (v  <  4)  =  10"17  T1/2  exp  (E  -E)/RT  (42) 


and 


k  (v  >4)  =  10'17  T1/2  (43) 

where  E  =  37  kcal/mole  and  Ev  =  vibrational  energy  of  the  v*h  state  as 
claimed  by  the  Russian  authors,  the  contribution  of  VI  to  the  overall 
reaction  becomes  unimportant.  However,  if  one  assumes  a  rate  constant 
of  1 0 ~  1  2f  the  program  predicts  an  output  energy  twice  as  large  and  a 
pulse  duration  of  about  two  times  shorter  than  that  using  the  slower  rate. 
The  calculation  with  the  fast  chain  branching  reaction  does  give  a  result 
(laser  duration  and  output  energy)  in  better  agreement  with  the  data 
(Table  V),  but  cannot  explain  the  variation  of  output  energy  with  coupling. 

1 .  3.  6.  5  The  Assumption  of  Constant  Cavity  Loss  and  Gain- 
Equal-Loss  Condition 

As  we  have  stated  at  the  very  beginning,  the  modeling 
calculation  assumes  that,  during  the  laser  action,  gain  equals  loss  and 
that  the  optical  loss  is  a  constant.  From  the  light  scattering  experiment, 
we  have  clearly  demonstrated  that  the  beam  can  te  seriously  distorted 
inside  the  laser  cavity.  Therefore,  the  assumption  of  constant  loss  during 
the  laser  duration  could  be  questioned.  In  order  to  test  the  effect  of 
possible  beam  distortion  on  the  performance  of  the  HF  laser,  a  qualitative 
analysis  was  made. 
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where 

gQ  is  the  small  signal  gain, 

Lt  is  the  total  losses  (L^  =  1  -  (1  -  L>c) ( 1  -  LD))» 

I  is  the  intensity  in  the  limit  L,j,— Lc— -0, 

Lc  is  the  output  coupling, 

£  is  the  cavity  length, 

Ln  (=1-R)  is  the  dissipative  losses  or  in  our  case  possible  losses 
caused  by  beam  distortion. 

Figure  37  shows  i/I  vs  LD  for  various  values  of  small  signal 
gain  for  two  experimental  values  of  output  coupling.  The  fact  that  twice  as 
much  intensity  is  observed  experimentally  for  the  highei  output  coupling 
case  indicates  that  L„  >  .4  since  l/lQ  must  certainly  be  less  than  one 
half  for  the  lower  coupling  case.  In  any  event,  this  analysis  has  clearly 
indicated  that  for  high  gain  lasers,  if  the  beam  becomes  distorted  due 
to  density  variations  within  the  cavity,  the  laser  energy  extracted  from 
a  high  output  coupling  cavity  is  not  necessarily  smaller  than  that  of 
smaller  output  coupling.  This  conclusion  is  consistent  with  our  experi¬ 
mental  observations,  but  does  not  prove  that  beam  distortion  (or  a  high 
value  of  Lp)  is  a  dominant  effect  in  our  laser. 
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1.4  SUMMARY  ON  THE  PERFORMANCE  OF  FLASH  LAMP  INITIATED, 
PULSED  HF  CHEMICAL  LASER 

1.4.1  On  the  Chemical  Efficiency 

In  order  to  evaluate  an  upper  limit  on  the  amount  of  chemical 
energy  that  is  available  for  laser  action,  the  initial  energy  distributions  of 
both  the  F  +  H2  and  H  +  F2  reactions  have  to  be  carefully  studied.  Re¬ 
cently  Polanyi  et  al  (  ^  )  have  given  detailed  reports  on  the  experimentally 
measured  vibrational/rotational  energy  distributions  of  these  two  reactions. 
Their  results  are  reproduced  in  Figs.  38  and  39.  They  found  that  the  cold 
reaction  (F  +  H2)  favors  formation  of  the  v  =  2  state  and  is  capable  of  pro¬ 
ducing  vibrational  states  of  HF  up  to  v  =  3.  For  the  hot  reaction,  (H  +  F2) 
they  found  that  the  reaction  favors  formation  of  the  v  =  6  state  and  is  capable 
of  producing  vibrational  states  of  HF  up  to  v  =  8.  Besides  the  formation  of 
vibrational  excited  HF,  these  data  also  indicate  significant  rotational 
non- Boltzmann  distributions  as  well.  Particularly,  the  cold  reaction 
preferentially  populates  higher  rotational  states  as  the  vibrational  state 
population  decreases,  while  for  the  hot  reaction  all  vibrational  states  have 
about  the  same  rotational  distribution  peaked  at  about  J  =  8. 

Using  this  data,  the  chemical  efficiency  of  our  laser  system  can 
be  analyzed  as  follows: 

TOTAL  CHEMICAL  ENERGY  LIBERATED  FROM  THE  REACTION  OF 


F  +  H2  -  HF  +  H  -  AH  =  31.  5  kcal/mole  0.66  0.08  0.26 

H  +  F2  -  HF  +  F  -  AH  =  102  kcal/mole  0.53  0.03  0.44 

H2  +  F2  -*  2HF  -  AH  =  133.  5  kcal  or  66.  75  kcal/mole  of  HF  formation 

where  fy,  fj^,  and  Tp  are  the  fraction  of  AH  that  end  up  in  the  HF  vibra¬ 
tional,  rotational,  or  translational  mode  respectively. 

TOTAL  ENERGY  INTO  VIBRATION  =  20.8  +  54.  0  =  74.  8  kcal  or  37.4  kcal/ 

mole  of  HF  formation 

=  1290  joule/100  CC  of  Hp  +  F2  at  1  atm 
or  64  5  joule/100  CC  of  HF  formation 

TOTAL  VIBRATIONAL  ENERGY  AVAILABLE  IN  100  CC  MIXTURE  OF 
H2/F2/He  =  0.  04/0.  04/0.  92 

=  1290  x  0.  04  =  51. 6  joule/100  CC  of 
standard  H2/F2/He  mixture  at  1  atm 
or  51.6  joule/8  CC  of  HF  formation 
at  1  atm. 
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RELATIVE  POPULATION  N(V)  RELATIVE  POPULATION  N(J) 


TOTAL  LASER  ENERGY  OBSERVED  =  4  JOULE  IN  3  /iSEC 
(LASER  ENERGY  (/(AVAILABLE  VIBRATIONAL  ENERGY)  «  8% 
EXTENT  OF  CHEMICAL  REACTION  ESTIMATED  WITH  F/F2  =  0.  01 
FOR  3  /iSEC  PERIOD  as 30% 

TOTAL  CHEMICAL  ENERGY  LIBERATED  INTO  VIBRATION  WITHIN 


3  /iSEC  DURATION  =  15  JOULES 

PERCENT  OF  VIBRATIONAL  ENERGY  CONVERTED  INTO  LASER 
ENERGY  DURING  THE  3  /iSEC  PERIOD  =  4/15  ss30% 

THE  MAXIMUM  POSSIBLE  CHEMICAL  EFFICIENCY  =  59% 

CHEMICAL  EFFICIENCY  =  LASER  ENERGY/TOTAL  CHEMICAL  ENERGY 


RELEASED  a  4% 


Therefore,  the  chemical  efficiency  of  our  system  for  0.04/ 

0  04/0.92  mixture  is  4%.  The  chemical  efficiency  for  the  0.04/0.02/0.94 
mixture  and  the  0.  16/0.  08/0.  76  mixture  are  also  4%.  A  substantial  fraction 
of  the  available  energy  of  the  reaction  that  goes  into  vibration  has  been 
recovered  as  laser  energy.  The  major  reason  for  the  termination  of  laser 
action  is  probaoly  the  heating  of  the  gas  mixture  as  the  reaction  proceeds 
with  subsequent  loss  of  gain. 


The  overall  electric-to-laser  energy  efficiency  obtained  from 
a  mixture  of  0.  16/0.  08/0.  76  was  calculated  to  be  about  1.  3%  which  is  about 
twice  as  much  as  that  reported  in  our  previous  contract  report.  -1  '  Since 
the  active  volume  of  our  laser  device  is  only  0.  1  liter  at  a  pressure  of 
1.  1  atm,  the  final  laser  energy  density  in  units  of  joule  s/liter -atm  is 
80  ioule/liter-atm  with  a  composition  of  16%,  H2,  8%  F2  and  76/>  He.  For 
higher  concentrations  of  F2  and  H2  in  He,  the  output  energy  would  be 
expected  to  increase  proportionately. 


\.4'2  On  the  Assumption  of  Equilibration  of  Rotational  Energy  Levels 

Among  almost  all  the  known  vibration- rotational  laser  systems 
reported  in  the  literature,  the  equilibration  of  the  rotational  energy  levels 
has  always  been  assumed  to  have  been  established  before  the  laser  transi¬ 
tions.  This  assumption  is  probably  true  for  lasing  molecules  having  small 
rotational  energy  level  spacing  and  hence  faster  rates  of  rotational  relaxa-  { 
tions  than  HF.  Since  the  molecule  HF  has  a  rotational  constant  B  -20  cm  , 
the  spacing  between  J  =  10  and  J  =  9  levels  is  2BJ  =  400  cm'1  which  is  about 
twice  kT  at  room  temperature.  The  time  for  the  rotational  relaxation  of 
HF  (J  =  10)  state  could  be  very  long  and  take  many  collisions  to  reach 
equilibration.  Polanyi  et  al  (18)  have  also  made  an  experimental  estimate 
of  the  probability  of  rotational  relaxation  of  HC£.  The  results  are  replottec, 
in  Fig.  40  as  the  probability  per  collision  of  a  rotational  transition  as  a 
function  of  AE,  the  energy  spac:rg  between  rotational  states.  It  was  found 
that  for  low  J  states  AJ  =  +  2  and  AJ  =  +  1  transitions  are  all  probable  But 
for  high  J  levels,  the  AJ  =  4_  1  transitions  dominate.  If  one  assumes  that 
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HF  ,r^!;atl0nal  encrsy  levels  relax  with  a  similar  energy  dependence  (E,  -  r. , 
as  II C  f  and  tKnrofnvn  fUn  «  t..  j _ i  i  i  i  .  j.  .  ^  .1 


E 


- -  i  tidA  wiui  d  similar  energy  aepenaence  [v^i  -  J 

ns  C(  and  therefore  the  probability  depends  only  on  the  rotational  energy 
spacing,  then  the  -elaxation  of  the  J  =  10  state  should  require  100  collisions, 
i  he  particular  J  state  of  HF  corresponding  to  AE  has  been  indicated  at  the 
top  of  Fig.  40.  It  is  evident  that  100  collisions  is  quite  slow  in  comparison 
with  the  rate  of  stimulated  emission  and  even  slower  than  the  rate  of  vibr,  - 
tional  relaxation  of  HF. 


Thus  the  assumption  of  a  Boltzmann  distribution  among  all  the 
rotational  energy  levels  before  the  lasing  accion  is  indeed  questionable.  The 
fact  that  (i)  far-IR  laser  emission  was  observed  in  our  H^/F,  system  and 
(ii)  two  distinctive  pulses  within  each  individual  vibrational- rotational 
transition  are  observed  for  high  J  states,  strongly  supports  the  argument 
that<  at  the  beginning  of  HF  laser  action,  the  rotational  energy  equilibration 

must  have  not  been  established.  Unfortunately,  the  double  resonance - 

absorption  experiment  to  measure  the  HF  rotational  relaxation  rate  was  not 
finished  and  was  inconclusive.  With  the  existing  data  we  can  only  state  that 
the  rates  of  rotational  relaxation  of  low  J  states  are  relatively  fast.  For  a 
detailed  knowledge  of  the  mechanism  and  rate,  more  experimental  work  is 
needed. 


Again,  it  is  not  clear  what  effect  the  rotational  energy  non¬ 
equilibration  has  on  extracting  IR  laser  photons.  For  one  thing,  it  helps 
in  building  up  population  invei  sion  among  certain  levels.  It  also  decreases 
the  tendency  of  losing  laser  gain  through  the  process  of  rotational  relaxation. 
The  fact  that  (i)  for  the  same  gas  mixture,  different  output  coupling  windows* 
(with  the  same  optical  property  at  2.  7  jli  but  considerably  different  behavior 
in  the  far  IR  region)  extract  different  amounts  of  IR  laser  energy,  and  (ii) 
adding  SF^  to  the  H2/F2  mixture  significantly  decreases  the  amount  of  IR 
laser  output  energy,  also  supports  the  contention  that  the  rotational  energy 
distribution  is  playing  an  important  part  in  the  process  of  HF  IR  laser  action. 

1.4.3  On  the  Possibility  of  Parasitic  Oscillation 

Scientists  at  the  Aerospace  Corporation  (19)  have  claimed  that 
parasitic  oscillation  may  be  playing  an  important  role  in  their  HF  laser. 

Such  oscillations,  if  consisting  of  modes  that  bounce  off  the  walls  of  the 
laser  tube,  could  result  in  the  output  laser  be  urn  appearing  at  larger  angles 
to  the  axis,  and  thereby  possessing  poor  optical  quality.  In  some  respects 
our  laser  performance  is  consistent  with  their  arguments.  For  example, 
we  get  only  half  as  much  energy  out  by  using  a  50%  output  coupling  (Si) 
window  compared  to  a  90%  output  coupling  (Sapphire).  However,  this 
phenomenon  could  also  be  attributed  to  the  differences  in  optical  properties 
o  these  two  output  window  materials.  However,  the  output  energy  observed 
by  using  a  NaCl  window  is  roughly  the  same  as  that  obtained  by  using  Si 
flat  although  these  materials  have  different  %  transmission.  Furthermore, 
the  output  energy  obtained  from  KBr  output  coupling  window  is  less  than 
that  of  NaCl  flat  although  these  materials  have  almost  identical  transmis¬ 
sivity.  Therefore,  the  parasitic  argument  is  far  from  conclusive.  Never¬ 
theless,  the  close  agreement  between  the  burn  pattern  and  the  rough 
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calculation  of  spot  size  and  multimode  pattern  strongly  indicates  that 
parasitic  oscillations  are  not  dominating  in  our  HF  laser.  Of  course, 
parasitic  oscillations  are  still  an  important  consideration  in  a  intrinsically 
high  gain  laser  such  as  this  ^2^1  system. 

1,4,4  On  the  Validity  of  Model  Calculations 

The  rough  agreement  between  the  experimental  observations 
(such  as  output  energy  and  pulse  duration)  and  predictions  (Table  V)  re¬ 
sulting  from  the  theoretical  calculation  have  indicated  that  the  model  used 
is  a  reasonable  one.  At  least,  the  model  is  capable  of  giving  a  qualitative 
picture  of  the  performance  of  the  HF  laser  studied.  However,  there  are  a 
few  fundamental  drawbacks  to  the  model  that  could  be  the  basis  for  further 
research: 

(a)  The  model  is  constructed  based  on  the  assumption  that  a 
constant  cavity  loss  is  maintained  during  the  laser  action, 
which  may  not  be  true. 

(b)  The  model  does  not  take  care  of  the  effect  of  rotational 
relaxation  and  the  possibility  of  vibration-to- rotational 
energy  transfer  which  could  seriously  change  the  output 
spectra  and  energy. 

(c)  There  is  a  contraction  in  logic  between  the  constant  gain 
assumption  and  the  rapidity  of  relaxation  rates.  Since 
the  constant  gain  assumption  is  constructed  on  the  basis 
that  any  of  the  relaxation  times  should  be  long  compared 
with  the  transient  time  of  photons  inside  the  cavity,  it 
clearly  requires  that  any  kind  of  relaxation  process  must 
be  slow  and  take  a  long  time  to  reach  equilibrium  if  it  is 
not  produced  with  an  equilibrated  distribution.  The  results 
from  Polanyi' s  experiments  have  shown  that  the  initial 
rotational  distribution  is  highly  non- Boltzmann.  Thus,  if 
one  uses  a  Boltzmann  distribution  of  rotational  energy 
levels,  one  must  assume  th.' t  the  rate  of  rotation  relaxation 
is  much  faster  than  the  transient  time  of  photons  in  the 
cavity.  Thus,  from  a  logical  point- of- view,  the  usefulness 
of  constant  gain  assumption  is  quite  questionable  for  this 
type  of  lase  r. 

(d)  For  the  gain-equals-loss  condition  to  be  true,  the  beam 
inside  cavity  must  be  regenerative,  which  may  not  be 
true  in  our  HF  experiments.  However,  little  work  was 
done  under  this  program  to  try  to  improve  the  fluid  mixing 
and  uniformity.  This  problem  appears  to  be  solvable 
experimentally  and,  hence,  need  not  be  a  fundamental 
limitation. 


1.  5  THEORETICAL  CALCULATION  ON  THE  VIBRATIONAL  ENERGY 
TRANSFER  RATE  CONSTANTS 

1.5.1  Introduction 


It  has  been  a  long  goal  of  this  program  to  develop  fast  yet  accu¬ 
rate  computer  programs  which  contain  the  essential  physics  for  the  calcula¬ 
tion  of  the  vibration- relaxation  cross  section  in  an  atom -molecule  collision. 
These  programs  in  a  simplified  version  have  previously  been  applied  to  the 
V-R  relaxation  of  OH  by  argon.  A  more  realistic  version  of  these  programs 
is  being  developed  and  described  in  the  following  section  (1.  5.  2).  It  is 
planned  that  such  theory  will  be  applied  to  problems  of  interest  in  the  HF 
laser  in  the  immediate  future. 

1  5.  2  Cross-Section  Theory 

The  number  of  collisional  transitions  per  second  or  the  transition 
probability  P.^  from  state  i  to  state  f  in  the  Born  approximation  is  given  by 


P 


B  lit 
if  Y 


<f  I  V  I  i  > 


2  6  (Ef  -E.) 


(45) 


where  V  is  the  interaction  potential  between  collision  partners  A  and  B;  E. 
and  Ef  are  initial  and  final  energies;  Ef  =  +  Eg1  +  E^1  +  Efj1,  t  and  i 

denote  the  translational  and  internal  energies  and  the  integration  is  over 
final  wave  vectors  of  translational  motion.  We  express  initial  and  final 
wave  functions  as 


i  =  (2tt) 


-3 


&A 


£a> 


*-B 


-B) 


B 


(46) 
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and 


-3  i(^A 
i  =  (lit)  e 


e 


-B 


X' 


B 


(47) 


where  Y  are  the  wave  functions  of  the  internal  motion,  r  .  and  r_  are  the 
center  of  mass  coordinates  of  A  and  B,  and  the  final  state  is  denoted  by 
primes.  We  now  assume  that  the  interaction  potential  between  A  and  B  is 
the  sum  of  the  interaction  potential  of  the  constituent  atoms,  i.  e. 


V  =  N  ^  V  (r  ) 

/Lt 


(48) 


where  the  summation  is  over  all  atoms  fl  constituting  A  and  all  atoms  v  con¬ 
stituting  B  and  we  take  the  potential  between  the  constituent  species  to  be 
a  central  one.  In  addition  we  assume  that  matrix  tlement  in  Eq.  (45)  is  the 
algebraic  sum  of  such  elements  from  individual  scattering  centers.  This 
enables  us  to  write 


<  f  |  V  |  i  > 


<t|v|(v’ j> 

M.  v 


(48a) 


Now  we  introduce  the  Fourier  transform  of  the  wave  functions  for  A  and  B 


(2 


7rf3/2  y.  exp[ik.  •  rj]  s  =  (2tt)"3'  2  J  dk^  exp[ikv  •  r  ^  ] 


X4>  lLy-Lv_v  *  *  > 


(49) 


where  i  =  A,  B.  Clearly  (r_  .  .  .  £_  k  v_j,  .  .  )  is  the  amplitude 

that  atom  v  has  momentum  k  simultaneous  with  other  atoms  being  at  r  ., 

—  v  —  1 

•  •  • .  £v_i<  £v+r  ‘  ‘  ‘  ' 
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The  matrix  elements  <  f 
can  be  written  as 


needed  to  evaluate 
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(50) 


where  dr  =  II  dr  and  dr  _  =  n  dr  . 

^  H-  H-  Li  V  —  v 

,  N*°W  WG  introduce  center  of  mass  and  relative  position  and 
vectors  for  atoms  /j.  and  v  by  the  relations 
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m  k  -  m  k 
k  =  v  ~U  u-v 

—  tiv  m  +  m 

^  [1  v 


This  enables  us  to  write 
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(53) 


where  prime  in  dT '  .  means  exclusion  of  dr  from  the  volume  element. 
r  A  -  ji 

At  this  point  we  restrict  our  considerations  to  the  case  where  A 
is  an  atom.  In  keeping  with  the  spirit  of  our  study  we  take  the  atom  to  be  a 
structureless  point  particle.  For  this  (from  Eq.  (49))  case  v.  =  1  and 

<r\ 


$.=S(k  -k  ) 
—  A  ~n' 


(54a) 


and 


Substituting  Eqs.  (54)  in  (53) 
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Assume  that  the  integral 


V  (rav)  exp 


a  m  +  m  ^av^  I 
a  v 


is  slowly  varying  function  of  /  and  using  the  relation 
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where  q  k'a  -  k&  is  the  momentum  gained  by  the  incident  atom.  Writing 


L  =  r  n  ■>  r 
v  — B  — , 


where  r_B  is  the  distance  of  CM  of  B  from  our  space  fixed  coordinate 
system  and  is  the  distance  of  atom  v  from  CM  of  B,  we  get 


-89- 
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The  delta  function  in  Eq.  (59)  is  simply  the  result  of  conservation  of  linear 
momentum.  Equation  (55)  can  row  be  rewritten  as 


,-3 
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Substitution  of  this  result  in  Eq.  (45)  gives 
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We  will  evaluate  Eq.  (64)  under  the  approximation 


K  =  k 
—  av  —a 


k  =*  k  +  k  n 
—  v  —a  — B 


(65) 


i.  e.  ,  internal  degrees  of  freedom  make  no  contribution  to  K  .  We  can 
then  write  av 
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(66) 


where  k  is  the  initial  relative  momentum  and  K  is  the  momentum  of  center 
of  mass  of  atom  and  molecule.  The  approximation  introduced  in  Eq.  (65) 
is  not  fundamental  to  the  formulation  presented  here.  To  simplify  the  cal¬ 
culations  further  we  approximate  the  left-hand  side  of  Eq.  (66)  by  k^,  i.  e.  , 
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‘( 


m  +  m 
a  v 


- T - W 

a  B  / 


In  addition  using  the  relation 


(67) 
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we  get 
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IMHHSi 


where  q  +  =  k'  +  k,  £  ~ (2g  +  1 ),  is  Legend  re  Polynomial  and  is  the 

sphericaT  Bessel  function.  Equation  (69)  uses  a  straight  line  trajectory 
during  a  collision.  Now  we  switch  to  distorted  wave  approximation  and 

,  iOfOO 

replace  (kr.^)  by  (kr^)  Rg  (kr,^)  e  where  is  the  solution  of 

the  equation 
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and  (k)  is  the  asymptotic  phase  shift  for  partial  wave  C  and  energy  k  in 
potential  V  (r).  Appendix  I  discusses  a  method  for  fast  numerical  evaluation 

OO 

of  the  integral  f  R^(kr)  R^(k'r)  V  (r)  dr  for  the  case  when  V  (r)  is  Morse 
*1) 

potential,  i.  e.  , 
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where  D  is  the  wall-depth  at  RQ  and  P  is  the  range  parameter. 
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Using  the  expression  for  <  X  p  e 


Xg  >  derived  in 


Appendix  II,  we  get  the  distorted  wave  cross  section  a.f 
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The  computation  of  cross  section  using  Eq.  (70)  is  being  carried 
out  first  for  atom  molecule  situation  and  it  is  later  proposed  to  extend  these 
computations  to  molecule-molecule  situations. 

1.5.3  Conclusion 

At  present  the  theory  is  still  at  the  developing  stage.  The 
research  will  be  continued  under  the  support  of  ARPA  and  AMICOM  in  Contract 
No.  DAH01  -72-C-0995. 

In  the  future  the  cross  sections  for  inelastic  non- reactive  scat¬ 
tering  processes  will  be  calculated.  The  theory  will  eventually  be  applied 
to  reactions  between  HF  molecules.  Reactions  like: 
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HF  (v  =  n)  +  HF  (v  =  o)  -  HF  (v  -  n  -  1)  +  HF  (v  =  1) 


and 


HF  (v  =  n)  +  HF  (v  =  o)  -*  HF  (v  =  n  -  1)  +  HF  (v  =  o) 

will  be  investigated.  It  is  hoped  that  this  work  will  determine  which  are 
the  most  important  deactivation  processes  in  the  HF  laser.  Since  (1) 
these  types  of  reactions  have  been  found  to  be  the  most  important  relaxation 
processes  involved  in  HF  laser  system  and  (2)  the  measurements  on  these 
types  of  rate  constants  involving  higher  vibrational  levels  are  quite  difficult, 
the  theory  of  this  kind  should  be  of  extreme  value  in  furthering  our  under¬ 
standing  of  the  HF  laser. 
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1  6  PULSED  HF  CHEMICAL  LASER  GENERATED  BY  OTHER  CHEMICAL 
REACTION  AND/OR  INITIATION  TECHNIQUES 

16.1  Motivation 

Our  theoretical  and  experimental  studies  on  the  performance  of 
the  HF  laser  generated  by  the  flash  initiated  H2/F2  reaction  have  clearly 
drnionst’-'.ted  that  the  HF  chemical  laser  may  offer  considerable  potential 
for  future  high  power  military  applications.  The  laser  energy  obtained  in 
the  H2  '¥ 2./  He  mixture;  80  joule/f-atm  indeed  represents  the  current  state- 
of-art  among  all  the  chemical  laser  systems  known  to  date.  However,  here 
are  a  few  drawbacks  that  are  worth  noting: 

( 1)  Due  to  the  reactivity  of  the  H  ,/f  5  mixture ,  the 
gases  cannot  be  premixed  and  stored  in  room 
temperature  containers. 

(2)  The  overall  electrical  efficiency  is  seriously  limited 
by  the  ineffectiveness  of  converting  electrical  energy 
into  useful  ultraviolet  photons  for  dissociating  F£  to 
initiate  the  reaction. 

In  order  to  overcome  these  two  problems,  two  separate  approaches 

were  taken: 

(1)  Other  energetic  chemical  reactions  (besides 
H2/F2  or  D2/F2)  were  studied.  These  other 
reagents  are  in  most  cases  quite  stable  and 
easier  to  handle  or  to  premix  compared  to  H2/F2. 

(2)  Other  techniques  that  might  directly  dissociate  F2 
molecules  without  going  through  the  route  of  gener¬ 
ating  ultraviolet  photons  were  investigated,  namely 
direct  electric  discharges  or  high  energy  electron 
beam  initiation. 


1.6.2  New  Possible  Cha  in  -  B  ra  nching  Chemical  Laser  Systems 
1 .  6.  2.  1  Introduction 

Manv  Soviet  laser  chemists  have  proposed  "chain-branching" 
chemical  lasers'^  )  and  have  claimed  successful  operation  of  these  lasers.  (  ) 

However,  only  minor  efforts  have  been  made  by  scientists  in  the  United  States. 
Krogh  and  Pimenlel^O)  recently  reported  a  study  on  possible  cha  in  -  branching 
(C-B)  laser  action  originating  from  chemical  reactions  in  a  C1F3,  H2/Ar  mix¬ 
ture.  Gensel,  Kompa  and  Wanner(^l)  have  also  observed  C-B  laser  action 
from  a  mixture  of  IF5  and  H2.  Thus,  it  is  reasonable  to  reinvestigate  these 
systems  and  other  similar  reactions  that  might  lead  to  the  generation  of 
efficient  C-B  chemical  lasers. 


From  the  chemical  kinet'cs  point  of  view,  there  are  several  ad¬ 
vantages  in  the  performance  of  a  C-B  chemical  laser: 

(1)  In  the  C-B  system,  the  population  inversion  maxi¬ 
mum  is  only  weakly  dependent  on  the  initiation 
energy. 

(2)  The  reaction  rate  grows  exponentially  with  time. 

Advantage  (  l)  would  indie  ite  a  possible  characteristic  of  high  efficiency  and 
advantage  (2)  would  short;n  the  laser  pulse  and  raise  the  laser  power.  For 
most  chemical  l?ser  sys.ems  studied,  the  most  competitive  process  reducing 
laser  emission  id  the  relaxation  process.  A  faster  chemical  reaction  would 
definitely  favor  a  higher  percentage  of  chemical  energy  converted  into  laser 
photons. 

In  the  following  section  the  performance  of  chemical  lasers  gen¬ 
erated  from  the  reactions  between  hydrogen  and  halogen  fluorides  (such  as 
CIF3,  BrFci,  IFc^ ,  and  CfF)  are  presented.  This  type  of  reaction  is  of  ex¬ 
treme  interest  for  many  reasons: 

(1)  It  is  one  of  th~  few  systems  that  might  proceed  by 
either  a  chain  or  chain-branching  (C-B)  reaction, 
which  enhances  the  potentiality  of  future  develop¬ 
ment  as  a  possible  high-power,  repetitively  pulsed 
chemical  laser. 

(2)  The  overall  reaction  is  very  exotherr.  ic  (as  compared 
with  Hz  +  F 2  system)  and  capable  of  producing  both 
vibrationa lly  and  rotationally  excited  species  which 

is  desirable  as  far  as  the  production  of  the  popula¬ 
tion  inversion  is  concerned. 

(3)  Despite  (1)  and  (2)  described  above,  the  reagents 
are  quite  stable  as  compared  to  mixtures  of  Hz/F z- 
Although  the  chemicals  will  react  slowly  when  mixed 
the  reaction  rates  are  slow  enough  for  a  complete 
mixing  in  a  regular  slow  flow  system.  By  adding  a 
small  amount  of  inhibitor  like  O2,  the  stability  has 
been  further  improved. 

(4)  Being  produced  by  chain  or  C-B  reactions,  the  output 
laser  energy  thould  not  be  very  sensitive  to  the  ini¬ 
tiation  energy,  which  suggests  that  in  principle  only 

a  relatively  small  triggering  energy  is  required  for 
the  laser  ignition. 

Table  VI  lists  the  important  thermodynamic  data,  physical  and 
chemical  properties  of  the  halogen  fluorides  that  we  have  studied.  Note  that 
all  of  the  chemicals  are  commercially  available  and  can  be  stored  at  room 
temperature. 
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TABLE  VI 

IMPORTANT  CHEMICAL  AND  PHYSICAL  PROPERTIES  OF  XFr 


Property 

ci  f3 

Br  F5 

^5 

SF6 

N2  F4 

Melting  Point 

-76.  32° C 

-60. 5° 

9.43° 

-50. 5° 

— 

Boiling  Point 

11.75° 

41. 30° 

100. 5° 

-63. 8° 

-7  3° 

Vapor  pressure 
at  room  temp. 

1 . 65  atm 

0.  53  atm 

10  Torr 

20  atm 

6  atm 

Corrosion  on 
Copper  or  Steel 

very 

slightly 

/ery 

slightly 

no 

no 

no 

Etching  Pyrex 
or  Quartz 

no 

no 

no 

no 

no 

Automatic  igni¬ 
tion  in  H2  or 
hydrocarbons 

yes 

no,  but 
will  ig¬ 
nite  by 
a  flame 

no 

no 

no 

Information  obtained  from  Halogen  Chemistry,  Edited  by  V.  Gutmann, 
Academic  Press,  1967. 
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1.6.2 .2  Reaction  Scheme 

Following  the  flash  photolysis,  the  possible  reactions  involved 


tn  the  mixture  of  H^  and  Cf  F3  can  be 

postulated 

a  s : 

ClF3+hv—  CIF,  4  F 

AH 

60  kcal/mole 

(74a) 

—  C1F  +  2F 

AH  = 

64 

(74b) 

—  C1F  4  F2 

AH 

26 

(74c) 

-  F2  +  C l  4  F 

AH 

85 

( 74d) 

CIF3  absorbs  strongly  at  wavelengths  beyond  2600  ft,  and  quanta 
of  this  wavelength  carry  110  kcal/mole  of  energy;  thus  dissociation  of  CIF3 
could  proceed  by  any  of  the  above  mentioned  reactions.  The  fact  that  HF 
laser  action  was  observed  indicates  that  reaction  (74c)  is  not  important. 
Besides  reaction  (74c),  all  the  other  reactions  lead  to  the  production  of  one 
or  two  F  -  a  tom  s .  Subsequent  photoly  s  i  s  of  C1F,  of  C  1  F  2  and  of  F^  is  also 
pos  s ible . 

The  fluorine  atoms  will  initiate  the  pumping  reactions,  which 
produce  H-atoms  and  introduce  the  reaction  chain  sequence: 


r  +  h2  - 

HF  +  H 

AH  =  - 3  1 .  6  kca l/mole 

(75) 

H  +  C1F3  - 

HF  +  CIF? 

AH  -  -76 

(76a) 

— 

HF  +  C  IF  +  F 

AH  =  -72 

(76b) 

— 

HF  +  Cl  +  2F 

AH  -12 

(76c) 

— 

HF  +  C l  4  F2 

AH  =  -50 

( 76d) 

— 

HC1  +  F2  +  F 

AH  =  -17 

(76e) 

— 

HC1  +  3F 

AH  =  +20 

(76f) 

Except  for  reactions  (76e)  end  (76f),  which  generate  only  new  F -atoms,  all 
the  other  reactions  (76)  proceed  by  generating  both  vibrationa lly  excited  HF 
molecules  and  new  F-atoms,  Cl-atoms,  and  radicals.  Atoms  will  continue 
their  reaction  with  by  either  (75)  or 

C1F2  +  H2  —  HC1  +  H  AH  =  +  1  kcal/mole  (77) 

Although  reaction  (77)  is  thermoneutral,  it  will  produce  an  additional  H-atom 
for  reaction  (76)'  s.  The  CIF^  radical  will  probably  react  with  H ^  by 
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-  HF  +  C1F  4  H 

AH 

-28  kcal/mole 

(78a) 

—  2  HF  +  Cl 

AH 

-66 

(78b) 

—  HF  +  HCl  4  F 

AH  ^ 

-32 

( 7  8cl 

A  pain  in  either  reaction  (78)  both  vibrationally  excited  HF  molecules  and 
new  atoms  will  be  generated  and  recycled.  The  reaction  between  C1F2 
radical  and  H  would  add  more  F-atoms  by  processes; 


H  +  CIF^  -  HF  4  C1F 

AH 

-  1  32  kcal/mole 

(79a) 

—  HF  4  Cl  *  F 

AH 

7? 

(79b) 

—  HCl  4  F2 

AH 

-77 

(79c) 

—  HCl  4  2F 

AH 

-39 

(79d) 

Of  course,  C1F  would  also  react  with  H  by 

H  4  C1F  —  HF  4  Cl 

AH  = 

-43.  4  kcal/mole 

(80a) 

—  HCl  4  F 

AH  ^ 

-75.  9 

(80b) 

and  liberate  new  atoms. 

As  one  can  clearly  see,  the  production  of  vibrationally  excited 
HF  and  HCl  molecules  can  be  very  efficient  and  chain  propagation  may  occur 
much  faster  than  that  of  the  H2/F2  system  if  reactions  (76)  proceed  with  a 
reasonable  rate  (as  compared  with  reaction  (75),  which  we  know  is  fast). 

And  as  Professor  Pimentel  pointed  out,  insofar  as  reactions  (76),  (78)  or 
(79)  are  concerned,  they  contribute  constructive  chain-branching. 

1.6.  2.  3  Experiment  and  Results 

The  same  apparatus  (see  Fig.  1  )  and  operational  procedures 
identical  to  that  used  in  H2/F2  system  were  employed.  The  halogen  fluorides, 
ClF-i,  BrFc,  IFc,  or  C1F  were  used  in  place  of  F2.  Measurements  of  the 
laser  output  energy  were  made  with  the  Avco  calorimeter  and  the  laser  dura¬ 
tion  or  pulse  shapes  were  monitored  with  an  InSb  PEM  detector  viewing 
scattered  radiation  from  the  calorimeter  surface.  Figure  41  shows  a  typical 
oscillogram  of  a  laser  pulse  observed  in  a  H2/C1F3  mixture.  The  total  gas 
pressure  was  1.  1  atmospheres  with  4%  H2,  4%  CIF3  and  92  %  He.  The  laser 
pulse  duration  is  about  1/isec  and  the  output  energy  is  over  1  joule.  I  he 
mixture  of  H?/IF7  was  found  to  be  extremely  stable.  It  could  be  premixed 
and  flashed  many  times  to  get  laser  pulses.  Figure  42  shows  the  laser  pulse 
observed  in  a  H?/CIF  mixture.  The  total  gas  pressure  is  also  1.  1  atmo¬ 
spheres  with  4%  H2,  4%  C1F  and  92%  He.  The  laser  duration  is  about 
1.0/isec  and  the  output  energy  is  over  .7  joule.  Table  VII  lists  the  output 
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Fig. 


41  Oscillogram  of  the  Laser  Pulse  Observed  in  H^/Cf  F 
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TABLE  VII 


fIF  CHEMICAL  LASERS  GENERATED  BY  REACTIONS 
OTHER  THAN  THE  H2  +  F£  SYSTEM 


Output  Energy 

Duration 

1) 

H2 

+  C1F3 

1  joule 

-  1  f_L  sec 

2) 

H2 

+  BrF  3 

0.  3  joule 

-  1  fisec 

3) 

H2 

+  C1F 

0. 7  joule 

(both  HC1  and  HF  lase) 

-  10/isec 

4) 

H2 

+  BrF5 

0.  5  joule 

~  1  fjisec 

5) 

H2 

+  IF5 

0.  5  joule 

~  1  (isec 
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energies  and  pulse  durations  observed  In  various  H^/halogen  fluoride.  He 
mixtures  studied.  These  results  compare  favorably  with  our  early  experi¬ 
ments  on  the  H2/F2  system,  and  the  relatively  high  energy  obtained  also 
suggests  that  these  reactions  may  operate  on  a  chain  mechanism.  Future 
HF  laser  development  based  on  these  reactions  is  potentially  feasible. 

1.6.3  Pulsed  HF  Chemical  Laser  Initiated  by  an  Electric  Discharge 

1 .  6.  3.  1  Introduction 

Although  the  electron-pii  discharge  technique  could  not  be  used 
to  volumetrically  initiate  chemical  reactions,  it  is  relatively  simple  and  easy 
to  construct.  The  object  of  this  experiment  was  to  evaluate  the  efficiency  of 
electrical  initiation  of  H2/F2  mixtures  to  produce  laser  action.  It  was  hoped 
that  enough  information  could  be  obtained  so  that  the  overall  electrical  effi¬ 
ciency  can  be  derived. 

1.6.  3.2  Experimental 

The  gas  mixing  technique  and  optical  cavity  used  is  identical  with 
that  used  in  the  fla  sh  -  lamp  -  initiated  H2/F  laser  experiment.  The  two  Xe 
flash  lamps  were  replaced  by  two  pin-discha rg  i  teflon  cavities.  Each  cavity 
contains  90  1  kV  resistors  located  within  a  15  mm  x  250  mm  rectangular 
cross  section.  An  8  mm  diameter  brass  bar  was  used  as  an  anode.  The 
high  voltage  power  supply  and  pulsing  unit  orijinally  used  for  Avco  N2  laser 
(which  puts  out  100  nsec  pulses)  were  used  to  drive  and  trigger  the  electric 
discharge.  Single  0.  0l2jjf  capacitors  or  a  pair  of  0.  012/Lif  capacitors  coupled 
in  series  were  used  for  varying  the  discharge  energy.  The  laser  output 
energies  were  measured  with  an  Avco  energy  meter,  and  the  pulse  shape 
was  monitored  by  an  InSb  PEM  detector. 

1 .  6 .  3 .  3  Re  suits 

Measurement  of  both  laser  output  energy  and  pulse  shape  on 
various  H2/F2/He  mixtures  were  performed  under  two  experimental  condi¬ 
tions: 

(1)  The  electric  input  energies  were  varied  by  changing 
the  discharge  voltage  across  the  capacitor. 

(2)  The  electric  input  energies  were  varied  by  changing  the 
capacitance  at  a  fixed  charging  voltage. 

Figure  43  shows  a  plot  of  the  laser  output  energy  versus  capacitor 
charging  voltage  observed  in  a  0.  08/0.  08/0.  84  He  mixture.  It  is  clearly 
demonstrated  that  the  relation  between  the  electric  input  energy  and  laser 
output  energy  is  not  linear.  By  decreasing  the  capacitance  from  0.012/i 
to  0.  006^f,  at  a  charging  voltage  of  10  kV,  the  laser  output  energy  was  found 
to  drop  only  slightly  as  indicated  in  the  figure.  Figure  44  shows  a  typical 
HF  laser  pulse  observed  in  the  mixture  of  H2/F2/He  =  0.  04/0.  04  0.  92.  The 
charging  voltage  was  17.  5  kV  and  capacitance  was  0.  012/Ltf.  An  output  energy 
of  0.  13  joule  was  observed. 
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ENERGY  (JOULE) 


CHARGING  VOLTAGE 


Fig.  43  The  Measured  Laser  Output  Energy  Versus  Capacitor  Voltage 
Observed  in  the  H^/F^/He  =  0.  08/0.  08/0.  94  Mixture 
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!•  6.  3.  4  Onthe  Overall  Electrical  Efficienc 


In  order  to  evaluate  the  total  elect  ronic -to- la  se  r  energy  effi¬ 
ciency,  the  exact  amount  of  electric  energy  that  was  deposited  into  the  gas 
mixture  has  to  be  carefully  calculated.  Let  -  total  resistance,  R„ 
resistance  across  the  gas  mixture,  Rr  =  resistance  across  resistor,  t 
capacitance,  V  charging  voltage,  EjnpUt  total  electric  input  energy,  Er 
total  energy  deposited  in  resistors,  E„as  =  total  energy  deposited  in  pas  and 
I  (t)  transient  current.  Then 


Rr  780  =  ohm  (81) 


R.  =•  R  +  R  (t) 
t  r  g' 


E 

r 


[  I(t)]  2  Rr  dt 


(82) 


(83) 


where  transient  current  is  assumed  to  be  I(t)  =  I  exp  (-t/RtC). 

or 

exp  (-2t/RtC)  dt  (84) 

o 

Since  Rt  is  also  a  function  of  t  and  it  is  not  clear  what  is  the  furmulation  of 
Rg,  the  exact  integration  is  difficult.  However,  if  Rr  »  R  on_*  can  approxi 
mate  Rt  w  Rr  then  K 


E 

r 


R 

r 


RtC 


E  _ 1 

input  l  +  ( R  /R  ) 
R  r' 


If  Rg  «  Rr, 


Er  ~  Einput- 


(85) 


(86) 

(87) 
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Unfortunately,  due  to  the  shortage  of  time,  the  exact  value  of 
Rt  or  Rp  was  not  accurately  measured.  If  one  assumes  that  Rg  «  Rr,  the 
lower  limit  on  the  overall  electrical  efficiency  for  a  5  kV  discharge  is 


Laser  output  ene rgy  _  56  mjoule 

Electrical  input  energy  150  mjoule 


40% 


where  most  of  the  electrical  energy  was  deposited  in  the  resistors  and  not 
the  gas. 


If  one  assumes  that  R^  ~  Rr,  the  ratio  of 

Laser  output  energy  Laser  output  energy  _  „  , 

Energy  deposited  in  gas  EL  ~7Z 

nj  r  n  input' 

Another  emperical  factor  that  is  important  but  cannot  be  accur¬ 
ately  measured  is  the  active  discharging  volume  of  gas.  If  one  assumes  that 
the  discharge  does  not  disperse  too  far  in  a  2  cm  distance,  an  active  excita 
tion  volume  of  3  cm  ^  is  obtained.  Thus  the  output  energy  is  0.  14  joule,  3  cm^ 
ss  50  joules/f  in  the  mixtures  of  H2/F2/He  =  8%/8%/84%.  The  technique  of 
electric  discharge  does  look  more  efficient  than  flash  lamp  initiation  of  HF 
chemical  lasers. 

1.6.4  Monoe  "rgetic  Electron  Beam  Dissociation  of  Fluorine 

1.6.4.  1  Int  roduction 

Initiation  of  chemical  laser  action  by  c  electron  beam  has  been 
reported  by  Gregg,  et  al.  (  22)  Since  no  values  of  lase  r  pulse  energy  were 
quoted  by  these  workers,  it  is  impossible  to  estimate  how  efficient  the  elec¬ 
tron  beam  is  as  an  initiation  source,  other  than  qualitatively  it  appeared 
better  than  flash  lamp  initiation.  Gregg,  et  al.  ,  used  fluorine  containing 
compounds  (I^F^,  NF3)  as  their  source  of  fluorine  atoms,  and  it  is  not 
known  whether  a  chain  or  single  step  reaction  resulted  from  the  initiation. 
Mixtures  of  fluorine  and  hydrogen  do  combine  to  form  HF  in  a  chain  reaction 
and,  therefore,  require  less  initiation  energy  to  go  to  completion  than  a 
simple,  one  step  reaction.  Laser  action  from  F mixtures  has  been 
previously  demonstrated.  (23) 

In  this  study,  a  high  energy  electron  beam  was  injected  into  a 
gaseous  mixture  of  F^,  He,  CO;>,  and  HCl.  Particular  attention  was  focused 
on  the  determination  of  the  number  of  fluorine  atoms  produced  by  means  of 
a  chemical  titration.  The  measurement  of  the  fluorine  atom  concentration 
was  preferred  to  a  measurement  of  laser  pulse  energy,  simply  due  to  the 
fact  that  the  latter  depends  on  cavity  coupling  and  mode  volume  as  well  as 
the  reaction  rate  of  the  atoms  produced,  and,  therefore,  is  less  direct. 
Previous  detailed  studies  on  the  efficiency  and  mechanism  of  HF  and  HCl 
chemical  lasers(24,  25)  have  been  frequently  hampered  by  the  uncertainty  of 
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the  atom  concentration  involved  in  the  initiation  step.  The  quantitative  mea¬ 
surement  of  fluorine  atom  concentration  produced  by  high  energy  electrons 
reported  here  should  be  important  to  the  further  development  of  laser  science 

and  technology. 

1 ,  6.4.  2.  Expo  rimental 

The  goal  of  this  experiment  was  to  use  a  high  energy  electron 
beam  to  dissociate  F2  and  then  to  quantitatively  measure  the  F-atoms  pro¬ 
duced.  From  a  knowledge  of  the  energy  deposition  of  the  electron  bean  in 
the  gas,  the  energy  necessary  to  produce  an  F -atom  can  then  be  determined. 
Since  the  amount  of  F2  dissociation  produced  was  expected  to  be  small, 
e.  g.  of  the  order  of  a  percent  or  so,  a  sensitive  and  direct  measurement 
of  the  F-atom  concentration  was  desirtd.  It  was  decided  to  use  a  chemical 
titration  technique  based  on  the  known  •ast  reaction 


F  +  HC1  —  HF  +  Cl  (88) 

The  product  HF  produced  in  this  reaction  could  be  quantitatively  measured 
in  absorption  via  a  probe  HF  laser.  However,  since  F2  and  HC1  react 
spontaneously,  but  slowly,  a  flow  system  to  mix  the  reagents  and  flow  theni2 
rapidly  into  the  reaction  chamber  must  be  employed.  Such  a  flow  system , 
constructed  for  experiments  involving  mixtures  of  H2  and  F2  to  investigate 
laser  action,  was  already  available  at  AERL.  It  was  decided  to  employ  this 
mixing  system  and  couple  it  to  an  existing  electron  beam  device  available 
from  Simulation  Physics,  Inc.  (SPI).  This  particular  electron  gun  had 
several  advantages:  (1)  it  was  available  and  easily  transportable  to  the  ex¬ 
periment  since  it  employs  a  dry  dielectric  energy  storage  system;  (2)  the 
size  and  dimensions  of  the  emitting  surface  were  reasonably  compatible 
with  the  flow  system  geometry;  (3)  the  gun  was  of  the  plasma  diode  type, 
producing  a  high  current  density  to  maximize  the  dissociation  of  F-atoms; 
and  (4)  SPI  had  suitable  calibration  techniques  and  calculation  procedures 
available  to  determine  the  energy  deposition  of  the  electron  beam  in  the  gas. 
Although  this  combination  of  existing  flow  system  and  electron  gun  seemed 
the  most  expeditious  and  feasible  approach  and  did  result  in  experimental 
data,  it  will  be  obvious  that  a  more  optimum  mating  of  E-beam  and  gas 
mixture  could  have  been  achieved. 

A  schematic  of  the  experimental  apparatus  is  shown  in  Fig.  45. 

The  fluorine  and  hydrogen  chloride  are  premixed  with  He  diluent  in  separate 
storage  tanks.  Both  storage  tanks  are  pressurized  well  above  one  atmosphere. 
The  remainder  of  the  flow  system  is  pressurized  to  about  l  atm  with  the 
diluent.  During  a  run,  fast  acting  valves  open  to  admit  the  gases  to  a  mixing 
manifold  which  distributes  the  gases  in  alternate  sheets  of  F2  and  HC1  about 
2  mm  thick.  The  gases  mix  as  they  flow  through  the  mixing  region,  the  re¬ 
sidence  time  being  typically  20  milliseconds.  The  mixed  gases  then  pass  into 
the  interaction  chamber,  where  the  electron  beam  is  fired  transverse  to  the 
direction  of  the  gas  flow.  The  spent  gas  exhausts  through  traps  for  removing 
the  HF.  The  probe  HF  laser  enters  the  interaction  chamber  transverse  to 
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both  thi'  direction  of  the  electron  beam  and  the  gas  flow.  A  particular 
vibration  rotation  transition  of  HF  is  selected  by  a  monochromator  and  mea  - 
sured  with  an  infrared  detector. 

The  electrons  were  produced  by  a  SPI  electron  beam  generator, 
capable  of  producing  electrons  with  a  beam  energy  ranging  from  20  to  140 
keV,  and  a  current  density  at  the  anode  from  100  to  600  amperes/cm^  Pulse 
widths  of  the  electron  beam  were  typically  250  nanoseconds.  Figure  46  shows 
an  actual  photograph  of  the  SPI- PULSE  2500  electron  beam  generator  used 
in  this  experiment.  Figure  47  shows  the  details  of  the  working  section  of  the 
gun  and  its  mating  to  the  interaction  chamber.  The  electron  gun  consists  of 
a  DC  charged,  low  impedance,  coaxial  structure  which  is  switched  into  a 
field  plasma  emission  tube  via  a  triggered  gas  gap.  The  1  ohm,  15,000  pf 
coaxial  energy  storage  has  a  characteristic  transit  time  of  30  nanoseconds. 
Charged  to  250  kilovolts  by  an  electrostatic  generator,  the  system  stores 
approximately  450  joules.  After  about  one  minute  charging,  a  command 
pulse  initiates  the  discharge  of  the  line  into  the  field  emission  tube.  Electrons 
emitted  from  the  cathode  fall  through  the  instantaneous  cathode -anode  poten¬ 
tial  difference,  strike  and  pass  through  the  transmission  anode,  composed  of 
a  0.0005"  thick  aluminum  foil  and  tungsten  wire  supporting  screen,  and  into 
the  experimental  chamber.  The  diode  operates  at  a  pressure  of  10'^  torr 
while  the  gas  channel  contains  the  gas  mixture  at  760  torr.  After  passing 
through  the  aluminum  foil  window,  the  electron  beam  drifted  through  a  4  cm 
long  X  l  cm  wide  "neck"  into  the  ‘  cm  wide  interaction  chamber.  The  dimen¬ 
sions  of  the  emitter,  the  transmission  anode,  and  "neck"  were  all  1  X  10  cm. 

To  diagnose  the  electron  beam  character  during  each  discharge  of 
the  generator,  both  diode  current  and  voltage  were  monitored.  The  diode 
current  monitor  is  a  low  inductance  resistive  belt  which  is  placed  circum¬ 
ferentially  at  the  cathode -anode  gap  in  the  wall  of  the  diode  (see  Fig.  47). 

The  voltage  monitor  is  a  capacitive  divider  formed  by  a  metal  band  and  the 
wall  of  the  diode  (see  Fig.  47).  The  voltage  monitored  is  that  of  the  cathode 
shank.  Inductive  corrections  to  this  signal  are  also  made.  Typical  diode 
current  and  voltage  data  are  shown  in  Fig.  51. 

In  addition  to  diode  current  and  voltage,  separate  experiments 
were  performed  to  assess  analytical  techniques  used  to  determine  the  degree 
of  ionization  created  in  the  neck.  These  experiments  consisted  of  measuring 
current  transmitted  through  the  neck  and  the  energy  per  gram  or  dose  created 
by  the  electrons  in  a  dosimeter  also  at  the  back  of  the  neck.  The  transmitted 
current  was  measured  using  a  Faraday  cup.  The  dose  measurement  was 
made  with  blue  cellophane  dosimetry.  A  typical  transmitted  current  mea¬ 
surement  is  also  shown  in  Fig.  51.  Figure  52  is  a  blue  cellophane  map  of 
the  current  density  distribution  as  determined  from  a  calibration  against 
dose  at  the  back  of  the  neck,  which  clearly  indicates  the  uniformity  of  the 
electron  beam  that  the  SPI-PULSE  machine  produced. 

For  monitoring  the  quantity  of  hydrogen  fluoride  produced  through 
the  rapid  reaction  (88),  the  absorption  of  the  1  P^  transition  of  a  HF  pulsed 
laser  was  measured.  The  orobe  HF  laser  was  an  Avco  C950  nitrogen  laser, 
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Experimental  Geometry  for  Gas  Excitation  Experiment. 
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Fig.  48  Typical  1 P4  Laser  Pulses,  Monitored  by  a  GeAu  Photoconductor. 

The  upper  trace  is  the  intensity  of  IP4  line  before  the  HF  ab¬ 
sorption.  The  lower  trace  is  the  intensity  of  that  after  the  IR 
absorption  of  HF.  The  optical  pathlength  is  10  cm.  A  ratio  of 
I/I  =  0.  03  is  observed  which  represents  a  partial  pressure  of 
0.  Ob  torr  of  HF  in  the  mixture  studied.  The  charging  voltage 
was  100  kV.  The  gaseous  composition  was  4%  F  1%  HC1, 

1%  CO2  and  94%  He  at  a  total  pressure  of  800  torr. 
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modified  to  operate  on  SF^  and  H2  to  produce  HF  laser  pulses.  Figure  48 
shows  typical  1  P4  laser  pulses ,  monitored  by  a  GeAu  detector  through  a 
0.  5-meter  Jarrel  Ash  monochromator,  showing  the  absorption  of  the  HF 
molecules  in  the  interaction  chamber  before  and  after  firing  the  E-beam, 
hue  to  the  possibl  mismatch  between  the  line  center  of  IR  absorption  and 
’aser  emission,  the  absorption  coefficient  of  the  1  P4  line  of  HF  was  cali¬ 
brated  against  the  actual  absorption  of  a  known  HF/He  mixture  in  a  separate 
experiment.  Figure  51  shows  the  measured  F/F2  ratio  as  determined  from 
the  1  P4  absorption  for  various  electron  beam  discharge  voltages.  The  scatter 
in  F/F2  ratio  observed  is  caused  mainly  by  fluctuations  in  power  of  the  1  P4 
probe  laser  transition  from  shot  to  shot. 

A  diagram  illustrating  the  overall  sequence  of  timing  and  trigger¬ 
ing  is  shown  in  Fig.  50.  Initially  the  fast  acting  valves  are  opened  and  the 
gases  begin  to  flow.  Approximately  100  msec  later  the  mixed  gases  are  en¬ 
tering  the  interaction  chamber.  The  electron  gun  is  fired  at  about  200  msec 
after  the  valves  opened;  about  150/isec  later  the  probe  HF  laser  is  fired. 

This  latter  period  is  needed  to  assure  that  the  HF,  produced  vibrationally 
excited  by  reaction  (88),  has  relaxed.  This  point  is  discussed  further  in  the 
next  section.  During  this  150  msec  the  gas  in  the  interaction  chamber  has 
moved  only  1-?  mm,  so  that  for  all  purposes  the  gas  sampled  by  the  E-beam 
and  laser  probe  are  identical. 

1.  6.  4.  3  Dissociation  and  Analysis  Scheme 

A  primary  high  energy  electron  incident  upon  a  gas  molecule 
may  lose  energy  in  three  ways:  ionization,  excitation,  and  dissociation  of 
the  target  molecule.  The  debris  produced  by  a  high  energy  electron  beam 
then  consists  of  a  swarm  of  electrons  and  ions,  together  with  excited  mole¬ 
cules  and  molecular  fragments.  Since  our  gas  is  mainl/  helium,  presumably 
most  of  the  impacts  will  lead  to  the  ionization  and/or  excitation  of  helium 
atoms  ,  i.  c .  , 


e  + 
P 


He 


He 


+  e  + 
P 


(»9) 


He 


(90) 


where  Op  and  eg  are  the  primary  and  secondary  electrons,  respectively, 
and  He"  stands  for  electronically  excited  species,  which  could  be  any  of  the 
(2  ^S,  2  ^S,  2  lp,  2  ^P)  states.  Landshoff  and  Magee(27)  have  shown  that, 
in  the  case  of  air,  the  primary  electron  mainly  undergoes  ionizing  collisions, 
rather  than  dissociative  collisions,  resulting  in  secondary  electrons  with 
approximately  22  eV  energy.  A  90  keV  primary  electron  would  be  capable 
of  producing  4  X  10  5  secondary  electrons,  which  clearly  indicates  that  the 
number  of  secondary  electrons  must  be  10^  times  that  of  the  primary  elec¬ 
trons.  Thus  for  the  F2/H0  mixtures  studied,  dissociation  of  F2  will  probably 
be  caused  mainly  by  the  secondary  electrons.  The  processes  which  could 
lead  to  the  ionization  and/or  dissociation  of  the  F2  molecule  during  the 
electron-molecule  collisions  are 
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Fig.  50  Schematics  of  Overall  Experimental  Timing  and  Triggering. 
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Representative  Diode  Current,  Voltage  and  Transmitted  Current 
Data.  The  charging  voltage  was  120  kV  and  cathode  area  is 
1  c  m  x  10  cm. 
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Fig.  52  Current  Density  Distribution  Behind  Gas  Channel.  The  charging 
voltage  is  120  kV,  and  cathode  area  is  1  cm  x  10  cm. 
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F2  i  e  -  f2‘  -  F  +  F" 

-  F2+  +  2e 
—  F+  +  F  +  2o 

—  F+  +  F"  +  e 

-  2F  +  c 


Ao 

-1.9  eV 

(91) 

Ao 

16.  6  eV 

(92) 

> 

o 

ii 

19.0  eV 

(9  3) 

> 

o 

ii 

15.  4  eV 

(94) 

Ao  = 

1.  5  eV 

(95) 

where  An  is  the  minimum  energy  required  for  the  forward  reaction,  i.  e.  , 
the  minimum  energy  that  would  have  to  be  supplied  by  the  colliding  electrons. 
For  processes  (92),  (93)  and  (94),  impacts  with  electrons  of  moderate  energy 
are  required.  For  processes  (89)  and  (95),  collisions  with  electrons  of  only 
low  energy  are  necessary.  Since  process  (91)  is  1.9  cV  exothermic  and  re¬ 
quires  small  (or  no)  activation  energy  it  is  expected  to  be  relatively  fast 
Using  a  shock  tube  technique,  Mandl<28)  has  measured  the  associative  detach¬ 
ment  rate  for  the  reaction  F  +  F"  -  F2  +  e  -  1 . 9  eV ,  i.  e.  the  reverse  of  (91). 
A  rate  constant  of  10-12  cm3  scc-l  at  4000PK  was  observed.  The  reaction 
rate  estimated  from  Mandl' s  data  for  reaction  (91)  at  room  temperature  indi¬ 
cates  that  reaction  (91)  must  be  fast  and  has  a  cross  section  of  the  order  o 
10-15  cm2.  Fox(29)  and  coworkers  have  actually  measured  the  cross  section 
for  dissociative  electron  attachment  of  the  I2  molecule.  A  value  of  3  *  10 
cm2  was  observed,  comparable  to  our  estimate  for  reaction  (91). 

Processes  like  charge  transfer  or  Penning  ionization  from  He 
could  also  be  of  importance  for  dissociation,  since  for  '"'cry  secondary 
electron  liberated  from  helium  there  is  also  an  accompanying  ion. 


He+  + 

F^  -  F/  +  He 

A q  =  -8.  0  eV 

(96) 

'i' 

He  + 

F2  —  F2+  +  He  +  e 

Ac  =  -3.  1  eV 

(97) 

He"  + 

F2  -  F2 -*  2F  +  He 

A q  =  -18.  2  eV 

(98) 

sta  nds 

for  metastable  Helium 

atom.  The  reaction  of  ionic 

re  - 

Here,  He 
combination  or  neutralization, 


+  F'  —  F  +  F+F  (99) 

could  generate  additional  fluorine  atoms.  It  is  rather  unfortunate  that  none 
of  fhese  reactions  have  been  studied.  At  the  present  stage,  we  can  only  sus¬ 
pect  that  reactions  (91),  (92)  and  (99)  are  most  likely  the  important  processes 
involved  in  making  F-atoms  via  high  energy  electron  bombardment  of  F  2- 

The  reaction  between  an  F-atcm  and  an  HC1  molecule  has  been 
thoroughly  studied  by  many  investigators.  A  recent  report  by  Kompa  and 
Wanneri 30)  has  indicated  that  the  reaction 
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F  +  HC1  — ■  HF  +  Cl 


(88) 


is  fast  and  possesses  a  rate  constant  of  1.5  X  101  3  cm  5  mole"  1  sec"  1  at 
300°K.  Thus  for  gaseous  mixtures  that  we  have  studied,  the  titration  reac¬ 
tion  should  take  only  a  few  hundred  nsec  to  reach  completion.  Furthermore, 
the  recombination  of  F-atoms  is  too  slow  and  unimportant  under  our  condi¬ 
tions  . 


For  a  quantitative  and  accurate  measurement  of  the  concentration 
of  F-atoms  originating  from  the  electron  beam  dissociation,  any  side  reac¬ 
tions  that  might  lead  to  the  formation  of  HF  must  also  be  considered  or 
corrected  for.  Since  there  are  by  necessity  a  significant  number  of  HC1 
molecules  present  in  the  mixtures  studied,  ionization  or  dissociation  of  HC1 
would  lead  to  further  complication  of  our  titration  scheme.  The  cross  sec¬ 
tion  for  dissociative  electron-attachment  for  HC1  has  been  carefully  studied 
by  Christophorou,  ( 3  *)  et  ai.  It  is  found  that  the  overall  cross  section  for 
the  reaction 

HC1  +  e  — -HC1"  -*■  H  +  Cl",  AQ=0.81eV  (100) 

is  7.4  X  lO"1^  cm2  eV  and  peaks  a t  a n  impact  energy  of  0.81  e V  with  a  maxi¬ 
mum  cross  section  of  1.  85  X  10"17  cm2.  If  this  cross  section  is  much  larger 
than  that  of  reaction  (91)  for  F2,  reaction 

H  +  F2  —  HF  +  F  (101) 

plus  reaction  (88)  would  lead  to  the  formation  of  two  additional  HF  molecules 
for  every  dissociative  electron-attachment  of  HC1.  Although  reaction  (  101) 
is  about  seven  times  slower  than  that  of  reaction  (88),  the  measurements 
were  made  on  samples  enriched  with  molecular  F2  ([  F2]  /|  HC1]  a?  4),  there¬ 
fore  it  is  not  all  that  clear  whether  contributions  from  reaction  (100)  to  the 
overall  formation  of  HF  molecules  are  negligible.  Other  side  reactions 
like  charge  transfer,  Penning  ionization  of  HC1,  and  ionic  recombinations, 

Hc+  4  HC1  —  HC1+  +  He  A0=-ll.8eV  (102) 

He  +  HCl  —  HC1+  +  He  +  3  AQ  =  -6.  9  eV  (103) 

HC1+  +  F"  —  HF  +  Cl  or  H  +  F  +  Cl  (104) 

would  further  hamper  the  interpretation  of  the  e.vperimenta  1  data.  Since 
reaction  (91)  is  likely  to  be  a  hundred  times  faster  than  any  of  the  above  side 
reactions  mentioned,  we  conclude  that  the  analysis  should  be  simple  and 
straightforward. 

Since  HF  molecules  generated  from  either  reaction  (88)  or  (101) 
are  vibrationally  excited,  CO2  gas  was  added  to  the  mixture  to  shorten  the 
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vibrational  relaxation  time  of  HF.  Then,  the  1 absorption  measurement 
could  be  performed  on  ground  state  HF.  Hancock  and  Grcen(^)  have  re 
cently  measured  the  rate  of  vibrational  relaxation  of  HF  by  CO^.  They  find 
that  the  reaction 

HF  (v  =  1)  +  COz  —  HF  ( v  =  0)  +  CO^  (n,  1°,  m)  (105) 

has  a  c  ross  section  of  0,27  x  10-1^  cm2  or  a  rate  constant  of  5.9  X  104  sec" 1 
torr-^.  For  mixtures  with  8  torr  of  CO^,  HF  vibrational  relaxation  should 
take  only  a  few  psec.  Since  we  have  delayed  our  absorption  experiment  for 
about  150  /nsec  after  the  E-beam  pulse  (see  Fig.  52),  the  amount  of  IP4  absorp¬ 
tion  should  give  us  a  direct  measurement  of  the  concentration  of  ground  state 
HF  molecule  produced. 

1.6.4.  4  Results  and  Discussion 


Table  VIII  lists  the  composition  of  the  gaseous  mixtures  studied 
and  the  amount  of  F-atoms  (or  HF  molecules)  observed  for  various  F^/HCl/ 
CG2/He  mixtures  for  five  different  electron-impact  energies.  For  a  fixed 
gas  composition  of  F2/HCl/C02/He  =  l/T  25/0.  25/23.  5,  the  amount  of  F- 
atoms  produced  varies  linearly  with  the  amount  of  electron  energy  deposited 
in  the  gas  as  indicated  in  Fig.  51.  In  order  to  obtain  a  gross  figure  on  the 
efficiency  of  F2  dissociation  by  electron  impact,  the  number  of  ionizing  colli¬ 
sions  produced  in  the  gas  by  the  electron  beam  has  to  be  estimated. 

For  a  monoenergetic  electron  beam,  if  the  stopping  power  and 
the  average  energy  loss  per  ionizing  collision  for  the  gas  sample  are  known, 
one  can  correlate  the  number  of  ionizing  collisions  per  cm^,  Ncon  by 


where  j(t)  is  the  total  current  density  (in  electrons/cm2  sec)  as  a  function 
of  time,  t.  T  is  the  pulse  duration  of  the  electron  beam.  W  is  defined  as 
the  average  energy  loss  per  ionizing  collision  that  the  primary  electrons 
suffer.  The  energy  degradation  per  unit  path  (dE/dx)  can  be  expressed  as 

v'M  ,107) 

i 

where  S-  (E)  is  the  stopping  power  for  molecule  i  at  an  electron-impact- 
energy  of  F.  p-  is  the  density  for  the  gas  molecule  i  and  the  summation  is 
over  all  gaseous  species. 

To  determine  the  electron  beam  energy  spectra  and  electron 
energy  deposition  profile  for  these  experiments,  a  separate  calibration  ex¬ 
periment  and  detailed  computation  were  performed.  For  simplicity,  the 
approach  used  is  summarized  as  follows: 
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ELECTRON  BEAM  DISSOCIATION  OF  FLUORINE  IN  F2.  MCI.  C02  AND 
He  MIXTURES  AT  A  TOTAL  PRESSURE  OF  «00  TORR 
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The  average  values  with  a  possible  error  up  to  _+  25%. 


(1)  Measure  diode  current  and  voltage. 

(2)  Calculate  the  electron  beam  spectrum  using  the 
EBSPEC  code.  (33) 

(3)  Calculate  the  energy  deposited  in  the  gas  mixture 
and  the  current  density  transmitted  through  the 
mixture  by  the  electron  beam  using  the  Monte 
Carlo  transport  code  ELTRAN.  (33) 

( -P  Measure  the  transmitted  charge  and  compare  it 
to  that  calculated  to  determine  the  fraction  of 
current  lost  to  the  walls  by  scattering  in  passing 
through  the  neck  and  into  the  channel. 

(5)  Use  this  fraction  to  correct  the  calculated  energy 
deposited  for  the  wall  losses. 

(6)  Measure  the  energy  deposited  at  the  back  of  the 
neck  using  blue  cellophane  dosimetry  and  compare 
to  the  corrected  calculation  to  verify  the  analysis. 

Figures  53a,  53b  and  53c  show  a  typical  experimental  diode 
voltage  (corrected  for  inductance)  and  diode  current  and  the  calculated 
electron  beam  spectrum  determined  by  EBSPEC  for  the  generator  charging 
voltages  of  100,  120  and  140  kV  employed  in  this  experiment. 


These  spectra  were  then  used  to  calculate  both  (1)  the  electron 
beam  energy  deposition  profiles  through  the  experimental  package  by  using 
ELTRAN,  and  (2)  the  number  of  ionizing  collisions  per  cm3,  Ncon  by  Eq. 
(106).  As  is  clearly  shown  in  Figs.  53a,  53b  and  53c,  the  electron  beam 
energy  injected  into  the  gas  sample  is  not  monoene rgetic.  The  distribution 
function  f(E)  is  thus  inserted  into  Eq.  (106)  which  becomes 
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f  (E)  j  (t) 
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pi 

S .(E)  ^  dt  dE 


(108) 


Figures  54a,  54b  and  54c  show  the  calculated  electron  energy  deposition  pro¬ 
files  for  the  charging  voltages  of  140,  120  and  100  kV,  respectively.  All 
these  profiles  have  been  normalized  to  unit  energy  fluence  at  the  anode. 

Table  VIII  lists  the  values  of  Ncqjj  obtained  through  numerical 
integration  of  Eq.  (108)  by  using  ELTRAN  and  by  adopting  the  corrected 
calculated  energy  deposition  as  shown  in  Table  IX.  The  values  of  Sj_je  and 
S (j (p,  were  taken  from  Berger  and  Seltzer,  (  3^)  and  the  values  of  Wj_je  and 
from  Bartner  and  Hurst.  (^)  The  value  of  Wj_j£q  was  taken  from 
Chr istopho rou.  01)  For  he lium,  due  to  pos sible  excitation  effects ,  0  1  )  the 
value  of  W^j  used  is  the  value  appropriate  for  helium  with  added  gas,  not 
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Typical  Electron  Beam  Spectrum  for  a  Charging  Voltage 
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Fig.  54c  Typical  Electron  Beam  Energy  Deposition  Profile  for  the 
Experiment  with  a  Charging  Voltage  of  100  kV. 


T ABLE  IX 


DATA  SUMMARY  FOR  ELECTRON  REAM  EXCITATION  EXPERIMENT 


Pa  ramete  r 

100  kV 

Char^inR  Voltape 

120  kV 

Avcragi*  Electron  Energy  in  Diode 
from  EBSPEC 
koV 

39 

43 

Fluence  Before  Sc  reen 
(  al/tm^ 

0  93 

1.11 

Flue  nee  Alter  Sc  reen  (75%) 
t  al/ctn^ 

0.  63 

0.  75 

Avvraiiu  Dose  to  "Mis  -  Calculated 
by  ELTRAN 

'  al/pr 

45,8 

64.  1 

ETOfAL  ln  Dlod*' 
joule  s 

38.  7 

43.  3 

qtotal  in  n,otU' 

1000 

1075 

^TOTAl  A^ter  Sc  reen  (75?o) 

7  50 

806 

Q  Transmitted  from  Channel- 
Calculated  by  ELTRAN 

I  32 

187 

Q  Transmitted  from  Channel  - 
Measured 

5.  1 

13.  5 

Q  Transmitted  Measured 

Q  Transmitted  Calculated 

0.  39 

072 

Dose  ”Mi\"  = 

°TM  .  „  ,  , 

pr -  ^  Dose /cal 

UTC 

t  al/«r 

1  .  8 

4. 6 

Average  Dose  to  "Mix"  -  Measured 
c  al  /iir 

1 .  9 

7.  3 

N 


i  oil 


14 


14 


140  kV 

S3 

1 . 72 

l  ib 

US 

72.0 

1  3S0 

1012 

383 

21.9 

057 

7,  1 

8  8 


clot  Iron  collisions/t  in ^ 


3X10 


7.  6  x  10 


11  8X10 


the  value  for  pure  helium.  The  value  of  Wp  .  is  not  known,  but  the  fluorides 
BFj  and  SF ^  have  a  value  of  36  eV,  as  do  omer  diatomic  molecules,  so  this 
value  was  assumed  for  fluorine.  Both  Sp  and  S |_j values  are  not  known 
either,  the  values  corresponding  to  Ne  ana  Ar  were  Chosen,  respectively. 

Table  IX  summarizes  the  results,  most  of  which  are  self- 
explanatory.  As  discussed,  the  approach  was  to  get  an  empirical  correction 
factor  for  the  calculated  energy  deposited  in  the  mixture  by  ratioing  the  mea¬ 
sured  to  the  calculated  transmitted  charge.  Applying  this  correction  to  the 
calculated  energy  deposited  brought  the  value  close  to  that  measured  as  is 
shown  in  Table  IX.  This  agreement  gave  confidence  to  the  analytical  approach 
taken.  Table  VIII  lists  the  values  of  Ncqjj  calculated  and  the  corresponding 
number  of  F-atoms  observed.  A  consistent  ratio  of  [  F]  /N  „  «  3  is  indicated. 
Thus,  the  comparison  between  the  experiment  and  the  calculation  suggests 
that  about  three  fluorine  atoms  are  produced  for  each  electron-molecule 
collision.  If  one  assumes  that  high  energy  electrons  lose  35  eV  per  inelastic 
collision,  it  then  costs  approximately  12  eV  for  every  F-atom  produced. 

Since  the  dissociation  scheme  is  obviously  quite  complicated, 
at  the  present  time  with  our  limited  information,  it  is  not  possible  to  draw 
any  conclusions  regarding  the  mechanism  of  dissociation  or  the  percentage 
of  contribution  from  each  individual  reaction. 

1 .  6.4.  5  Conclusion 

The  experiments  described  have  determined  the  efficiency  for 
electron  beam  dissociation  of  F2  molecules  in  HC1,  CO^  and  He  gaseous 
mixtures.  An  efficiency  of  three  F-atoms  per  electron-molecule  collision 
is  observed.  If  one  assumes  that  high  energy  electrons  lose  35  eV  per 
inelastic  collision,  it  costs  approximately  12  eV  for  every  F-atom  produced. 
Thus,  an  electron  beam  does  apnea  r  to  be  an  efficient  means  for  initiation 
of  HF  chemical  lasers.  A  few  possible  processes  that  might  lead  to  F^ 
dissociation  we  re  suggested.  However,  due  to  the  lack  of  knowledge  on 
rate  constants,  it  i3  not  possible  to  give  any  firm  or  quantitative  conclusion 
on  the  mechanism  or  the  relative  importance  of  each  individual  process  in¬ 
volved  in  the  reaction.  Nevertheless,  the  information  obtained  from  the 
experiment  should  be  useful  to  the  further  development  of  electron  beam 
dissociated  chemical  lasers,  and  hopefully  will  stimulate  additional  detailed 
studies  of  high  energy  electron  kinetics  in  various  gases. 
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1.7  CONCLUDING  REMARKS 


At  the  beginning  of  this  Contract  the  future  of  the  high  energy, 
high  pressure,  pulsed  HF  laser  did  not  appear  bright.  Although  mixture 
stability  had  been  achieved,  there  were  problems  of  reproducibility. 

Although  laser  action  had  been  demonstrated,  the  efficiencies  were  very 
poor,  and  more  discouraging  was  the  apparent  lack  of  understanding  of  the 
data  obtained.  Certainly  the  progress  achieved  during  the  past  year  and 
reported  in  this  Final  Technical  Report  is  impressive  and  definitely  suggests 
the  potential  of  tnis  type  of  device  for  lutvre  military  applications. 

The  lase  "  performance  demonstrated  on  a  small  scale  -  8  joules 
in  a  100  cc  flash  initiated  device  for  a  mixture  of  10%  F2/  10%  H2/80%  He 
at  a  total  cavity  pressure  of  1.  1  atm  -  is  a  milestone  in  pulsed  chemical 
laser  technology.  This  performance  scales  to  80  j/  1-atm  and  a  specific 
energy  of  100  kj/lb  which  is  certainly  as  good  or  letter  than  demonstrated 
by  any  other  laser;  electric,  gas  dynamic,  or  othe  r  chemical  system. 

Another  important  milestone  that  has  been  achieved  is  the  develop¬ 
ment  of  a  kinetic  model  that  has  provided  considerable  insight  into  the  detailed 
chemistry  and  physics  and  performance  of  this  laser  system.  There  are 
still  a  number  of  important  areas  of  disagreement  between  model  prediction 
and  observed  laser  characteristics,  but,  considering  the  complexity  of  this 
system,  the  significance  of  the  agreements  far  over  shadows  the  areas  of 
uncertainty.  There  are  obvious  areas  of  future  research  that  should  be 
pursued  in  theoretical  model  development. 

The  third  significant  achievement  under  this  program  was  the 
first  experimental  measurement  of  the  efficiency  of  producing  fluorine  atoms 
from  F2  via  a  high  energy  electron  beam.  A  value  of  12  eV/F-atom  was 
obtained.  This  measurement  is  important  because  it  indicates  that  electric 
discharge  techniques  can  be  attractive  volumetric,  high  efficiency  initiation 
techniques  for  large  scale  devices.  This  statement  is  supported  by  pin 
discharge  initiated  laser  experiments  done  under  this  Contract  and  else¬ 
where  that  have  achieved  high  electric  efficiency  (-  100%)  in  small  volumes. 

The  future  developments  in  this  area  seem  clear.  In  order  to 
obtain  the  high  power  levels  necessary  for  eventual  applications,  it  will  be 
necessary  to  scale  to  larger  volumes.  The  inherent  high  gain  of  this  type 
of  laser  will  make  the  extraction  of  high  powers  from  large  devices  a  non¬ 
trivial  task.  In  order  to  achieve  high  efficiency  it  will  be  necessary  to  find 
a  volumetric,  scalable  discharge  scheme  for  initiation  of  the  laser.  Finally, 
since  this  device  operates  on  a  pulsed  basis,  for  high  power  applications, 
it  will  be  necessary  to  repetitively  fire  the  laser.  This  requirement  will 
necessitate  quenching  the  chemical  reaction  in  the  cavity  volume,  and 
refilling  with  fresh  gas  without  preignition.  In  addition,  the  acoustic  waves 
generated  from  the  initation  pulse  and  reaction  must  be  damped  for  reas¬ 
onable  medium  homogeneity.  A  recent  RFP  has  been  issued  by  AFWL  to 
address  some  of  these  considerations. 

In  conclusion  it  is  our  opinion  that  a  new  area  of  chemical  laser 
development  is  beginning,  and  the  work  and  achievements  reported  herein 
have,  in  large  measure,  contributed  to  that  birth. 
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SECTION  II 


PROPAGATION  TASK 


2.  1  PULSE  PROPAGATION 
2.  1.  1  Introduction 


Evaluation  of  the  propagation  characteristics  of  a  laser  beam 
through  an  absorbing  fluid  is  a  problem  in  non-linear  optics.  The  inter¬ 
action  of  the  laser  with  the  medium  is  through  the  index  of  refraction  which 
has  both  a  real  and  imaginary  part  giving  rise  to  absorption  and  phase 
changes  associated  with  the  density  changes  of  the  atmosphere.  The  inter¬ 
action  is  non-linear  because  the  index  of  refraction  depends  upon  the 
intensity  of  the  propagating  wave  through  the  hydrodynamic  equations.  The 
non-linear  interaction  process  which  has  been  treated  in  this  investigation 
is  absorption  of  a  pulsed  laser  beam  in  a  medium  in  which  the  pulse  time  is 
of  the  order  of  the  hydrodynamic  time,  i.  e.  ,  the  time  required  for  a  signal 
traveling  at  the  sound  velocity  to  cross  the  beam. 

/  ?  /  \ 

Earlier  work'  '  has  been  devoted  to  a  study  of  the  steady-state 
situation  which  is  obtained  when  a  transverse  wind  is  present  to  remove  the 
absorbed  energy.  This  present  work  is  an  extension  to  the  time  regime 
before  the  wind  has  an  effect.  The  next  step  in  the  development  is  the 
inclusion  of  wind  in  the  time  dependent  problem  in  order  to  study  the  inter¬ 
mediate  times  between  beam  turn  on  and  the  steady-st  ..te  regime. 

This  study  is  restricted  to  collimated  Gaussian  beams  in  the 
limit  where  linearized  hydrodynamics  applies.  Comparison  will  be  made 
with  analytical  theory  for  situations  where  the  theory  is  valid.  In  addition, 
results  for  cases  for  which  the  analytic  theory  is  known  to  be  a  poor  approxi¬ 
mation  will  be  presented. 

2.  1.  2  Interaction  of  a  Pulsed  Laser  Beam  with  an  Absorbing  Atmosphere 

In  this  section  the  dependence  of  the  gas  density  on  the  ir radiance 
distribution  of  the  laser  beam  is  calculated.  The  density  response  is  deter¬ 
mined  from  the  linearized  equations  of  mass,  momentum,  and  energy. 
Viscosity  and  heat  conduction  are  neglected  and  times  of  interest  are 
restricted  to  be  less  than  the  time  required  for  free  convection  to  set  in. 

This  latter  restriction  allows  the  problem  to  be  solved  in  cylindrical 
symmetry,  giving  the  following  set  of  equations  for  the  changes  in  the 
hydrodynamic  variables, 


(109) 


134 


v  V 

y  t 

-  Vp 

(HO) 

d  h 

d  p 

(in) 

TT  ' 

Tt  =  Ql 

£ 

P 

(112) 

where  p,  p,  v,  a,  I,  h  are,  respectively,  the  pressure,  density,  and 
velocity  changes,  total  absorption  coefficient,  Learn  intensity  and  enthalpy, 
and  where  subscript  zero  refers  to  amuient  values. 

The  above  set  of  equations  can  be  reduced  to  the  following  single 
equation  involving  only  the  density: 


4^  Mv  -  1)  Q  v2  i 


(113) 


where  cg  yPo/Po  the  s°und  velocity.  Equation  (113)  is  subject  to  the 
following  initial  conditions 


p  =  0,  3p/9t  =  0,  d2p/dtZ  =  0 


(114) 


which  are  also  the  initial  conditions  of  Eqs.  (109),  (1’'  and  (111), 
respectively. 


For  a  pulse  length  tp  which  is  much  shorter  than  the  hydrodynamic 
time  tpj  (tpj  =  Rm/cs  where  Rm  is  the  laser  beam  radius)  the  density  de¬ 
pendence  is  given  by, 


p  =(  y  - 1 )  Q  1 


(115) 


For  tp  much  longer  than  tp^,  the  density  dependence  is  given  bv 
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„  (Y-l)a  /  T 

P  "  2  /  dt  •  (116) 
Cs  o 

It  the  true  situation  is  modeled  by  assuming  that  the  heating  of  the  gas  is 
done  by  a  non-varying  beam  having  a  Gaussian  shape,  Eqs.  (115)  and  (116) 


p  =  (Y-l)  at3  V2  Iq/6 


for  the  short  time  behavior  and 

(Y-l)  Qtl 

P1  = - 2 -  (118) 

c 


for  the  long  time  behavior  where  I0  denotes  the  Gaussian  beam  intensity  in 
a  vacuo.  It  is  clear  that  for  short  enough  times  the  density  will  always  change 
according  to  Eq.  (117)  for  beams  oC  arbitrarily  high  intensity,  but  it  is 
only  for  small  values  of  a  In  tp  that  Eq.  (117)  will  be  obeyed  out  to  times 
as  long  as  tH.  '  If  the  density  as  given  by  Eqs.  (  1  17)  and  (118)  are  equated 
and  sowed  for  t  one  can  thereby  define  an  intermediate  time  tc  which 
separates  the  two  time  regimes  for  those  cases  where  the  beam  is  not 
highly  distorted  throughout  the  early  regime.  This  time  is  tc  =  n/I^T 
Thus,  for  those  cases  which  result  only  in  a  slightly  modified  profile  at 
times  greater  than  tc,  the  b*»am  follows  Eq.  (117)  for  times  close  to 
tn  and  follows  Eq.  (1  18)  for  times  greater  than  about  3  tH.  In  the  inter¬ 
mediate  regime  at  about  2  tH  the  detailed  behavior  of  the  beam  and  the  fluid 
is  best  studied  self-consistently  by  numerical  procedures.  Furthermore, 
interesting  cases  will  include  those  for  which  the  beam  is  highly  distorted 
from  its  initial  Gaussian  profile  for  t  «tn-  For  these  situations  self- 
consistent  numerical  techniques  are  necessary. 

In  the  optical  frequency  range,  Maxwell's  equations  can  be 
approximated  by  a  parabolic  equation.  The  electric  field  is  taken  to  be  of 
the  form 


E  =  A  (r,  z,  t)  exp  (iknQz  -  az/2)  k  =  2  tt/A  (1  19) 


and  the  index  of  refraction,  n,  is  related  to  the  gas  density  changes  by 
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(120) 
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no  +  ~~  (P  -  PQ) 
Ko 


The  diffracted  amplitude  A  is  complex  and  the  irradiance  I  is  defined  to 
be  AA"\  The  amplitude  obeys  the  following  paraxial  approximation  to  the 
scalar  wave  equation(37) 


2ik 


9  A 


r 


a 

37 


/  2  2  -QZ  . 
(n  -  nQ  )e  A 


=  0 


(121) 


The  derivatives  in  time  (2  ik/c)  3 A/at,  in  Eq.  (121)  are  neglected  since 
the  time  dependence  of  the  index  of  refraction  is  on  the  scale  of  t j_j. 

Since  the  initial  beam  is  assumed  to  be  collimated  at  the  laser 
face,  the  appropriate  boundary  conditions  for  the  diffracted  amplitude  are: 


2 

A  (r,  0,  t)  ^  e  " r 
A  (r,  z,  t)  -  0  as  r  —  °o 


(122) 


A  uniform  square  pulse  in  time  is  studied. 

2.  1.  3  Numerical  Results 

Equation  (121)  has  been  programmed  using  the  stable,  explicit, 
numerical  technique  of  Ref.  (36).  For  r  /  0,  derivatives  in  the  r-direction 
are  replaced  by  the  central  difference  approximation 


+  7  TT~  =  ^7  |  Ai  +  1(  k’ s  ^  +  Trj 


-  2  A 


l,  k,  s 


+  Ai  -  1,  k,  s  l1  "  ir 


■'! 


(123) 


where  Ax  ^  s  =  A  (r,z,  t)  and  r  =  iAr,  z  =  kAz,  and  t  -  sAt.  At  r  =  0, 
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OA/dr  0  and  the  central  difference  approximation  becomes 


l 

r 


9A  -  4  ,  (A 

TT  ~^T  (Al,k,  s 


‘  A0,  k,  s  ^ 


(124) 


The  derivative  in  the  propagation  direction  z  is  replaced  by  the  symmetric 
difference  quotient 


9  A  (Ai,  k  +  1,  s  ~  Ai,  k  -  1,  s  ^ 
S  z  2  Az 


(125) 


where  Az  is  the  width  of  the  mesh  in  the  propagation  direction.  This  dif¬ 
ferencing  scheme  for  Eq.  (12  1),  is  stable  under  the  condition 


|  i  +  _^!_  k2  e-“(n2  -  no2)  j  £  1.0  (126) 


The  initial  value  hydrodynamic  problem  can  be  completely 
formulated  in  terms  ot  a  Green's  function  but  computationally  it  is  far 
simpler  to  use  the  standard  finite  difference  equations  for  a  hyperbolic 
system  of  equations  rather  than  to  complete  the  Green's  function  quadratures 
numerically.  Differencing  of  the  hydrodynamics  wave  equation,  Eq.  (113), 
in  the  spatial  coordinate,  r,  is  identical  to  Eqs,  ( 123)  and  (  124).  Differenc¬ 
ing  in  the  time  direction  is  accomplished  by  central  differences 


dp  1 
dp 


2  At  ^i,  s+1  '  Pi,  s-  1 


f  =-7-2  'Pi,  S+l-2Pi,s+Pi,  .-!> 


(127) 


3t  (At) 


where  plt  s  =  p 
the  condition  that 


(r,t)  and  r  =  iAr,  t  =  sAt.  This  differencing  is  stable^)  under 


c  At 
s 
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Both  the  set  of  equations,  Eqs.  (109),  (110),  (111),  (112),  and 
the  resultant  Eq.  (113)  have  been  programmed  successfully  using  central 
differences  in  time,  and  give  identical  results.  If  forward  differencing  in 
time  is  used,  the  equations  are  unstable  unless  the  concept  of  artificial 
viscosity  is  introduced. 

There  are  two  methods  for  determining  the  spatial  and  time 
variation  of  the  irradiance  distribution.  The  first  method  integrates  in 
propagation  distance  z  at  a  fixed  time  arid  repeats  this  process  until  times 
of  interest  have  evolved.  The  second  numerical  procedure  computes  the 
index  change  as  a  function  of  time  and  space  for  each  propagation  level 
z  =  kAz  starting  with  the  known  initial  profile  at  z  =  0.  The  governing 
equation  for  the  electromagnetic  field  is  computed  for  all  time  and  space 
at  z  =  (k  +  l)Az  using  the  previous  computed  values  Aj  s.  The  procedure 
is  then  repeated  for  all  propagation  levels  and  a  complete  spatial  and  time 
history  is  obtained. 

This  latter  method,  which  was  chosen  for  this  work,  is  used  to 
advantage  when  fewer  steps  in  time  are  required  than  steps  in  range.  This 
can  be  seen  by  reference  to  Table  X  where  the  storage  locations  required 
for  each  method  are  enumerated. 

The  table  assumes  that  central  differencing  in  both  range  and 
time  is  used.  This  explains  the  entries  of  the  number  3  in  Table  X. 

The  two  methods  are  seen  to  require  comparable  storage  when 
nz  =  nt.  This  condition  is  met  when  studying  pulses  of  a  few  hydrodynamic 
times  at  ranges  which  are  appreciable  fractions  of  the  far  field  range.  For 
times  much  longer  than  this  the  inclusion  of  wind  becomes  necessary. 

Wind  introduces  an  asymmetry  which  effectively  replaces  nr  by  nr^/2  in 
the  above  enumeration  making  storage  requirements  prohibitively  large. 

2.1.4  Results 


Figure  55  shows  the  time  evolution  of  the  gas  density  perturba¬ 
tion  for  a  case  where  the  laser  beam  is  never  highly  distorted  throughout 
the  pulse.  A  compressive  wave  propagates  away  from  the  beam  and  its 
strength  decreases  as  r"v^.  The  dashed  curves  in  Fig.  55  represent  the 
approximate  analytic  relations  for  the  short  and  long  time  limit.  The  early 
time  solution  is  valid  for  times  shorter  than  tn  and  the  long  time  limit  is 
valid  for  times  longer  than  3tH»  as  discussed  above. 

A  series  of  plots  of  normalized  intensity  vs  normalized  range  at 
five  different  times  after  pulse  initiation  for  four  ranges  is  shown  in  Fig.  56 
Times  are  given  in  units  of  tH-  Ranges  are  given  in  units  of  the  Rayleigh 
range  kR2M.  As  expected,  the  central  irradiance  decreases  with  time  and 

increases  to  a  peak  at  the  edge  of  the  beam.  This  enhancement  can  be 
understood  from  a  geometrical  optics  argument  which  predicts  that  the  rays 
near  the  center  of  the  beam  will  deflect  outward  to  the  region  of  higher 
index. 
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TABLE  X 


ENUMERATION  OF  STORAGE  REQUIREMENTS 
FOR  ALTERNATIVE  COMPUTING  METHODS 


Variable 


Method  I 

(march  in  z, 
fixed  t) 


Method  II 

(march  in  t, 
fixed  z) 
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number  of  samples  in  r,  z,  t  respectively). 


140 


Rm  =  MIRROR  RADIUS 
Qm  =  SPEED  OF  SOUND 


Time  Dependence  of  the  Normalized  Density  Distribution 
f  ,r  a  Gaussian  Irradiance 


Fig.  56 


Time  Dependence  of  the  Irradiance  Distribution  for 
q  Pr.  =  1.0  and  =  .  133  R„/a 


A  comparison  of  two  very  different  cases  at  the  same  Rayleigh 
range  is  shown  in  Fig.  57.  The  upper  plot  is  taken  from  Fig.  5b  while  the 
lower  plot  depicts  a  case  whcre  the  power  and  beam  size  were  chosen  to 
increase  the  average  irradiance  by  a  factor  of  4.  The  power  was 
decreased  by  a  factor  of  100  and  the  beam  radius  was  decreased  by  a 
factor  of  20.  Thus,  for  the  same  elapsed  real  time  the  smaller  beam 
pulse  is  twenty  times  longer  measured  in  units  of  tH-  In  spite  of  this  higher 
power  per  unit  area  and  longer  time  evolution  with  respect  to  tj-j  the  lower 
beam  is  clearly  less  distorted.  The  difference  can  be  attributed  to  the^ 
large  difference  in  energy  which  has  been  deposited  in  the  volume  7rRm  /a. 
This  energy  is  responsible  for  the  heating  of  the  gas  and  the  changes  in  the 
index  of  refraction.  For  the  case  considered  the  upper  beam  received 
100  times  as  much  energy  as  the  lower  beam  at  the  end  of  their  respective 
pulses,  and  2000  times  the  energy  of  the  lower  beam  at  .  13  of  the 
respective  hydrodynamic  times.  The  energy  delivered  per  unit  volume 
per  un'c  (hydrodynamic)  time  was  in  the  ratio  of  5  to  1,  upper  beam  to  lower 
beam,  using  original  beam  volumes  for  the  calculation. 

The  behavior  of  the  time  averaged  irradiance  is  plotted  in 
Fig.  58  for  three  cases  where  tp  was  varied  but  the  product  atpp  was  held 
constant.  The  time  averaged  results  for  the  longest  pulse  duration 
corresponds  to  a  loss  of  30  percent  of  the  average  far  field  brightness 
(after  correction  for  absorption  effects).  Up  to  .  07  tpj,  the  time  averaged 
irradiance  effectively  defines  a  beam  whose  radius  is  only  slightly  larger 
than  the  undistorted  average  radius.  The  .  1  3  tH  curve  corresponds  to  the 
profiles  shown  in  Fig.  56. 

2.  1.  5  Conclusions 


The  nonlinear  problem  of  laser  propagation  through  an  absorbing 
gas  has  been  studied  using  a  stable,  convergent  nume  rical  algorithm.  The 
study  was  restricted  to  collimated  Gaussian  beams  in  the  limit  where 
linearized  hydrodynamics  applies.  The  computer  results  agree  in  detail 
for  those  cases  where  both  long  and  short  pulse  analytical  theory  is  known 
to  be  valid.  Results  for  situations  for  which  the  analytic  theory  is  known 
to  be  poor  approximation  have  been  presented  (i.e.,  tp  ~  Rm/am).  The 
distortion  of  the  beam  is  a  strong  function  of  the  energy  deposited  in  a  time 
of  the  order  of  the  hydrodynamic  time.  Irradiance  distributions  which  have 
been  averaged  over  a  pulse  length  indicate  that  even  for  beams  which  are 
highly  distorted  by  the  end  of  the  pulse,  the  average  degradation  of  pulse 
shape  can  be  relatively  small. 
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2.2  PRELIMINARY  STUDY  ON  PROPAGA'l  ION  CHARACTERISTICS 
OF  UNSTABLE  RESONATORS 

2.  2.  1  Introduction 

Unstable  resonators  at  high  output  coupling  are  an  efficient 
means  of  extracting  high  power  in  lasers.  Thus  a  study  on  their  propaga¬ 
tion  characteristics  was  undertaken.  This  study  is  preliminary  in  nature 
and  was  undertaken  using  present  finite  difference  methods.  The  problem 
with  sharply  truncated  beams  is  that  they  introduce  infinite  second  deriva¬ 
tive  at  the  edges  into  the  analysis.  It  is  not  clear  on  how  appropriate  or 
accurate  these  schemes  are  in  treating  truncated  beams.  Nonetheless, 
we  proceeded  to  make  an  analysis  using  the  stable  explicit  finite  differenc¬ 
ing  scheme.  However,  future  work  in  this  area  will  require  a  totally  new 
numerical  approach  to  resolve  the  computational  uncertainties. 

Unstable  resonators  have  a  relatively  flat  irradiance  (Fig.  59) 
distribution  and  this  irradiance  distribution  should  reduce  somewhat  the 
blooming  transverse  to  the  wind.  No  matter  how  flat  the  irradiance 
distribution  is  in  circular  unstable  resonators  these  resonators  have  an 
astigmatic  component  due  to  the  combination  of  a  circular  irradiance 
distribution  and  a  flow.  The  density  distribution  at  z  =  0  associated  with 
a  circular  unstable  resonator  is 


x 


The  linear  term  in  x  is  simply  a  ray  deflection  term  which  does  not 
introduce  any  significant  defocusing  effect.  However  the  y-dependence 
introduces  an  astigmatic  component  and  we  can  expect  crescents. 

If  we  have  a  square  or  rectangular  unstable  resonator  at  high 
output  coupling  the  density  distribution  is  simply  linear  without  a  gradient 
in  y.  At  first  glance  in  the  near  field  this  should  introduce  just  a  bending 
term  without  much  irradiance  degradation  in  the  near  field.  This,  how¬ 
ever,  is  not  borne  out  in  the  detailed  numerical  calculation.  In  fact,  the 
diffraction  effects  of  the  edges  introduce  such  peaking  of  the  irradiance 
distribution  that  rectangular  and/or  square  unstable  resonators  have 
density  distributions  with  more  severe  gradients  than  either  circular 
unstable  resonators  or  Gaussian  beams. 

2.  2.  2  Theoretical  Analysis 

In  the  optical  range,  the  circular  frequency  of  oscillation  of 
the  electromagnetic  field  is  very  great  and  Maxwell's  equation  admits  an 
asymptotic  series  expansion  in  powers  of  l/kR^.  The  first  term  in  this 
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asymptotic  expansion  gives  geometric  optics.  Geometric  optics,  by 
replacing  Maxwell's  equations  by  two  first-order  partial  differential 
equations,  is  a  considerable  simplification.  Geometric  optics  gives  in¬ 
finite  irradiance  at  certain  singular  points  or  curves  such  as  foci  and 
caustics;  at  these  singularities,  diffraction  effects  must  be  considered. 

To  include  diffraction  effects,  we  put 


E  =  A  (x,  y,  z)  exp  (iknw  z  +  at  z/2);  k  =  2  ti/\,  (128) 


where  A,  the  diffraction  amplitude,  has  both  a  real  and  imaginary  part. 

The  irradiance  I  (x,  y,  z)  is  defined  to  be  the  A  times  the  complex  conjugate 
A::.  Substitution^  Eq.  (128)  into  Maxwell's  equation  and  neglecting  terms 
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X  A  |  j"  AA  '1  |  1  -  6  exp  ^  -  (x  -  x'  )/p  J  j  dx'j  =  0,  (129) 


where 


€0 


z0  An=o 


(130) 


The  nondimensional  quantity  kRw  /z0  is 


related  *o  the  Fresnel  number. 


'Vo  assume  that  we  have  a  light  beam  whose  rays  are  focused 
at  z  =  zf.  The  appropriate  boundary  conditions  xor  the  diffracted  amplitude 
are 


A  (x,  y,  0)  =  A  i  (x,  y),  (131a) 

A  (x,  y,  z)  =  0  as  x,  y  —  +  «  (13  lb) 

The  governing  equation  has  been  programmed  using  the  stable,  explicit 
numerical  te  chnique  of  Ref. (38).  We  replace  derivatives  in  the  s,  y  direction 
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by  the  central  different 
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where  A  is  the  width  of  the  mesh  perpendicular  of  the  propagation  direction. 
The  derivative  in  the  propagation  direction  z  is  replaced  by  the  symmetric 
difference  quotient 


A 

z 


^Ai,  j,  k+  1  '  Ai,  j,  k-  1  ^ 
2  A  z 


(133) 


where  A  z  is  the  width  of  the  mesh  in  the  propagation  direction.  This 
differencing  scheme  is  stable  under  the  condition 


4e  A  z 

72 


<1.0 


(134) 


The  second  derivative  in  Eq.  (128)  at  the  edge  is  very  large  and 
gives  rise  to  a  solution  which  oscillates  extremely  fast  near  z  =  0.  This 
requires  an  extremely  small  grid  structure  to  accurately  calculate  the 
complex  amplitude.  Thus,  calculations  with  relatively  small  Fresnel 
number  are  the  only  ones  which  can  be  computed  with  any  degree  of 
confidence. 


2.  2.  3  Discussion  and  Results 


We  now  present  calculations  for  a  circular  unstable  resonator 
focused  at  1  km  and  compare  with  an  infinite  Gaussian  beam  also  focused 
at  1  km.  The  unstable  resonator  at  the  focal  point  z  =  1.  0  has  a  far  field 
brightness  reduced  by  a  factor  of  2  over  a  Gaussian  beam.  The  peak  ir- 
radiance  is  a  maximum  for  both  beams  at  z  =  .  50f  (minimum  waste).  The 
phase  changes  associated  with  the  heating  simply  does  not  allow  the 
beam  to  be  focused  any  smaller.  The  physical  reason  for  the  unstable 
resonator  having  less  irradiance  at  the  focal  point  is  due  to  the  fact  that 
diffraction  for  a  truncated  beam  is  much  greater  as  we  approach  z  =  l. 

This  reduces  the  far  field  brightness  over  the  infinite  Gaussian  beam  as 
shown  in  Fig.  60. 


Summarizing,  for  those  calculations  which  are  accessible 
(relatively  low  Fresnel  number)  using  the  present  explicit  finite  difference 
scheme,  unstable  resonators  do  not  propagate  as  well  as  infinite  Gaussian 
beams.  Since  neither  shape  is  acceptable  as  far  as  far  field  brightness 
is  concerned,  new  emphasis  must  be  put  on  beam  shaping  which  minimize 
these  effects.  Also,  further  computational  procedures  must  be  developed 
in  order  to  calculate  propagation  characteristics  of  beams  with  high  Fresnel 
numbers.  High  Fresnel  number  beams  (i.  e.  large  radius)  are  one  other 
way  to  minimize  these  refractive  effects. 
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APPENDIX  I 


THE  METHOD  USED  FOR  NUMERICAL  EVALUATION 

00 

OF  THE  INTEGRAL  J  R£  (kr)  R{  (k'r)  V  (r)  dr 

0 


In  this  appendix  we  describe  a  method  for  obtaining  approximate 
solutions  Rji  to  the  equation 


d  Rf  (kr) 


dr 

where  we  take 


k2  _  l  {i  t.lj.  _  2ji  v  (r) 
r2  K 


Rf  (kr)  =  0  (135) 


V  (r)  =  D 


[e‘2Mr'R0)_2e'P(r'R°)] 


(136) 


i.  e.  ,  assume  a  Morse  type  potential  with  well  depth  D  at  r  =  Rq  and  range 
parameter  The  solutions  of  inte  est  to  us  satisfy  two  boundary  conditions 


lim  Rj  (kr)  =  0 


(137) 


and 

R^  (kr)  ~  cos  (kr  +  6)  for  r>  r^ 
where  r^  is  determined  by  the  condition  that 


(138) 
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We  also  point  out  that  we  are  interested  in  solution  of  Eq.  (135)  to  be  able 
to  evaluate  ti  e  integrals  of  the  type 
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J  (kr)  R^  (k'r)  f  (r)  dr  (140) 

0 

where  f  ( r)  is  a  smooth  function  of  r  which  goes  to  zero  at  large  r.  If 
k  <  k*  and  R^  and  R^  are  defined  such  that 


(141) 


(142) 


the  integral  (140)  can  be  approximated  as 


R£  (kr)  Rfi  (k'r)  f  (r)  dr  (143) 


where  is  the  number  of  the  order  of  1  A  while  rnay  equal  several  ang¬ 

stroms  depending  upon  k,  k'  and  f  (r)  The  reason  is  that  for  r  <  Rj  *  e  j  > 

Rg  (kr)  and  R^  (k1  r)  are  fast  decaying  functions  while  for  r  >  R£  +  t^le 
product  Rjj,  (kr)  R^  (k'r)  oscillates  approximately  as  cos  ((L  -  k')  r).  In  the 
latter  case  the  Expression  (143)  is  good  approximation  to  integral  (135)  if  (i) 
f  (R£  +  €3)  <<.  f  (R^),  i.e.,  iff  (r)  goes  to  zero  rapidly  enough  with  increasing 
r  or  (ii)  IK-  k’  I  r  >>  1  so  that  the  product  R^  (kr)  R*  (k1  r)  oscillates  rapidly 
and  |  f  (r)  I  is  monotonically  decreasing  function  of  r  for  r  >R2  +  ^2'  this 
appendix  we  restrict  ourselves  to  the  case  £  =  £'  and  compare  our  approxi¬ 
mate  results  with  those  obtained  by  numerical  integration  of  Eqs.  (135)  and 
(141)  for  the  case  that  f  (r)  =  V  (r). 

To  obtain  the  approximate  solution  of  Eq.  (135)  we  write 


k  -  2.46  x  1017  e’2pr  +  4.  74  x  1016  e'^r  +  4.  33  x  1  01 4  (144) 

r 


-154- 


a: 


Figure  61  shows  a  graph  of  right-hand  side  of  Eq.  (144)  and  shows  that  the 
above  approximation  is  a  reasonable  one  for  2  x  10'^  <  r  <  5  x  10"^  =  cm. 
Now  we  define 


?  n  7 

a„  =  e  -  2.46  x  lo'  £  (£  +  1) 

c  K2 


(145) 


p.  =  e  0  -  4.  74  x  1016  £  (£  +  1) 

£  X 


(146) 


k  =  k2  -  4.  33  x  1014  £  (£  +  1) 


(147) 


We  will  restrict  our  approximate  method  to  the  case  for  which 

a  f  —  0- 


(148) 


Equation  (135)  can  now  be  written  as 


kc  +  “t  e"2Pr  +  pt  e' 


>r  R,  =  o 


(149) 


Restriction  (148)  makes  V  (r)«  e'^r  for  r«0  and  permits  us  to  alte: 
the  boundary  condition  (137)  to 

lim  (kr)  =  0. 

r  — -  oo 


(150) 


The  solution  of  (149)  which  satisfies  boundary  conditions  (148)  and  (150)  is 
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and 

.  =  |  e'* 

(153) 

r  in  Eq.  (150)  is  gamma  function  and  W  is  a  Whittaker  function. 


The  appropriate  form  of  integral  is 
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Equation  (154)  gives  the  approximate  value  of  the  integral  when 
the  classical  turning  points  (Eqs.  (141)  and  142))  are  less  than  5  R.  For 
values  of  f  and  k  such  that  one  of  the  two  turning  points  is  greater  than  5  R, 
we  approximate  the  solution  of  (135)  by  Bessel  functions,  i.e.  ,  we  ignore 
V  (r)  term  in  £q.  (135).  In  this  approximation 

(kr)  =  kr  (kr)  (16 

where  1  s  are  spherical  Bessel  functions.  The  integral  to  be  evaluated 
now  is 
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where  Q/>  is  Leeendre  function  of  the  second  kind,  Q'  ^  (z)  =  d/dz  (z)  and 
z  =  P^  +  k^  +  k'^/Zkk1.  Using  the  recursion  relation 


Q'  f  (z)  =  (z2  -  l)"1  (f  +  1)  [Q|  +  j  (z)  -  zQ(  (z)j 


one  can  write 
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where  z^  =  4P2  +  k2  +  k,t/2kk'. 


Using  Eqs.  (165)  and  (166),  we  can  write 
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Table  XI  compares  the  results  obtained  for  the  left-hand  si  Je  of 
Eq.  (154)  by  numerical  integration  (Col.  a)  with  the  approximate  expression 
given  on  the  right-hand  side  of  Eq.  (154)  (Col.  b).  Column  c  is  that  obtained 
by  Takayanagi's  modified  wave  number  approximation  while  Col.  d  is  ob¬ 
tained  using  Eq.  (167).  The  potential  parameters  used  are  (3  =  1.4  x  10®  cm"*, 
R0  =  4  x  10'°  cm  and  D  =  5  x  10' 14  ergs.  The  potential  parameters  are  con¬ 
sidered  appropriate  for  H-Ar  interaction.  The  results  show  that  the  numeri¬ 
cal  results  and  the  approximate  expression  on  the  L.  H.  S  of  Eq.  (154)  agree 
well  before  the  former  changes  sign.  After  the  change  of  sign  the  numerical 
integral  agrees  well  with  that  obtained  using  Eq.  (167).  The  physical  basis  of 
this  result  is  that  a  change  in  sign  of  the  integral  occurs  when  the  attractive 
term  in  the  Morse  potential  becomes  more  important.  At  that  point  the  tra¬ 
jectory  is  hardly  deflected  and  one  can  approximate  Rp  (kr)  by  kr  jp  (kr). 
Operationally  during  a  calculation  we  switched  ;vev  from  Col.  b  to  Col.  d 
when  Col.  d  changed  sign.  Takayanagi's  results  agree  well  with  the  numeri¬ 
cally  obtained  ones  when  kssk'.  When  k  /  k' ,  the  Takayanagi  modified 
wave  number  approximation  appears  to  be  valid  for  low  values  of  f . 

Table  XII  demonstrates  the  smoothness  of  the  transition  when 
change  from  Col.  b  to  d  is  made  when  d  changes  sign. 

For  f  =49  and  50  the  sign  of  the  results  in  Cols,  b  and  d  differs 
from  that  in  Col.  a.  This  disagreement  should  not  be  taken  seriously.  The 
absolute  val  -e  of  the  integral  for  these  partial  waves  is  about  an  order  of 
magnitude  smaller  than  the  nearby  ones  with  the  result  that  their  contribu¬ 
tion  to  the  cross  section  is  negligible. 
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TABLE  XII 

COMPARISON  ON  THE  NUMERICAL  VALUES  OBTAINED  BY 
METHODS  OF  EXACT  INTEGRATION,  BESSEL  APPROXIMATION  AND 

THE  PRESENT  THEORY 


k  =  I.  35  x  109  k'  =  1.  43  x  109 


i 

a 

b 

c 

40 

3.  65  x  10‘22 

3.  20  x  IQ’22 

2.  84  x  10"21 

42 

2.  75  x  10‘22 

2.  47  x  10"22 

1 .  66  x  IQ'21 

44 

1.  88  x  10"22 

1. 79  x  10'22 

9.  16  x  10‘22 

46 

1.  06  x  10*22 

1.  18  x  ID’22 

4.  53  x  10"22 

48 

3.  26  x  10'23 

6.  54  x  IQ'23 

1.  73  x  10"22 

49 

-1.  85  x  10-25 

4.  29  x  10‘23 

8.  06  x  10‘23 

50 

-2.  91  x  10-23 

2.  30  x  10-23 

1.  10  x  10"23 

51 

-5.  59  x  10'23 

5.  89  x  10'24 

-4.  03  x  10"23 

52 

-7.  83  x  10-23 

-8.  31  x  10’24 

-7.  73  x  10"23 

54 

-1.11  x  10-22 

-2.  76  x  10'23 

-1.  20  x  10"22 

56 

-  1.  29  x  10-22 

-3.  53  x  10'23 

-1.  35  x  10'22 

58 

-1.  34  x  10'22 

-3.  22  x  10-23 

-1.  34  x  10'22 

60 

-1.  28  x  10"22 

-2.  15  x  10'23 

-1.  25  x  10'22 

(a)  Numerical  integration 

(b)  Present  results 

(c)  Bessel  approximation 
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THE  DERIVATION  FOR  TrlE  EXPRESSION  OF  <xB  e 


-lS  •  I, 


-\9  •  Iv 

In  this  appendix  we  evaluate  the  matrix  elements  <Xf  e 
where  Xi  an£l  Xf  are  wave  functions  for  the  rotational  and  vibrational  motion, 
cj  is  the  momentum  exchanged  during  the  collision  and  _ry  is  the  distance  of 
the  vth  atom  of  the  molecule  from  the  center  of  mass  of  the  molecule.  We 
consider  only  diatomic  molecules.  In  this  case  _r^  =  rr^/m,  +  m2  r  and 
£2  =  _mj/rnj  +  m2  £>  where  m.  and  m2  are  the  masses  of  the  atom  1  and 
2,  respectively,  and  r  is  the  internuclear  distance.  Further  we  write 
£  =  _b  +  u,  where  b  is  the  equilibrium  internuclear  dist?»nce  and  u  is  the 
vibrational  displacement.  We  also  write  x^  =  (u)  x  Yjimi  (ft)'  whe  re  is 

harmonic-oscillator  wave  function  for  vibrational  level  i  and  are  the 

J  1 1* 

spherical  harmonics. 

I  -is  .  £,  I 


-is  •  li 

<Y.  <  <Af  e  <A.>  Y. 

jfm{  vf  j.m.> 


.  m2 


-i  T-f-  q  .  u  -i  -t-r-  q  .  b 

<Yjm  <*t  e  k>«  >  (168) 


l.m. 

1 


<<t> f\  e 


” 1  TT"  -S-u 


6.  >  =  47 T 

1 


(-)  1  Y„  (fi  )  Y*  (SI  ) 
£  ,m,  9  9  ,m  '  V 


/  **  h ,  (tt  “ u)  *1 


(169) 


In  Eq.  (169)  we  have  written  Y (fi,)  f°r  Y  .  (fi  ) 

Jfj  tn j  b  jfj  m?  '  u 

because  band  u  are  collinear,  j.  are  spherical  Bessel  functions  and 
M  =  m  j  +  m^. 
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The  vibrational  wave  functions  have  the  property 


<t>v  (-u)  =  (-)V  4>v  (u)  (170) 

Since  the  spherical  Bessel  Functions  have  the  same  parity  as  the  order  of 
the  Bessel  Function,  i.e., 


f  i 

).  (-*)  =  (-)  j  (x), 

1  1  1  1 


(171) 


the  integral  over  the  wave  functions  vanishes  unless  i  -I  f  +  p  is  an  even 
number.  We  can  therefore  write 


(a jU)  exp  (-u  )  (u)  X  Hf 


(u) 


X 


(174) 


Using  the  relation 


[-f] 

Hn  <x>  =  Z2 
m  =  0 


,  .m  . n-  2m 

(')  (2x> 

m !  (n  -  2  m)  ! 


(176 


where  [  n/2]  is  the  largest  integer  less  than  or  equal  to  n/2,  Eq.  (174)  can 
be  rewritten  as 


X 


m!  (i-  2m)  !  (n!)  (f  -  2n)  ! 


j  (a1u)e-u2  (u,l  +  *-2m-2n 

0 


where  R1  (/  j ,  i,  f,  q) 


where  jF 
function. 


<(f>[  e 


1 

=  T 


)  (f!)  2 


i+f 


1/2 


2 

(a;  b;  z)  =  1  +  -^  z  +  ^  -  -j-|-  -|y  +  .  .  .  *s  confluent  hypergeometric 


Substituting  Eqs.  (177),  (178),  and  (172)  in  Eq.  (169)  we  get 


-l 


m_ 

Tr3'"l*i>  =  4'  (V'  Vl(v 

/  1  f  f  +  h  \ 

X  1 1  +  (-)  1  1  R1  (;  r  i,  f, 


q) 


(179) 


Using  the  expansion  for  the  planewave  and  standard  relations  between 
spherical  harmonics  we  can  write 


168 


REFERENCES 


1.  Pimentel,  G.  C.  and  Kasper,  J.  V.  V  ,  "HCC  Chemical  Laser," 

Phys.  Rev.  Letters  14,  352  (1965). 

2.  Batovskiy,  O.  M.  ,  Vasilyev,  G.  K.  ,  Makarov,  Ye.  F.  and  Talroze, 

V.  L.  ,  "Chemical  Laser  Based  on  a  Branched  Chain  Reaction  of 
Fluorine  with  Hydrogen,"  Zhurnal  ekspe rimental'  noy  i 
teoreticheskoy  fiuiki,  Vol  9,  No.  6  p  341-343  (1969). 

3.  Polanyi,  J.  C.  and  Woodall,  K.  B.,  "Energy  Distribution  among 
Reaction  Prodictf  VI.  F+H2,  D2,  "  J.  Chem.  Phys.  57.  1574  (1972); 
Polanyi,  J.  C.  and  Sloan,  J.  J.  ,  "Energy  Distribution  among 
Reaction  Produces  VII.  H+F2,"  J.  Chem.  Phys.  Sl_,  4988  (1972). 

4.  Wilson,  J.  and  Stephenson,  J.  C.  ,  "Atmospheric  Pressure  Pulsed 
Chemical  Laser,"  Applied  Phys.  Letters  20^,  64  (1972). 

5.  Wilson,  J.  ,  et  al,  "Exploratory  Development  Work  in  Support  of  the 
High  Power/Energy  Laser  Program  (U),  "  Final  Report  under 
Contract  F29601  -70-C-0073  prepared  for  Air  Force  Weapons  Lab¬ 
oratory,  Kirtland  Air  Force  Base,  New  Mexico,  July  1971. 

6.  Ksander,  Yuri,  "Soviet  Chemical  Laser  Research,"  ARPA  Report 
189-1  (1971). 

7.  Kapralova,  G.  A.,  Trofimova,  E.  M.  ,  Shilou,  A.  E.,  Kinetika  i 
Kataliz,  6,  977  (1965). 

8.  Levy,  J.  B.  and  Copeland,  B.  K.  K.  ,  "The  Kinetics  of  the  Hydrogen 
Fluorine  Reaction  III  The  Photochemical  Reaction,  "  J.  Phys.  Chem. 
72,  3168  (1968) 

9.  Hancock,  J.  K.  and  Green,  W.  H.  ,  "Laser  Excited  Vibrational 
Relaxation  Studies  of  HF,  "  J.  Chem.  Phys.  56,  2474  (1972). 

10.  Hess,  L.  D.  ,  "HF  Chemical  Laser  Studies,  Use  of  MoF^  to  Increase 
Reaction  Rates  in  H2  -  F2  Mixtures,"  J.  Appl.  Phys.  4  3,  1  157  (1972). 

11a.  Fried,  S.  ,  Wilson,  J.  ,  and  Taylor,  R.,  "Temperature  Dependence 
of  HF  Vibrational  Relaxation,"  presented  at  3rd  Conf.  on  Chemical 
and  Molecular  Laser,  St.  Louis  (1972)  in  press. 

lib.  Polanyi,  J.  C.  ,  "Vibrational  Rotational  Pcpulation  Inversion,"  Appl. 
Opt.  Suppl.  2,  109  (1965). 


169 


n 


12.  Deutsch,  T.  F.  ,  "Laser  Emission  from  HF  Rotational  Transitions," 
Appl.  Phys.  Letters  1  1 ,  18  (1967). 

13.  Cohen,  N.  ,  Tech.  Rept.  TR-0172  (2779)-2  The  Aerospace  Corpora- 
•  tion,  Sept.  1971. 

14.  Airey,  J.  R.  ,  "C£  +  HBr  Pulsed  Chemical  Laser;  A  Theoretical  and 
Experimental  Study,  "  J.  Chem.  Phys.  j>2,  156  (1970). 

15.  Meredith,  R.  E.  ,  Spellieg,  R.  L.  and  Smith,  F.  G.  ,  "Strengths  and 
Collision  Widths  in  the  Second  Overtone  Band  of  Hydrogen  Fluoride,  " 
J.  Chem.  Phys.  57,  51  19  (1972). 

16.  Airey  and  Smith,  "Quenching  of  IR  Chemilumene scence ;  Rates  of 
Energy  Transfer  from  HF  (V  <  5)  to  CO?  and  HF,  and  from  DF  (V  <  3) 
to  C02  and  HF,  "  J.  Chem.  Phys.  57,  1669  (1972). 

17.  Rigord,  W.  W.  ,  "Saturation  Effects  in  High  Gain  Lasers,"  J.  Appl. 
Phys.  36,  2487  (1965). 

18.  Polanyi,  J.  C.  and  Woodall,  K.  B.,  "Mechanism  of  Rotational 
Relaxation,"  J.  Chem.  Phys.  56,  1563  (1972). 

19.  Suchard,  S.  N.  ,  Kerber,  R.  L.  ,  Emanuel,  G.  ,  and  Whittier,  J.  S.  , 
"Effect  of  H2  Pressure  on  Pulsed  H2+F7  Laser,  Experiment  and 
Theory,"  J.  Chem.  Phys.  57,  5065  (1972). 

20.  Krogh,  O.  D.  and  Pimentel,  G.  C.  ,  "Chemical  Lasers  from  the 
Reaction  of  CfF  and  CCF3  with  H2  and  CH4.  A  Possible  Chain 
Branching  Chemical  Laser,"  J.  Chem.  Phys.  56,  969  (1972). 

21.  Gensel,  P.  ,  Kompa,  K.  L.  and  Wammer,  J.,  "IF5  -  Hz  Hvdrogen 
Fluorine  Chemical  Laser  Involving  A  Chain  Reaction,  "  Chem.  Phys. 
Letters  5,  179(1970). 

22.  Gregg,  D.  W.  ,  Krawetz,  B.,  Pearson,  R.  K.  ,  Schleicher,  B.  R.  , 
Thomas,  S  J.  ,  Huss,  E.B.,  Pettipiece,  K.  J.,  Cheighton,  J.  R.  , 
Niver,  R.  E  ,  and  Pan,  Y.  L.  ,  "Electron  Beam  and  Flashlamp 
Initiation  of  a  Pulsed  Hydrogen  Fluoride  Chemical  Laser,  "  Chem. 
Phys.  Letters  8,  609  (1971). 

23.  Wilson,  J.  ,  and  Stephenson,  J.  C.  ,  "Atmospheric  Pressure  Pulsed 
Chemical  Laser,"  App.  Phys.  Letters  Vol  20,  pp  64-66  (1972). 

24.  Airey,  J.  R.  ,  "C£  +  HBr  Pu  sed  Chemical  Laser:  A  Theoretical 
and  Experimental  Study,  "  J.  Chem.  Phys.  Vol  52,  pp  156-167  (1970) 

25.  Airey,  J.  R.  and  McKay,  S.  F.,  "A  Supersonic  Mixing  Chemical 
Laser,"  Appl.  Phys.  Letters  J_5,  401  (1969). 


r 


1 70 


26a,  Wilson,  J.,  Stephenson,  J.  C.,  and  Northam,  D.  ,  "Exploratory 

Development  Work  in  Support  of  the  High  Power /Energy  Laser  Pro¬ 
gram  (U),  "  Final  Report  under  Contract  F29601-70-C-0073  prepared 
for  Air  Force  Weapons  Lab.  and  Advanced  Research  Projects  Agency, 
July  1971. 

26b.  Wilson,  J.  ,  Northam,  D.  and  Lewis,  P.  ,  "HF  Laser  Action  above 
the  Second  Explosion  Limit,  "  paper  presented  at  3rd  Symposium  on 
Chemical  and  Molecular  Lasers,  St.  Louis,  April  1972. 

27.  Landshoff,  R.  K.  and  Magee,  J.  L.  ,  "Thermal  Radiation,"  Lockheed 
Report  3-27-67-1  Vol  4,  Nov.  1967,  DASA  1917-4. 

28.  Mandl,  A.,  AVCO  Everett  Research  Laboratory,  private  communica¬ 
tion. 

29.  Fox,  R.  E.,  "Dissociative  Attachment  of  Electrons  in  I;,  "  Phys. 

Rev.  209,  2008  (1958). 

30.  Kompa,  K.  L.  and  Wanner,  J.  ,  "Chem.  Phys.  Letters  5,  179  (1970). 

31.  Christophorou,  L.  G.  ,  "Atomic  and  Molecular  Radiation  Physics," 
Wiley-Inter science,  New  York  (1971). 

32.  Hancock,  J.  K.  and  Green,  W.  H.  ,  "Las^r  Excited  Vibrational 
Relaxation  Studies  of  HF,  "  J.  Chem.  Phys.  J36,  2472  (1972). 

33.  Wilson,  J.  ,  Chen,  H-L,  Fyfe,  W.  ,  Taylor,  R.  L.,  Little,  R.  and 
Lowell,  R.,  "Electron  Beam  Dissociation  of  Fluorine,  J.  Chem. 

Phys .  to  be  published. 

34.  Berger,  J.  J.  and  Seltzer,  S.  M,  ,  "Tables  of  Energy  Losses  and 
Ranges  of  Electrons  and  Positrons,"  Paper  10,  NAS-NRC  Publication 
1133  (1964),  also  NASA  SP-3012  (1964). 

35.  Bortner,  T.  E.  and  Hurst,  G.  S.,  "Ionization  of  Pure  Gases  and 
Mixtures  of  Gases  by  5-Mev  Alpha  Particles,"  Phy.  Rev.  Vol  93, 
pp  1236-1241  (1954). 

36.  Hayes,  J.  N.  ,  Ulrich,  P.  B.,  and  Aitken,  A.  H.  ,  "Effects  of  the 
Atmosphere  on  the  Propagation  of  10.6/i  Laser  Beams,"  Appl.  Opt. 

U_,  257  (1972). 

37.  Hayes,  J.  N.  ,  "Thermal  Blooming  of  Laser  Beams  in  Fluids,"  Appl. 
Opt.  1J_,  455  (1972). 

38.  Harmuth,  H.  F.,  "On  the  Solution  of  the  Schroedinger  and  the  Klein- 
Gordon  Equations  by  Digital  Computers,"  J.  Math.  Phys.  36,  269 


